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J. Hadamard's ideas about the correct statement of the problems of 
mathematical physics have been analyzed. In this connection various inter- 
pretations of the directly related Banach theorem about the inverse operator 
has been touched. The contemporary apparatus of mathematical modeling 
is shown to be in a drastic contradiction with concepts of J. Hadamard, S. 
Banach and a number of other outstanding scientists in the sense that the 
priority is given to the realization of algorithms, which actually imply that 
incorrectly stated problems are adequate to real phenomena. 

A new method is developed for solving problems traditionally associated 
with the Fredholm integral equation of the first kind = /, X G [0, 1]. 
It is based on the representation of the integration error in the form Sf = 
ip — XBip, X G [0, 1], where B is the integral operator with limits —1,1 and 
Poisson kernel; A is parameter. Incompletely continuous perturbation of 
operator A with I — XB, provided that 5f — 0, makes it possible to change 
the statement of the problem. This involves (i) the extension of the problem 
6f = 0; 6f = jiAip — jif {ji is parameter) onto x E [—1, 0) and (ii) the use 
of equations with similar structure and the same function ^, x G [0, 1]. The 
essence of this is the practical realization of the condition / + G R {A). 
A key point here is to interrelate the components of the above systems of 
equations to enable their mutual conversion. In the case when function 
is harmonic the problem is reduced to a Fredholm integral equation of the 
second kind with properties favorable in computational respect. 

The solution to this equation is given in the form of Fourier series with 
coefficients depending on parameters of particular problem and also param- 
eter < r < 1. The class of possible may be extended to L2 by the 
limit transition r — 1. In the second approach, the belonging ip E L2 was 
assumed from the very beginning. Accordingly, condition 6f = needs to 
be addressed in terms of generalized functions. In comparison with the first 
approach, this one is more formal; relatively simple transformations result 
in second-order Fredholm integral equation with properties most favorable 
for the numerical realization. 

The general concept in this book is as follows. There is one and only one 
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function ip = ip^: such that Aip^: = f , but the problem to restore it from this 
equation with known A and / is incorrect. Along with this, it is not difficult 
to imagine a Predholm integral equation of the second kind with such a free 
term that ■0* — IJ'Zijj^ + F, x E [0, 1]. The essence of this text is to show 
how to construct this equation starting from Alp = / . In other words, the 
problem to find a function that satisfies the Predholm integral equation of 
the first kind is stated correctly. 

The possibility is shown to extend the approach suggested for a wide 
circle of problems that may be reduced to two-dimensional Predholm inte- 
gral equations of the first kind; these are linear boundary-value and initial- 
boundary-value problems with variable coefficients, non-canonical domain 
of definition and other peculiarities complicating their solution. The elabo- 
rated algorithm is shown to be directly applicable to them. Note that this 
may be used for examining the above problems for solvability. 

In discussing the statement of problems of mathematical physics, con- 
siderable attention is paid to methodological aspects. Conclusions about 
cause-and-effect relations are argued to be essentially illegitimate when the 
solution of a problem is traced in the long run to a primitive renaming of 
known and unknown functions of a corresponding direct problem. The aim 
of this work is a constructive realization of J. Hadamard's opinions that 
physically meaningful problems always have correct statements. 

E-Mail: eperchik@bk.ru 
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Chapter 1 
Introduction 



At the beginning, we should explain the title of the work and, in the first 
place, the meaning of the employed notions. In this regard, we assume 
the availability of information allowing us to formulate a mathematical 
model of a certain phenomenon in a traditional way. Correspondingly, 
the determination of unknown functions using the data of the problem 
is implied. If the dependence of the solution to the problem on these 
data with respect to the norm of the chosen space is continuous, such a 
problem, as a rule, belongs to the domain of analysis or, in other words, 
its formulation is direct. 

However, the investigation of a concrete phenomenon in a variety 
of the determining factors with the aim of obtaining, as a final result, 
of qualitatively new information (the synthesis of knowledge) also en- 
visages the realization of problems in their inverse formulation,i.e., the 
restoration of data using the hypothetically known solution: In other 
words, the restoration of the cause using its consequence, which is usu- 
ally identified with the necessity of solving ill-posed problems. 

The purpose of the present investigation consists in the justification 
of the illegitimacy of this statement and, on the contrary, in a construc- 
tive development of J. Hadamard's ideas of the existence of well-posed 
problems, adequately describing real processes and phenomena. Note 

that the difference between these two notions in the context of the book 
is unessential. However, the term "process" accentuates a time factor. 

In the focus of the attention is a natural, to our mind, issue that. 
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as an example, can be explained by the evaluation of the integral 

1 

(A^) {x) = Jk {x, ^ iO d^-fix), xe [0, 1] , (1.1) 



which amounts to the determination of the function / (,t) using given 
k{x,^) and ip (x) (from the space L2). This procedure can be easily 
associated with a lot of physical, as well as other, interpretations. Its 
realization, at least in the case of the bounded integrand, does not pose 
any problem. 

On the other hand, if the kernel k{x,^), which is assumed to be 
bounded, and the function / (x), evaluated beforehand from Eq. (1.1), 
are given, the function ip {x) is objectively existent and unique. Thus, 
the questrion is whether it is legitimate to restore this function by 
means of the solution of the Fredholm integral equation of the first 
kind (1.1), just renaming the known and unknown components in the 
formulation of the direct problem, i.e., by assuming that the function 
/ {x) is given and ip {x) is to be determined. And, generally speaking, 
what is the basis to argue that mathematical formulations of the direct 
the inverse problems can be absolutely identical? 

The very fact of mechanical renaming of the known and unknown 
functions, without any additional corrections, raises objections. Thus, 
we put forward a thesis that an adequate approach to the formulation of 
inverse problems should differ from the approach that became common. 
This position predetermined the presence in the title of the work the 
notion of methodology. 

As a matter of fact, we hope to find reserves of the synthesis of the 
whole complex of knowledge about the phenomenon by investigating it 
from different sides, using formulations whose mathematical represen- 
tations arc not identical. Although Eq. (1.1) is absohitcly sufficient for 
the evaluation of / (x), the restoration of the function ip (x) does not 
necessarily consists in the solution of the Fredholm integral equation of 
the first kind, which is an ill-posed problem. 

However, if there exists an alternative to the above-mentioned re- 
naming of the known and unknown components, one can assume that 
corresponding formulations of the inverse problems may possess much 
more attractive properties in a computational sense. From this point of 



9 



view, the arguments of Hadamard acquire a rather concrete meaning, 
stimulating a search of correct and, at the same time, appropriate to 
the nature of the considered phenomena formulations of problems of 
mathematical physics. A realization of the outlined orientation seems 
to be possible in the context of the following considerations. 

The reasons for the difficulties related to the solution of ill-posed 
and essentially mathematically senseless problems are, in principle, well 
understood. In the Predholm integral equation of the first kind (1.1), 
there exists a mismatch between the function / (x) and a solution of 
the corresponding direct problem (the result of integration) , which is a 
result of errors in the determination of the data as well as of rounding 
of digits in arithmetical operations. 

As a consequence, considerable attention is paid to the phenomenon 
of (as it is sometimes called) smoothing of information about the func- 
tions in the process of their integration. At the same time, the data 
of the problem, i.e., the free term f (x) and also the kernel k{x,^), 
are usually determined experimentally, which inevitably incurs a con- 
siderable error in Eq. (1.1). In this regard, we should point out the 
dominance of the methodology of A. N. Tikhonov that is based on 
objective incorrectness of the formulation of most problems of mathe- 
matical simulation. 

There appears a rather obvious, as it seems, question: Why not take 
into account in practice the above-mentioned errors in the formulation 
of problems, instead of merely bearing them in mind when identifying 
the reasons for computational discrepancies? One can assume that an 
adequate simulation of the error may contribute to a correct formulation 
of the inverse problems. 

Here, the adequacy implies, in the first place, the functional struc- 
ture of the representation of the error. In this regard, let us turn to the 
procedure of integration (1.1). On the basis of general considerations, 
it is logical to represent the loss of information about the function ip (x) 
in the evaluation of / (x) in the form 

1 

{Sf){x)=i;{x)-Xj h{x,0i^{0d^, xe[0,l]. (1.2) 
-1 

Here, the function ip {x), x & [—1, 0), the kernel h {x, ^) and the param- 
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eter A should satisfy the requirement of the reahzation of the condition 

(Sf) (x) = (1.3) 

in the spaces C [0, 1] or L2 (0, 1) for ip{x), x G [0,1] from a rather 
representative class. 

Note that compared to the values of the sought and given functions 
the considered error is really small. Therefore, in the case of the con- 
struction of a stable algorithm of the evaluation of if) (x), its exclusion 
by condition (1.3) should not considerably influence the solution. 

The structure of the error (1.2) under the condition (1.3) embodies 
the difference between the function (a;), subject to integration, and 
its approximate expression that, in turn, appears as a result of the 
execution of an analogous procedure. One should note the absence of 
any a priori premises of self-sufficiency of (1.2) in achieving the goal, 
namely, a correct formulation of the problem of the determination of 
the function -0 [x] from the data (1.1). 

As a matter of fact, we put forward a hypothesis about the priority 
of a qualitative side of the phenomenon of smoothing of information 
in modelling the error of integration as well as, in general, about the 
expediency of the suggested "measures" for the realization of a correct 
formulation of the problem that is inverse to the procedure (1.1). 

On the basis of (1.2) and (1.3), instead of the ill-posed problem 
(1.1) for the determination of the function -0 (x), the following system 
of equations will be employed: 

/i(A^) (x)=/i/ {x) + {5f) (x); 

(5/)(x) = o, xe[o,i], (1.4) 

where is a parameter analogous to A. 

We have worked out two versions of the solution of the formulated 
problem. In the first version, the Fourier coefficients of the function 
■0 {x) G [0, 1] are represented in quadratures via the data of the prob- 
lem by means of the solution of the Predholm integral equation of the 
second kind that possesses rather favorable properties. This implies the 
absence of singularities or oscillations of the kernel that are not caused 
by k (a;, ^) (i.e., those that are enforced by the employed algorithm) as 
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well as the absence of a small factor explicitly multiplying the sought 
function ip (x). 

In the second version, the computational procedure reduces to a 
consecutive solution of two Predholm integral equation of the second 
kind, i.e. of the above-mentioned one and of another one that differs 
from the former one only by the form of the free term. We also demon- 
strate a possibility of the determination of the function ip (x) with the 
use of the solution of just the first of these two equations. 

The basis of the effectiveness of the performed transformations is 
formed by the following factors: 

1) An incompletely continuous perturbation of the operator A in 
combination, of course, with the condition (1.1) that led to the deriva- 
tion of the system (1.3) instead of (1.4). 

2) The extension of (1.4) to x e [—1, 0), which allows one to use a 
typical peculiarity of the solution of the Predholm integral equation of 
the second kind 



where k {x) is an undefined function stipulated by the form of the free 
term. 

3) The use of an equation, which is an analogue of (1.5), that pos- 
sesses the same solution on x E [0, 1] and a free term that goes to zero 
on the other part of the interval of definition [—1,0). 

4) The choice of the kernel h{x,C,) in such a way that be means 
of a linear change of the variables it could be transformed into the 
canonical Poisson's kernel with a parameter r, which allows one to do 
the following: 

- determine the function i(j (x) in the form of an expression that 
explicitly depends on r for the case when it is harmonic; 

- by proceeding to the limit r ^ 1, express the kernel and the free 
term of the above-mentioned Predholm integral equation of the second 
kind via the data (1.1); 

- as a result, extend the class of admissible belonging of the function 
ip (x) to the whole space L2 (0, 1) ( the first version of the solution of 
the problem). 




0, a;e[0,l]; 

k{x) , X e [—1, 0) , 



(1.5) 
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5) The consideration of the function ip (x) satisfying Eq. (1.1) as 
a generahzed function, which resulted in a considerable simplification 
of the procedure of meeting condition (1.3) (the second version of the 
solution of the problem). Simultaneously, the necessity of using the 
passage to the limit with respect to the parameter r lost its relevance. 

The objective of this work can also be explained by means of the 
following example. Consider a beam (bar) supported at the ends and 
subject to a transverse load: the problem consists in finding its defiec- 
tion (in the linear interpretation) . Correspondingly, using the notation 
of (1.1), we have: 

k{x,C,) is the deflection at the cross-section with the coordinate x 
caused by a unit force applied in the cross-section with the coordinate 

■0 (x) is the intensity of the distributed load; 

/ (x) is the defiection whose determination by integration according 
to (1.1) is successfully carried out by undergraduate university students 
taking a course in the strength of materials. 

However, the deflection of the beam is here to stay: it can be mea- 
sured; and the load does exist in reality. Therefore, fully justifled is 
the formulation of the inverse problem that consists in the determi- 
nation of t/-' (,t) from the given k{x,^) and f{x). Such a problem is 
considered to be of an incomparably higher degree of complexity, and 
the attempts of its solution is a preoccupation of not the students, but 
rather of scientists and, in particular, of their lecturers. In these at- 
tempts, the Fredholm integral equation of the flrst kind (1.1) is used, 
whose solution, in reality, cannot be realized. Moreover, even obtaining 
a palliative implied by this sohition requires application of great efforts. 

At the same time, it is reasonable to suggest that the difficulties arise 
due to the fact that the problem is ill-posed, as was explained above. 
Really, there exists, on the other hand, a very convenient, from the 
point of view of a numerical realization, object, namely, the Fredholm 
integral equation of the second kind that can be represented in the form 

1 

ij{x)=f^jK{x,Oi^{Od^ + F{x), xe[0,l], (1.6) 



where the kernel K [x, ^) and the free term F (x) are given; the function 
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ip (x) (the intensity of the load) is subject to determination; is a 
parameter that has to be chosen from the solvabiUty condition. 

However, the question arises: What are the reasons to beheve that 
the function (x), provided it is the same as that entering (1.1), should 
satisfy this equation, and what is understood under K (x, ^) and F (x)l 
On the other hand, if the function ip (x) is assumed to be known and the 
kernel K {x, ^) is given even in an arbitrary form from the space L2, one 
can always find a free term F (x) allowing one to satisfy Eq. (1.6). As a 
consequence, a Predholm integral equation of the second kind satisfied 
by the sought function ip (x) objectively exists.^ Moreover, the number 
of such equations is not limited. 

The construction of Eq. (1.6), simultaneously with Eq. (1.1) satis- 
fied by the function ip (x), is, in fact, is the objective of this work. In 
other words, it is devoted to the construction of the free term F (x) de- 
pending on the data (1.1) in such a way that the function ip (x) and the 
solution of (1.6) coincide. Thus, equation (1.6) is a well-posed problem 
for the determination of the function ip (x) satisfying (1.1). 

However, a broad class of linear boundary- value and initial-boundary- 
value problems of mathematical physics can be rather elementarily re- 
duced to the Predholm integral equations of the first kind. To this end, 
considering, for example, a problem described by the Laplace equation, 
one has to set d'^u = ip (or d'^u + j3u = ip, where j3 is a constant). 

By means of integration with respect to x the function u {x, y) and 
its derivatives are expressed via ip. Integration of the differential equa- 
tion with respect to y allows one to obtain a second representation of 
u via Ip. The one-dimensional functions of integration in these rep- 
resentations are expressed via ip and satisfy the boundary conditions. 
The elimination of u from the representations of the solution leads to 
a two-dimensional Predholm integral equations of the first kind with 
respect to the function ip {x, y). 

The outlined scheme is practically indifferent to the type and or- 
der of differential operators, the presence of variable coefficients, the 
configuration of the boundary of the domain of the function and some 
other factors that usually complicate the realization of numerical algo- 



^One can assume that analogous arguments formed the basis of J. Hadamard's 
statement. 
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rithms. The above-mentioned method of the solution of the problem 
(1.1) is directly extended to the determination of the function ip {x,y) 
(the variable y plays the role of a parameter). In this regard, there 
appears an interesting possibility to check the solvability of the prob- 
lems of mathematical simulation that can be represented by partial 
differential equations. 

Perhaps, the motivation of the proposed investigation could be of 
certain interest. The reason was the confusion caused by the absence 
in the specialized literature of a clear statement of the universality 
of the outlined method of the reduction of problems of mathematical 
physics to the Fredholm integral equations of the first kind. However, 
a placement of the whole lot of initial data in their kernels is neverthe- 
less rather attractive. Indeed, a conventional classification of problems 
according to the complexity of their numerical realization is, in fact, 
depleted, and the construction of an effective method of the determi- 
nation of the functions satisfying the stated type of equations comes to 
the foreground. 

In Chapter 2, we analyze J. Hadamard's arguments concerning the 
issue of correct formulation of problems for partial differential equa- 
tions. Both related and alternative positions on this issue of known 
specialists are illuminated. We also discuss Banach's theorem on in- 
verse operator that is closely related, in a contextual sense, with the 
methodology of correctness. Arguments are given that mathematical 
formulations of the direct and the inverse problems should not be iden- 
tical. 

Chapter 3 contains an analytical review of the methodological ap- 
proaches and methods of the solution of ill-posed problems (mostly, of 
Fredholm integral equations of the first kind) related to the concepts 
of A. N. Tikhonov and V. M. Fridman. Pointed out are some expert 
opinions about rational use of digital information and, on the whole, 
about priorities of the development of computational mathematics. 

The material of Chapter 4, in a sense, refracts principle difficul- 
ties, accompanying the solution of ill-posed problems by the prism of 
fundamental concepts of J. Hadamard and S. Banach. We present argu- 
ments for inconsistency of the methodology of the solution of ill-posed 
problems. General premises for a correct formulation of the problem 
of determination of the function satisfying the Fredholm integral equa- 
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tions of the first kind arc given. 

Chapter 5 is devoted to the construction of the method of the re- 
duction of the problems, usually associated with the Predholm integral 
equation of the first kind, to the solution of Predholm integral equa- 
tions of the second kind. This section is. in a constructive sense, basic. 
Exactly here we consecutively construct an algorithm that practically 
realizes the main factors ensuring the efficiency of the transformations 
pointed out above. The first version of the solution of the problem is 
presented. 

In Chapter 6, we emphasize the main points of the carried OTit trans- 
formations and also study a possibility of their variation. An interpreta- 
tion of the algorithm of the previous section is given from a generalized 
point of view. Furthermore, the second version of the considered prob- 
lem is given. Its correct formulation in terms of the Predholm integral 
equation of the second kind is presented. 

The material of Chapter 7 iUustrates the universality of the tech- 
nique of the reduction of linear boundary-value and initial-boundary- 
value problems to Fredholm integral equations of the first kind. An 
extension of the suggested algorithm of the solution of the equations 
of this type to a two-dimensional case is demonstrated. The issue of 
solvability of the problems of mathematical simulation is touched on. 

Chapter 8 develops the outlined orientation of the reduction of the 
problems to the Fredholm integral equation of the first kind involv- 
ing into the sphere of transformations sufficiently nontrivial applica- 
tions (including factors of nonlinearity, singular perturbations and some 
other). The presentation of the material has the form of sketches. 

In Chapter Conclusions, we summarize the main points of the work 
from the same position of priority importance of correct formulation of 
problems of mathematical simulation for the efficiency of their numer- 
ical realization. 

Mathematical techniques employed in the presentation of the ma- 
terial is comparatively simple: basics of the classical theory of integral 
equations; elements of functional analysis; general principles of formu- 
lations of problems of mathematical physics and of methods of their 
solution. When performing transformations, we often refer to the book 
by F. G. Tricomi, Integral Equations (Dover, New York, 1957). 

The literature to each chapter is given in reference order. (Note 
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that page numbers refer to the Russian edition of a corresponding ht- 
erature source.) Chapters and sections (chapters and paragraphs of the 
hterature sources) are referred to, respectively, as Chapter 1, section 
1.1, sections 1.1, 1.2. 

The numbering of formulas in the text is dual: the first numeral 
refers to the chapter number, whereas the second one refires to the 
formula number inside the chapter. 

The author is most grateful to I. Zhuravlev for pointing out a con- 
tradiction in the transformations of Chapter 5. As a result, this section 
was substantially revised. Some revisions were made also in other chap- 
ters. However, the methodology of the work and the basis of the method 
of the solution did not change. 

I also thank I. Stepanov for making a web-site in the Internet, M. 
Katchamanova for her assistance in the preparation of the manuscript, 
and S. V. Kuplevakhsky for translating the Russian version of the 
manuscript (see: www.pelbook.narod.ru) into English. 



Chapter 2 



The issue of the correct 
formulation of problems of 
mathematical physics 

2.1 Hadamard's definition of correctness 

J. Hadamard has defined two conditions that should be satisfied by a 
correctly formulated boundary-value (initial-boundary-value) problem 
for partial differential equations: existence and uniqueness of the solu- 
tion ( [1], p. 12).^ At the same time, the third condition of Hadamard's 
definition of correctness that concerns continuous dependence on the 
data of the problem is well-known. Indeed, he paid serious attention 
to the investigation of this issue with regard to Cauchy-Kovalevskaya's 
theorem concerned with the solution of the differential equation 

= f (t, xi, X2,..., Xn, dtu, dx^u, dx^u, d^^uj (2.1) 

(a system of analogous equations), where / is an analytical function of 
its arguments in the vicinity of the origin of coordinates, with initial 
conditions 

U (0, Xi, X2,..., Xn) = (fiO {Xi, X2,..., X„) ; 

^For the first time, the concept of correct formulation was put forward by 
Hadamard in his article of 1902. 
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dlu (0, Xi, X2, . . . , Xn) = ^Ps {X1,X2, ■ ■ ■ , X^) , 

s^l,...,k-l. (2.2) 

As is pointed out by Hadamard, the consideration of the problem 
(2.1), (2.2), named after Cauchy, raises three questions ([1], p. 17): 

1) Does it admit a solution? 

2) Is the solution unique? (In general, is the problem well-posed?) 

3) Finally, how the solution can be derived? 

Cauchy-Kovalevskaya's theorem (in its authors' interpretation) states 
that, except for some special cases, the above-mentioned problem ad- 
mits a unique solution that is analytical at the origin of coordinates. 
Moreover, the functions (po,..., (fk-i in (2.2) can be not only analytical 
but regular, i.e., continuous together with their derivatives up to a cer- 
tain order. A possibility of a uniform approximation of (fo,..., (fk-i by 
Taylor series expansions in powers of ,Xi, , retaining all operations 

on analytical functions, including differentiation up to a corresponding 
order, is implied. 

However, such an approach was strongly criticized by Hadamard. 
In his opinion, the question is not how such an approximation affects 
the initial data, but rather what is an effect on the solution? He em- 
phasized the non-equivalence of the notion of small perturbation for 
given Cauchy 's problem and of the solution to this problem ([1], p. 
39). In this regard, J. Hadamard presented his prominent example of 
a solution of the differential equation 

d^u + dlu = 0, (2.3) 

subject to the conditions 

u{x,0) — 0; dtu {0, x) — an sin (nx) , (2-4) 

where is a rapidly decreasing function of n. 

The expression on the right-hand side of (2.4) can be arbitrarily 
small. Nevertheless, the problem admits the solution 

u {x, t) — — sin (nx) sinh (nt) . (2.5) 

For an = 1/n or or e~^, this solution is rather large for any 

nonzero t, because of the prevailing growth of e"* and, correspondingly. 
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of sinh(nt). Thus, the function (2.5) does not depend continuously on 
the initial data and, as a result, the problem (2.3), (2.4) is ill-posed. 

Concerning the regularity of the right-hand side of (2.2), J. Hadamard 
remarked: " ...actually, one of the most curious facts of the theory is that 
equations, seemingly very close to each other, behave in a completely 
different way" ([1], p. 29). 

A large number of investigations devoted to the issue of the correct 
formulation of Cauchy's problems. The authors of these investigations 
concerned themselves with specification of corresponding classes of dif- 
ferential equations and with minimization of requirements imposed on 
the initial data (sec [2]). However, we are mostly interested in the 
actual character of the dependence of the solution on the data of the 
problem and, in this regard, the classic J. Hadamard's statement that 
"an analytical problem is always well-posed in the above-mentioned 
sense, when there exists a mechanical or physical interpretation of the 
question" ([1], p. 38). 

As was pointed out by V. Y. Arsenin and A. N. Tikhonov [3], the 
latter questioned the legitimacy of studies of ill-posed problems, speci- 
fied by the authors as the following: the solution of integral equations of 
the first kind; differentiation of approximately known functions; numer- 
ical summation of Fourier series whose coefficients are approximately 
known in the metric I2', analytical continuation of functions; the so- 
lution of inverse problems of gravimetry and of ill-defined systems of 
linear algebraic equations; minimization of functionals for divergent se- 
quences of coordinate elements; some problems of linear programming 
and of optimal control; the design of optimal systems and, in partic- 
ular, the synthesis of aerials. It is emphasized that this list is by no 
means complete, because ill-posed problems appear in investigations of 
a broad spectrum of problems of physics and engineering. 

In his talk at the meeting of the Moscow Mathematical Society de- 
voted to Hadamard's memory, G. E. Shilov said the following [4]: "Our 
time has brought aboTit corrections in Hadamard's instructions, be- 
cause it turned out that ctill-posed, according to Hadamard, problems 
could have meaning (as, e.g., the problem of restoration of a potential 
from scattering data). However, the studies of well-posed problems, 
proclaimed by Hadamard, was a cementing means for the formation of 
the whole theory" (functional analysis is implied). This quotation is 
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borrowed from a biographical sketch by E. M. Pohshtuk and T. O. Sha- 
poshnikova [5] , where it is also pointed out that in the course of time 
J. Hadamard's opinion about the importance for practice of exclusively 
well-posed problems was understood in a less absolute sense. 

At the same time, rather sharp statements were made: 

"And what is more, Hadamard put forward a statement that ill- 
posed problems had no sense at all. Since (as can be seen from a modern 
point of view) most applied problems, represented by equations of the 
first kind, are ill-posed, this statement of the outstanding scientist, 
apparently, strongly slowed down in 1920-1950's the development of 
the theory, methods and practice of the solution of problems of this 
class" ([6], p. 12). 

"Until quite recently, it was thought that ill-posed problems had no 
physical sense and that it was unreasonable to solve them. However, 
there are many important applied problems of physics, engineering, 
geology, astronomy, mechanics, etc., whose mathematical description is 
adequate although they are ill-posed, which poses an actual problem of 
the development of efficient methods of their solution" ([6], p. 225).^ 

"Prom the results of this work [of A. N. Tikhonov] followed a limita- 
tion of the well-known notion of J. Hadamard [1] of a well-posed prob- 
lem of mathematical physics, which was of indisputable methodological 
interest, and inconsistency of Hadamard's thesis, wide-spread among 
investigators, that any ill-posed problem of mathematical physics was 
unphysical." ([7], p. 3). 

"For a long time, activities related to the analysis and solution of 
problems called ill-poscd used to be relegated (by famous mathemati- 
cians too) to the domain of metaphysics" ([8], p. 126). "A prevailing 
number of mathematicians (including Hadamard) expressed. their atti- 
tude towards this problem in the following way: If a certain problem 
does not meet the requirements of correctness, it is of no practical 
interest and, hence, does not need to be solved" ([8], p. 127) (I. G. Pre- 
obrazhenskii, the author of the section "Ill-posed problems of mathe- 
matical physics"). 

Note that the latter paper most distinctively reveals the style that 
causes a principal objection. Thus, A. Poincare is accused of incon- 



^In the context of what follows, we draw attention to the "adequate description" . 
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sistency of methodological views on the nature of causal relationship 
([9]) ("The Last Thoughts"). Indeed, the text does not contain any 
evidence that he makes a fetish of the problem of restoration of the 
cause from the effect. On this basis, a conclusion is made about the 
great scientist's misunderstanding of the essence of instability of com- 
putation procedures inherent to ill-posed problems and, in particular, 
to integral equations of the first kind.^ 

The adequacy of employed models to considered concrete processes 
is not even touched on by the authors of ([8]). Thus, a quite legiti- 
mate question arises: How does one know that Poincare, if necessary, 
could not find a way of a mathematically correct formulation of the 
same physical problems? Anyway, is there any contradiction in general 
arguments for the existence of such a possibility, including the aspects 
of its constructive realization? 

By the way, exactly Poincare repeatedly mention Hadamard while 
establishing a relationship between the correct formulation of problems 
and a practical realization of employed models. Wc draw attention to 
an expressive thesis: "If a physical problem reduces to an analytical 
one, such as (2.3), (2.4), it will seem to us that it is governed by a pure 
occasion (according to Poincare, it means that determinism is violated) 
and it does not obey any law" ([1], p. 43). 

In light of the above, the arguments of I. Prigoginc and I. Stcngcrs 
[10] are of interest: "...one can speak of a 'physical law' of some phe- 
nomenon only in the case when this phenomenon is 'coarse' with respect 
to a limiting transition from a description with a finite accuracy to that 
with an infinite accuracy and thus inaccessible to any observer, whoever 
he may be" (p. 9). "Scientist in a hundred different ways expressed 
their astonishment that a correct formulation of the question allows 

^In particular, the exact statement reads: "However, one must remember that 
vagueness of philosophical positions of some scientists in the West, even rather 
renown, results in the fact that, based on correct starting points, they draw rather 
inaccurate conclusions, repeating old mistakes of, for example, A. Poincare, who 
writes: 'If two organisms are identical, or simply similar, this similarity could not 
occur by chance, and we can assert that they lived under the same conditions...' In 
other words, the fact of possible incorrectness of the inverse problem is completely 
ignored." However, one would hardly mention Poincare's mistakes... if modern 
" spontaneous supporters of the principle of determinism did not repeat them" ( [8] , 
p. 134). 
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them to solve any puzzle suggested by nature" (p. 44). 

Thus, underlined are: first, methodological importance of correct 
formulation of problems of mathematical physics; second, a leading 
role of the employed procedures and, finally, substantial influence of 
the quality of their realization on the degree of complexity of obtaining 
the final result. In other words, one and the same problem can be 
better or worse formulated. 

The above-mentioned statement or Hadamard's postulate, as it 
called by S. K. Godunov ([11], p. 113), as a matter of fact, imphes a 
possibility of a "good" (correct) formulation of any meaningful problem 
and, consequently, can be interpreted as having a global orientation.^ 

In this regard, one can establish an obvious relationship to D. 
Hubert's comments on his 20th problems that suggested a possibility 
of correct formulation of arbitrary boundary- value problems of mathe- 
matical physics by means of special requirements on boundary values 
of corresponding functions (a type of continuity or piecewise differen- 
tiability up to a certain order) and, by necessity, by giving an extended 
interpretation to the notion of the solution ([12], pp. 54-55). 

For the first time, the three conditions of the correctness of problems 
of mathematical physics were clearly pointed out by D. Hilbert and R. 
Courant ( [13], pp. 199-200): existence, uniqueness and continuous 
dependence of the solution on the data of the problem. Concerning 
the last, they say: "...it has crucial importance and is by no means 
trivial... A mathematical problem can be considered adequate to the 
description of real phenomena only in the case when a change of given 
data in sufficiently narrow limits is matched by an alike small, i.e. 
restricted by predetermined limits, change of the solution" . 

V. A. Steklov's position is quite analogous ([14], p. 62): "...if differ- 
ential equations with the above-mentioned initial and boundary condi- 

^The above mentioned reference contains the following definition: 
A problem is called well-posed if is solvable for arbitrary initial (or boundary) data 

belonging to a certain class, has a unique solution, and this solution continuously 

depends on the initial data. 

A problem is called ill-posed either if it is not solvable for arbitrary initial data 

or if it is impossible to choose such norms for the solution and for the initial data 

that continuos dependence of the solution on the data of the problem with respect 

to these norms be ensured. 
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tions are not constructed on erroneous grounds, are not in direct con- 
tradiction to the reality, they must yield for each problem a unique and 
completely definite solution...". Along the same hnes, I. G. Petrovskii 
writes ([15], p. 87): "The above-mentioned arguments for the correct 
formulation of Cauchy's problem show that other boundary- value prob- 
lems for partial differential equations are of interest for natural science 
only in the case when there is, in a sense, continuous dependence of the 
solution on boundary conditions" . 

S. L. Sobolev is less categorical ([16], p. 38): "The solution to an 
ill-posed problem in most cases has no practical value" . Of considerable 
interest is the opinion of V. S, Vladimirov ([17], p. 69): "The issue of 
finding correct formulations of problems of mathematical physics and 
methods of their solution (exact or approximate) is the main content 
of the subject of equations of mathematical physics" . 

V. V. Novozhilov, in fact, drew attention to the potential of vari- 
ation of the formulation of the considered problem with the aim of 
the simplification of the procedure of its numerical realization ([18], p. 
352): "The absence in the term "a mathematical model" of the indi- 
cation of its inevitable approximate character leaves way for a formal 
mathematical approach to models, disregarding those concrete prob- 
lems for whose solution they were intended, which is, unfortunately, 
wide-spread at present" . 

2.2 J. Hadamard's postulate and incor- 
rectness of "real" problems 

Thus, J. Hadamard and a number of other outstanding scientists thought 
that any physically interpretable problem could be well-posed. How- 
ever, a quite opposite point of view dominates in modern publications. 
Indeed, a visibly larger part of practically important problems con- 
sidered therein are incorrect. However, is the actual methodology of 
mathematical formulation of these problems and, correspondingly, the 
results of its refraction with respect to realities adequate? 

Here we will not elaborate on something like general principles of 
the construction of differential equations, and, generally speaking, it 
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is reasonable at the beginning to restrict the question to the follow- 
ing: What arguments allow one to conclude that an ill-posed problem 
adequately describes an observable phenomenon or a potentially real 
process? In this regard, let us turn to the procedure of the solution of 
the Predholm integral equation of the first kind 

1 

jk{x,i)i^{i)di^f{x), xe[0,l], (2.6) 



which is a classical incorrect problem: the square summable kernel 
k{x,^) and the free term / [x] are given; the function -0 (x) is to be 
restored. 

Let us assume that the kernel is symmetric and closed, i.e. k{x,^) = 
k (^,x) and its eigenfunctions (x), being nontrivial solutions of the 
integral equation 

1 

i^{x)^\j k{x,0^{OdC-f{x), [0,1] 



with characteristic numbers A = A„, n = 1,2,..., form a complete in 
L2 (0, 1) orthogonal system of elements. In this case, the solution to 
Eq. (2.6) exists and is unique under the condition (see, e.g., [19], pp. 
185-187) 

00 J, 

<oo, f{x)'ijjn {x) dx. (2.7) 

n=l ^ 

If all the above-mentioned conditions are fulfilled, there is still the 
third condition of correctness that, as it is known, is certainly not sat- 
isfied by Eq. (2.6). Numerous literature references clearly illustrate an 
inadequately strong infiuence on the solution of small perturbations of 
the data of the problem, in the first place of / (x). As a rule, this func- 
tion is determined experimentally and mismatch the kernel k{x,^), 
in particular, with respect to smoothness. Thus, Eq. (2.6), strictly 
speaking, looses sense. At the same time, a possibility of an equiva- 
lent description of the problems of mathematical physics by means of 
integral equations of the first kind is indisputably admitted at present, 
which is confirmed by their colossal list ([6], section 4.2). 
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Let us specify Eq. (2.6): 

(2.8) 

where m is an integer. For this choice, A„ = {mn) ; ipn (x) = \/2sin (mix); 
n = 1,2,... ([20], p. 149). 

Since the kernel k is symmetric and continuos, and all A„ > 0, the 
use of Mercer's theorem [19], according to which 

n=l '^n 

and a representation of ip (x) as a series expansion in terms of ipn (x) 
with undetermined coefficients allows one to find the solution to Eq. 
(2.6): 

ip (x) — sin (rmrx) . (2-9) 

However, the procedure of calculations turned out to be so simple 
owing to a special choice of the data of the problem. If this is not the 
case or in the case of the solution of Eq. (2.6) with the kernel and the 
free term (2.8) by means of one of numerical methods, the complexity 
of the realization of a an approximation of sufficiently high order is 
practically identical to the most general situation, characterized by an 
error in the determination of / (x) and k (x,^).^ As a matter of fact, 
even if the data are objectively compatible, the incorrectness of Eq. 
(2.6) appears as a result of rounding off the digits in the process of 
calculations. 

The factor of the incorrectness of Eq. (2.6) follows from a compari- 
son of the free term (2.8) with the solution (2.9). Indeed, by increasing 
m, the function / (x) may turn out to be arbitrarily small, whereas 
the bounds of the values of i/j (x) arc unchanged. Correspondingly, any 
error in the calculations with / {x) is projected onto the function ip (x) 
with the factor m^. The mechanism of this phenomenon of the smooth- 
ing of information about the function in the process of integration will 
be repeatedly discussed in what follows. 



^Here, complexity implies an ill definition of the system of linear algebraic equa- 
tions obtained as a result of some sort of discretization. 
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However, let us return to the question of the relation of an incorrect 
formulation to the reality. In this regard, we draw attention to the 
following. By considering (2.6) as the Predholm integral equation of the 
first kind , we mean the solution of the inverse problem (I). However, 
equation (2.6) can be used for the solution of the corresponding direct 
problem (D): the determination of the function f (x) from the data 
k{x,^) and ip (x). This procedure is correct and thus is radically simpler 
than the problem I. It is sufficient to note the absence of any principal 
difference between the evaluation of the integral (2.6) in an analytical 
form and its essentially numerical realization. 

Here we want to draw attention to an issue that seems to be of 
substantial importance. The problem D, as a rule, is transparent: in 
its categories, we adequately model realistic current processes and phe- 
nomena by, which should be emphasized, explicit means of linear su- 
perposition. Correspondingly, if, for instance, k {x, ^) is a characteristic 
of the system and ip (x) is intensity of external influence, a resulting 
effect in this or that subject sphere is to be elementarily summed up. 

The situation is diametrically different for the problem I. One could 
hardly point out any realistic process (phenomenon) for which it could 
be formulated in mathematical terms directly on the basis of the subject 
sphere. In other words, without any relation to the problem D, which 
commonly implies a transformation of the latter into the problem I just 
by means of renaming of known and unknown components. 

It seems that the methodology, which states the adequacy of the 
problem I, obtained by the above-mentioned renaming of the compo- 
nents, to the realities on the basis of a high-quality information about 
a concrete problem D, is profoundly deficient. Correspondingly, the 
opinion of experts who a priori reject J. Hadamard's argument for the 
existence of correct formulations of problems of mathematical physics 
should be considered unjustified. 

Let us turn to the problem D that describes some realistic process 
(phenomenon) (2.6). For this process, the determination of ip (x) from 
the data k{x,^) and f (x), i.e. the formulation of the corresponding 
inverse problem that will be denoted as I , is, of course, reasonable. 
Suppose that in this case Hadamard's argument holds, and, hence, the 
problem I is correct. However, the problem I, the solution of the Fred- 
holm integral equation of the first kind (2.6), is ill-posed by definition. 
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The conclusion is obvious: Mathematical formulations (represen- 
tations, expressions) of the problems I and I are non-identical. As a 
result, the formulation of the problem l' cannot be restricted to read- 
dressing the status of the unknown variable between the functions / 
and ip in the problem D. Note in this regard that a development of the 
methodology of the correct formulation of the problem that is inverse 
to D, i.e. I , is the main objective of the present investigation. 

The above arguments seem to be rather convincing, however, at 
this stage of our consideration, we can neither prove the correctness 
of Hadamard's postulate (argument) in the general case nor illustrate 
its constructive character with respect to separate classes of problems. 
One should also bear in mind that, using special methods, the solution 
of the ill-posed problem I (or what is understood under the solution), 
as a rule, can be obtained with accuracy that is considered to be prac- 
tically acceptable. In this regard, the question arises: Should one aim 
at the correct formulation I , if the algorithm of the calculation of the 
function ijj (x) in the formulation of the problem I in some way realizes 
its regularization? This implies a well-known deformation of the for- 
mulation I with the use of a small parameter that yields the property 
of correct solvability. 

Thus, can the algorithm to a full extent, including the efficiency of 
numerical realization, level off the principal difficulties inherent to the 
incorrectness of the problem I in the form (2.6)? It is clear that the 
answer is definitely negative: Otherwise, the deep-rooted differentiation 
between ill-posed problems and well-posed ones would make no sense. 

Furthermore, the indicated difference is of exceptional importance, 
because correctness of the formulation is a criterion of a qualitative 
level, whereas the efficiency of a method of the solution of the Fred- 
holm integral equation of the first kind can be estimated only in terms 
of quantitative factors of a palliative property. The latter is caused by 
a direct relationship between a degree of regularization and the defor- 
mation (distortion) of the problem I. 

What is, however, the actual difference in the interpretation of the 
formulations I and I ? The answer to this question is contained in sec- 
tions 4.5, 5 and 6. At this stage, we only note that a transformation of 
the formulation I into the formulation I will be realized by means of 
an incompletely continuos perturbation of the integral operator of the 
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problem (2.6) that simulates the phenomenon of smoothing of informa- 
tion. 

2.3 Banach's theorem on the inverse op- 
erator 

Let us quote ([5], p. 175): "First, Hadamard defined the correctness of 
the problem by the conditions of solvability and uniqueness and strongly 
insisted on continuous dependence of the solution on the initial data 
only in the consideration of Cauchy's problem. In the book 'The theory 
of partial differential equations', published in Peking a year after his 
death, he wrote: 'This third condition that we introduced in "Lectures 
on Cauchy's problem..." but did not consider as part of well-posed 
problems, was added, quite justified, by Hilbert and Courant [13]. Here, 
we accept their point of view.'" 

E. M. Polishuk and T. O. Shaposhnikova made the following com- 
ment on this text [5], pp. 175-176]: "Prom a mathematical point of 
view, the question of the necessity of the requirement of the continuity 
of the solution with respect to the data seems to be rather delicate. As 
a matter of fact, according to Banach's well-know theorem on closed 
graph, unique solvability of a linear problem leads to boundedness of 
the inverse operator and, thus, continuous dependence of the solution 
on the right-hand sides." It is pointed out that variations of the coef- 
ficients of differential equations and of the boundary of the considered 
domain can also influence the solution of the problem; hence, the use 
of the three conditions of the correctness is preferable. 

At the same time, Banach's theorem on the inverse operator ([21], p. 
34), being a consequence of the above-mentioned one, is more closely 
related to the considered issue. Its formulation, given by A. I. Kol- 
mogorov and S. V. Fomin, is the following ([22], pp. 259-260): Let A 
be a linear bounded operator that maps a Banach space Bi in a one- 
to-one fashion onto a Banach space B2. Then the inverse operator A~^ 
is unique. 

In addition, L. A. Lyusternik and V. I. Sobolev ([23], pp. 159- 
161) emphasized that a one-to-one mapping of the whole Banach space 
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Bi onto the whole Banach space B2 is imphed. Besides, a situation is 
discussed when " ...an operator, being the inverse of a bounded operator, 
although linear, turn out to be defined not on the whole space B2 but 
only on a certain linear manifold and unbounded on this manifold" . 

The formulation of the same theorem in ([24], p. 60) reads: If a 
linear bounded operator A that maps a Banach space Bi onto a Banach 
space B2 has an inverse A^^, then A~^ is bounded. It is pointed out 
that this statement becomes invalid if one gives up the requirement of 
completeness of one of the spaces. There is also a clarification: The 
existence and uniqueness of the solution of the equation Aip = f with 
an arbitrary right-hand side from B2 leads to continuous dependence 
of the solution ip = A'^^f on /. 

S. Banach himself made the following statement: If a hnear opera- 
tion realizes a one-to-one transformation of Bi onto B2 , the transforma- 
tion is mutually continuos. At the same time, in the formulation of the 
theorem on the closed graph, he pointed out that the transformation 
Bi is realized onto the whole space B2. 

L. V. Kantorovich and G. P. Akilov made a refinement concerning 
a mapping under the specified conditions onto a closed subspace of the 
Banach space B2 ([25], p. 454). The essence is that a closed subspace 
of a Banach space is itself a Banach space. 

S. G. Mikhlin gave a proof of the theorem ([26], p. 507): For the lin- 
ear problem Aip = / to be well-posed in a pair of Banach spaces Bi, B2, 
it is necessary and sufficient that the operator A~^ exist, be bounded 
and map the whole space B2 onto B^. At the same time, a clear dis- 
tinction is made between the category of the existence and uniqueness 
of the solution of the boundary- value problem and its correctness as a 
whole, which implies, as a result, continuous dependence on the data 
(the third condition according to Hadamard) . The following definition 
is given: "A boundary- value problem is called well-posed in a pair of 
Banach spaces Bi, B2 if its solution is unique in Bi and exists for any 
data from B2, and if an arbitrarily small change of the solution in the 
norm B2 corresponds to a sufficiently small change of the initial data 
in the norm Si" (p. 204). 

The author pointed out that the problem might turn out to be 
well-posed in one pair of spaces and ill-posed in another one. Besides, 
the fact that the Fredholm integral equation of the first kind (2.6) is 
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ill-posed follows from the contradiction: If the problem is well-posed, 
there exists a bounded operator A^^ and, hence, the identical operator 
I — A~^A is completely continues in the corresponding infinite dimen- 
sional space, which contradicts the fundamentals of the general theory 
[24]. Mikhlin also quite encouragingly pointed out the approach of an 
approximate solution of ill-posed problems headed by A. N. Tikhonov. 

In an analogous, as to its content, course ([27], pp. 169-170), 
Mikhlin reiterated the above-mentioned formulations. However, Tikhonov 
is not mentioned at all, whereas the discussion of Eq. (2.6) found a 
rather interesting continuation (p. 171). It is shown that the prob- 
lem of its solution becomes well-posed if the pair of spaces Bi, B2 is 
replaced with such one that the operator A is no longer completely 
continues. The general considerations are illustrated by the following 
example. Let k {x, ^) and / (x) satisfy the conditions of section 2.2, in- 
cluding (2.7). It turns out that if one retains L2 (0, 1) as Bi and for B2 
also takes a Hilbert space of functions normalized according to (2.7), 
the solution of Eq. (2.6) becomes a well-posed problem: the operator 
A is incompletely continues and ^4"^ is bounded. 

In this regard, one can point out that the operator A is restrictively 
invcrtiblc not only when it acts from Bi onto the whole Hilbert space 
B2. It is sufficient that the operator A be bounded from below and that 
its range R (A) be dense everywhere in B2. At the same time, R (A) is 
not necessarily closed ([28], p. 34). 

A decade later, Mikhhn, in fact, gave up the investigations related to 
the issue of correctness [29]: "The author adheres to the classical point 
of view, according to which the problem being solved by mathematical 
methods should be considered as well-posed. Of course, there are other 
opinions (p. 7)... Thus, we neglect the so-called incorrigible errors 
related to the formulation of the above-mentioned problem as a problem 
of natural science or of social studies (measurement errors, insufficient 
accuracy of basic hypotheses, etc.)" (p. 17). 

M. M. Lavrentiev and L. Y. Saveliev characterized investigations of 
the issue of the solvability of Eq. (2.6) on the basis of considerations 
of the type of [27] as trivial, because it is difficult to imagine that for 
experimentally determined / (x) the corresponding error may prove to 
be small in the norm of the space B2 ([30], p. 217). At the same time, 
it is pointed out that, generally speaking, for any operator equation. 
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one can choose pairs of spaces such that the problem of its solution will 
be well-posed. 

G. M. Vainikko and A. Y. Veretennikov draw attention to the com- 
plexity of the description of such spaces. Thus, even the Volterra inte- 
gral equation of the first kind 

X 

jk{x,i)i;{i)di^f{x), xe[o,i], 



which admits the regularization 

X 

^{x) + Jd,k{x,Oi^{Od^ = f'i^), xe[o,i] 



and is elementarily solvable by quadratures, for reasons of the norm for 
ip {x), as a rule, has to be considered as an ill-posed problem ([31], p. 
6). 

As regards the pair of spaces that realize the conditions of the cor- 
rect formulation, an original remark of K. I. Babenko is of interest ([32], 
p. 304): "Hadamard's well-known example (2.3), (2.4) that yields the 
solution of Cauchy's problem of the type (2.5) by no means tells of the 
absence of continuous dependence on the initial data, as it is usually 
interpreted. It rather tells of the fact that small changes of the initial 
data may result in leaving the totality of the initial data for which the 
solution of Cauchy's problem exists." 

By the way, R. Richtmyer demonstrated the correctness of the pro- 
cedure of a numerical realization of a rather complicated problem of 
the above-mentioned type with the representation of sought functions 
by two-dimensional power series and with the use of special methods 
of suppression of errors of arithmetical operations ([33], section 17.B). 

In the course of V. A. Trenogin ([34], p. 225), the following two 
theorems are given: 

Let El and E2 be infinite dimensional normed spaces, with E2 being 
complete. If A is a completely continuos linear operator from Ei into 
£'21 different from a finite dimensional one, its range R{A) is not a 
closed manifold in E2. 
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Let A he a completely continuos operator from an infinite dimen- 
sional normed space Ei into a normed space E2, with the inverse oper- 
ator A~^ existing on R (A). Then A~^ is bounded on R (A). 

2.4 The premises of the reahzation of the 
conditions of correctness 

Let us assume that f = f* {x) is an exact result of integration of the 
function [x] e L2 (0, 1) and of the symmetric closed kernel k {x, ^) by 
means of the formula (2.6): 

{A^){x) = U{x), a;e[0,l]. (2.10) 

However, in general, the right-hand side of this equation is actually 
the following: 

/ (x) = /, (x) - (5/0 (x) , (2.11) 

where 5f' is an admissible error. ^ Quite naturally, / does not belong 
to the range of the operator A defined by the condition (2.7). In what 
follows, the space of functions for which this condition holds will be 
denoted 1'2. In contrast to the usual Z2, the property of belonging to 
depends both on the function / (x) and on the operator A. 

Thus, ^2 is the space of functions obtained as a result of integration 
according to (2.6) of the given kernel k (x, ^) and of the whole set ip {x) 
from L2. As a matter of fact, ['2 and R {A) coincide in the considered 
case. At the same time, the notion of the space ^3 characterizes to a 
greater extent the form of the normalizing functional (2.7). Between I2 
and such an abstraction as as the range of the operator A is Picard's 
theorem that provides the condition of the solvability of Eq. (2.6) [20]. 
Moreover, it may prove to be useful to compare with the space I2 
whose close relationship with L2 is estabhshed by the Riesz-Fischer 
theorem [19]. 

On the contrary, the free term of Eq. (2.10) {x) G Sia(\, at 
the same time the only verification of the condition (2.7) may prove to 
be infeasible because of the accumulation of errors of the calculations. 



^Here, the prime is used to match the notation of section 4.5 and thereafter. 
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Specific "diffusion" of the space I2 is caused by the structure of its nor- 
malizing functional. In this sense, the L2 (0, 1) is much more tangible 
for the function / (x) . Nevertheless, the use of it incurs rather negative 
consequences. 

Indeed, in this case R (A) docs not belong to the closed space 
L2(0, 1), the operators A and A^^ become, respectively, completely 
continuos and unbounded. As a consequence, the procedure of a nu- 
merical realization of equation (2.6), in fact, turn out to be beyond the 
sphere of the application of Banach's fundamental theorem on the in- 
verse operator. Isn't it a too high price to pay for seemingly ephemeral 
clarity in the formulation of the problem under the conditions of map- 
ping inside the space L2? 

We draw attention to a known point of view that a choice of ap- 
propriate spaces for the solutions to problems of mathematical physics 
should be done on the basis of practical applications, which can hardly 
be disputed. As the same time, a wide-spread opinion that, for example, 
a sociologist should formulate a problem to be solved by mathematical 
methods with a specification of appropriate spaces for its data. This, 
as a rule, admits variety, which is a prerequisite for an increase in the 
efficiency of procedures of numerical realization. 

Are there any prospects to overcome the above-mentioned complex- 
ity in mating the free term of Eq. (2.6) with the adequate space In 
this regard, let us turn to Eq. (2.10) that by virtue of (2.11) takes the 
form 



However, there is a chance of a reduction of the given function / {x) 
to /* e ^2 by means of adaptive simulations of the error {6f') (x) . In- 
deed, it can be interpreted as the smoothing of information by the 
procedure of integration. From this point of view, it seems to be rea- 
sonable to represent 6f' as a difference between the explicit form of the 
sought function ip (x) and an integral over this function whose kernel 
would not impose any additional restrictions on the formulation of the 
problem. Simultaneously, given that the error of integration by means 
of (2.6) is objectively small, there appears a condition of the form 



(Ai:) (x) = f(x) + (Sf) (x), xe [0, 1] . 



(2.12) 



II^/IIl2(o,i) 



0. 



(2.13) 
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Thus, instead of traditional determination of the function ip (x) by 
means of the solution of the Fredholm equation of the first kind (2.6), 
we suggest to employ the perturbation 5f' which leads to the prob- 
lem (2.12), (2.13). In this way, a prerequisite is formed for ensuring 
f + 5f' G R (A) 7 As will be shown, by use of some additional consider- 
ations, the determination of the function tjj (x) that satisfies Eq. (2.6) 
can be reduced to the solution of a well-posed problem. 

Note that, in the case of a considerable mismatch between R{A) 
and the actually known function / (x), condition (2.13) can hardly be 
regarded as feasible. Nevertheless, the outlined approach still applies 
interpreting, figuratively, the reduction of the free term of Eq. (2.6) to 
a form which makes it solvable. 



'^A practical realization of the outlined orientation is a key aspect of the con- 
structive part of the present consideration (see sections 4.5, 5, 6). 
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Chapter 3 



The existing approaches to 
the solution of ill-posed 
problems 

3.1 A. N. Tikhonov's methodology 

The consideration of this subsection is based on the material of the 
monograph by A. N. Tikhonov and V. Y. Arsenin [1] that is, hterally, 
pierced by the concept of the adequacy of incorrect formulations and, in 
particular, of integral equations of the first kind to problems of math- 
ematical physics. As an illustration, we show that the solution to the 
Predholm integral equation of the first kind 

b 

{A^) {x) = jk (x, ^ (0 di-fix), xe [a, b] , (3.1) 

a 

with k {x, ^) and dxk (x,^) being continuous with respect to x, can un- 
dergo arbitrarily considerable changes both in the metric C and L2 for 
small in L2 (a, b) variations of the right-hand side in the form e sin (ujC,)- 
The situation with the perturbation of the kernel k {x, ^) is, in fact, 
analogous. In this regard, the authors pose the question: What should 
be understood by the solution of Eq. (3.1) when k and / are known 
approximately? In their opinion, a problem of this type should be 
considered "underspecified" , and, correspondingly, a choice of possible 
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solutions should be made taking into account "usually available" addi- 
tional qualitative or quantitative information about the function ijj (x). 
In this regard, we draw attention to N. G. Preobrazhenskii's consider- 
ations concerning a system of linear algebraic equations, obtained by 
the discretization of (3.1) ([2], p. 130): 

"An analysis shows that choosing sufficiently high order of an ap- 
proximation, we transform [the above-mentioned problem] into an ar- 
bitrarily ill-defined one... Under these conditions, it is necessary to add 
to the algorithm some a priori nontrivial information, only by the use 
of which we can expect to filter out veiling false variants and single 
out the solution, closest to the sought one. any purely mathematical 
tricks that do not employ additional a priori data are equivalent to 
an attempt to construct an informational perpetuum mobile producing 
information from nothing." 

The so-called method of the selection of the solution to ill-posed 
problems is based on a priori quantitative information. It is shown 
that if a compactum M of a metric space Ei is mapped in a one-to-one 
and continuous manner onto a set F of a metric space E2, the inverse 
map F onto M is also continuous. Correspondingly, an assumption 
that the solution, in particular, to Eq. (3.1) belongs to the compactum 
M allows us to consider the operator to be continuous on the set 
F = AM. 

A practical realization is reduced to an approximation of M by 
a series with parameters that change within certain limits (for M to 
represent a closed set of a finite dimensional space) and should be de- 
termined from the condition of the minimum of the error of closure of 
(3.1). Note the absence of any more or less general recommendation 
with respect to the choice of M. 

In light of the above, M. M. Lavrentiev has formulated the notion 
of correctness according to Tikhonov for an equation of the type (3.1) 
[3]: 

1) It is a priori known that the solution ip^, to the considered equation 
exists and belongs to a set M of the space Bi. 

2) The solution t/j^ is unique on the set M. 

3) The operator A"^ is continuous on the set AM of the space B2. 
If M is a compactum (this case is called " usual" ) the last condition 

becomes a consequence of the first two conditions. 
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Those problems in which the operator is unbounded on the 
set AEi and the set of possible solutions Ei is not a compactum are 
called substantially ill-posed. For such problems, Tikhonov has put 
forward an idea of a regularizing operator G, in a sense close to A~^, 
whose value domain for the map from E2 into Ei admits matching to 
the right-hand side of (3.1), known approximately. Moreover, G must 
contain a regularization parameter a that depends on the accuracy of 
the initial information. 

The operator G (/, a) is called a regularizing operator for Eq. (3.1) 
if it possesses the following properties: 

1) It is defined for any a > and / G £'2. 

2) For Atp^: = /^,, where i/j^, and are corresponding exact expres- 
sions, there exists such a (6) that for any < e < p^j (V',,'V'a) there is 
S (e) > pE2 (^. , ^a) • Here, iIj^ = G {f,a). 

It is implied that there is a possibility of a choice of a {5) such that 
for 5 — the regularized solution '0^ i/;*, i.e., e — > 0. At the same 
time, it is pointed out that the construction of the dependence a (5), for 
which the operator G (/, a {5)) is a regularizing one, is algorithmically 
complicated for classes of practically important problems. There are a 
lot of pubhcations of Tikhonov's followers devoted to the resolution of 
this difficulties, which will be discussed below. 

Namely in [1], the construction of G{f,a) is carried out by the 
use of techniques of calculus of variations that reduce the evaluation of 
ijj (x) to the minimization of the functional 

[/, V'] = pk (^V', /)' + »^ m ■ (3.2) 

For Eq. (3.1), its stabilizing component is recommended to be taken 
in the form 

b 

^[^]^ j {Po {x) 1/;^ {x) +pi {x) [d^ij; {x)f} dx, (3.3) 

a 

where po, pi > are given functions. 

In the case of a symmetric kernel k {x, C,), the procedure of the min- 
imization of (3.2) is equivalent to the solution of the integrodifferential 
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equation 

a {po (x) i/j (x) - [pi (x) d^i/j (x)]} + (Ai/j) (x) = / (x) , x e [a, b] , 

(3.4) 

under the conditions 

p,{x)d,iPix)v{x)\''^^0. (3.5) 

Here, v (x) is an arbitrary variation of ip {x) in the class of admissible 
functions. 

In the opinion of the authors of [4], an overwhelming majority of 
inverse problems are ill-posed, and attempts to solve them, in view 
of their great practical importance, were being undertaken for a long 
period. "But only as a result... of the appearance of fundamental 
publications of academician A. N. Tikhonov, the modern theory of the 
solution of inverse problems, based on the notion of a regularizing algo- 
rithm, was constructed" (p. 7). In what follows, the authors construct 
the procedure of a numerical realization of the Fredholm integral equa- 
tions of the first kind, related to the interpretation of astrophysical 
observations, by means of the selection of the compactum of possible 
solutions in the class of monotonically bounded functions. 

As is pointed out by O. A. Liskovets [5], "...the correctness accord- 
ing to Tikhonov is achieved at the expense of the reduction of the 
admissible manifold of solutions to the class of correctness" (p. 13). 
The following quotation from the above-mentioned monograph is also 
of interest: "In contrast to a previously prevailing opinion that all the 
problems describing physical reality are ill-posed, according to the mod- 
ern point of view any realistic problem can be regularized, i.e., it has 
at least one regularizer" (p. 14). 

Here is V. A. Morozov's conclusion ([6], p. 9): "A. N. Tikhonov's 
method of regularization turned out to be simple in practice, because 
it did not require actual knowledge of the compactum M that con- 
tained the sought solution to Eq. (3.1)... The main difficulty of the 
application of this method consists in the formulation of algorithmic 
principles of the selection of the parameter of regularization a". Ac- 
cording to his own monograph ([7], p. 4), "The importance of A. N. 
Tikhonov's paper [8] can hardly be overestimated. It served as impetus 
for a number of publications by other investigators in different fields of 
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mathematical analysis and natural science: spectroscopy, electron mi- 
croscopy, identification and automatic regulation, gravimetry, optics, 
nuclear physics, plasma physics, meteorology, automation of scientific 
research and some other spheres of science and engineering" . 

V. V. Voevodin's opinion ([9], p. 43) is as follows: "The success of 
the application of the regularization method to the solutions of unstable 
systems of algebraic equations is explained to a large extent by the fact 
that A. N. Tikhonov and his followers did not restrict themselves to 
an investigation of separate fragments of this complicated problem but 
considered the whole complex related issues. This, in the first place, 
concerns a clear formulation of the problem itself, the construction of 
a stable with respect to perturbation of the input data algorithm of its 
solution, the development of an efficient numerical method, estimates of 
a deviation of the actually evaluated object from the sought one taking 
into account a perturbation of the input data and errors of rounding" . 

A quotation from the preface to the collected volume by A. A. 
Samarsky and A. G. Sveshnikov [9] reads: "A clarification of Audrey 
Nikolaevich Tikhonov of the role of ill-posed problems in classic mathe- 
matics and its applications (inverse problems) is of fundamental impor- 
tance for the who;e modern mathematics. He proposed a principally 
new approach to this class of problems and developed methods of the 
construction of their stable solutions based on the principle of regular- 
ization" . 

3.2 A brief review of the development of 
the outhned concepts 

The results of investigations devoted to the determination of the regu- 
larization parameter a are summarized in [10]. Based on the assump- 
tion that errors in the determination of the free term / (x) and the 
kernel k {x, ^) of Eq. (3.1) are known, one uses different methods of the 
minimization of the error of closure of the type 

Aip^- f\\^ = li\\df\\^, fie (0,1). 

The evaluation of the parameter a as a root of the corresponding equa- 
tion does not pose any problem. However, a choice of n is, in fact. 
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related to considerable uncertainty. The main obstacle is that a reliable 
estimate of the error stipulated by the " measure of incompatibility" of 
the concretely considered equation Aip = / is rather questionable. 

Considerable efforts were undertaken to reduce the volume of infor- 
mation necessary for the evaluation of the parameter a. A noticeable 
step in this direction was made by A. N. Tikhonov and V. B. Glasko who 
suggested a criterion of the minimization of the functional ||Q;(i-?/'Q,/(iQ;|| 
with respect to a > [11] (see also [1], section 2.7). However, its theo- 
retical justification proved to be possible only for rather narrow classes 
of problems. A number of methods of the determination of a is related 
to the use of solutions to Eq. (3.1) for a special form of the functions 
fix). 

In [10], the status of the studies of estimates of the accuracy of 
methods of the solution of the integral equation (3.1) is also illumi- 
nated. If ip (x) belongs to a compactum, any serious complications, as 
a rule, do not arise, and the main interest is focused on the algorithm 
of regularization. If pi = in (3.3) and the parameter a is finite, 
Eq. (3.4) becomes a Fredholm integral equation of the second kind, 
to which, under the assumption that the error in the determination of 
k{x,^) and f (x) is known, the whole general theory of approximate 
methods of L. V. Kantorovich applies ([12], section 14.1).^ 

At the same time, as shown by V. A.Vinokurov [13], when a priori 
information about the solution to Eq. (3.1) is missing, the estimate 
of the error of the evaluation of (x) by means of regularization is 
impossible in principle. Justified is only a formulation of the question 
of the convergence of the procedure of computation or of a possibility 
of the regularization of the corresponding problem. 

In this regard, we note the arguments of A. B. Bakushinskii and A. 
V. Goncharskii ([14], p. 13): "Unfortunately, in the general case, it is 
impossible to estimate the measure of closeness of G (/, a) to A~^ (/*) 
without additional information about the solution to Eq. (3.1). This 
is a characteristic feature of ill-posed problems. In the general case, 
a regularization algorithm ensures only asymptotic convergence of an 



^Note that the definition of the type of the equation, i.e. "the second kind", in 
this case, because of the presence of a, is purely formal. This important issue will 
be repeatedly discussed in what follows. 
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approximate solution to the exact one for 5^0". 

The name of M. M. Lavrentiev is associated with a particular case 
of a practical reahzation of A. N. Tikhonov's method consisting in the 
reduction of the problem (3.4), (3.5) to the solution of the Predholm 
integral equation of the second kind 



where a > is a small parameter. 

It is shown that || V'a ~ || — >■ for 5 — >• 0, 7 — >• and (5 + 7) /a {6, 7) 
0. Here, 7 is an error in the determination of the kernel k {x,^), anal- 
ogous to 6 (see section 3.1). 

V. K. Ivanov's method [15] allows one to find the so-called quasiso- 
lution minimizing the error of closure of (3.1) for a class of functions 
ip (x) G Mr, where Mr G i?i is a compactum. The quasisolution to 
(3.1) on such a compactum has the form 



b 




X G [a, b] , 



(3.6) 



a 



00 



V' (2^) = X] c, 



■n 



(A + Xn) ipn (x), xe [a, b] . 



(3.7) 



n=l 



Here 



b 




a 



Xn and ipn (x) are, respectively, characteristic numbers and cigcnfunc- 
tions of the kernel k {x, ^); the parameter A = and represents a posi- 
tive root of the equation 




(3.8) 



under the conditions, respectively. 



00 00 



E«<^'; E«>^'- (3-9) 



n=l n=l 



Special methods of regularization are developed for the situations 
when considerable volume of information of statistical character (spec- 
tral densities, mathematical expectations, etc.) about the solution to 
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an equation of the type (3.1) is available. Thus, V. N. Vapnik [16] 
rather constructively employed the specifics of problems concerned with 
recognition of images, related to nonuniqueness and, as a result, to ex- 
treme behavior of the sought functions. We point out a definition in 
the above-mentioned monograph (p. 8) that, apparently, was implied 
by many authors but did not receive such a clear formulation: 

" The problem of the restoration of dependencies from empirical data 
was and, probably, will always be central in applied analysis. This 
problem is nothing but mathematical interpretation of one of the main 
problems of natural science: How to find the existent regularity from 
random facts." 



3.3 V. M. Pridman's approach 

Let k{x,^) be symmetric, positive definite kernel and Eq. (3.1) be 
solvable. Then, as shown by V. M. Pridman [17], a sequence of functions 
determined by iteration 



fix)- lk(x,0i^n(0d^ 



n = 0, 1,..., 



(3.10) 

converges in L2 (a, b) to the solution of Eq. (3.1) for an arbitrary choice 
of the initial approximation (x) G L2 {a, b) and < A < 2Ai, where 
Ai is the smallest characteristic number of the kernel k (x,^). 

M. A. Krasnoselskii [18] extended this result to an arbitrary solvable 
equation of the type (3.1) with a linear bounded operator A in a Hilbert 
space H. A theorem on the convergence of successive approximations 

il;n+i = {I -iyAi)iJn + i^fi (3.11) 

to the solution is proved. Here, Ai = A* A; fi (x) = {A*f ) (x); I is the 
identity operator; A* is the conjugate operator to A; < z/ < 2/ \\Ai\\; 

ipo {x) e H. 

Note that in the case of the integral operator (3.1) 
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where 

b 

a 

A number of procedures are known that improve convergence of 
iterations according to Fridman (see [10]). For example, under the 
conditions that are specified with respect to the procedure (3.10), 

V'n+1 (x) = — - ^ ipn (x) , (3.12) 

^ + n=0 

where (fo (x) e L2 (a, b); 

b 

(fin (x) = (fin-l (x) + f{x)- j k (X, ^) (pn-1 (0 71 = 1, 2, ... . 

a 

G. M. Vainikko and A. Y. Veretennikov [19] studied an iteration 
algorithm of an implicit type: 

b 

aipn+i {x)+ j k {x, (0 = oitpn (x) + f (x) , n = 0, 1, . . . , 

(3.13) 

where ■0o (x) & L2 (a, b); the parameter a> 0. 

Note that in contrast to the regularization of the type (3.6), based on 
the smallness of a, the considered approach is characterized by multiple 
iteration with, on the contrary, sufficiently large value of this param- 
eter. Moreover, one of the merits of the procedures (3.10)-(3.13) is a 
possibility of a constructive application of an a posteriori estimate of 
the error to accomplish the iteration. 

In the simplest case, one finds the number n for which for the first 
time 

llV'n+l - V'n|L,(a,6) < Ci5 + C27, 

where 6 and 7 are errors in the determination of f (x) and k{x,^), 
respectively; ci, C2 are constants meeting a number of requirements to 
ensure the stability of the procedures of computation. The infiuence of 
errors, small in a probabilistic sense, on the convergence of successive 
approximations is also investigated. 
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The authors of [20] gave arguments for usefulness of the combination 
of the regularization of the equation of the type (3.1), whose parameter 
is the number of iterations, with algorithms of the saddle-point type. 
This approach has its origin in the pubhcation by V. M. Pridman [21] 
and is realized, in particular, according to the scheme 

ijjn+l ='K- PnA* {A^n ' f) , (3.14) 

where 

_ \\A*{A^^-f)f 

^" \\AA*iA^r.-f)\r 
which is adequate to the choice of the step of the descent from the 
condition of the minimum of the error of closure 

A„+i = [[AV'n+l - /|lL2(a,&) ■ 

3.4 Inverse problems for differential equa- 
tions of mathematical physics 

The monograph by O. M. Alifanov, E. A. Artyukhin and S. V. Rumyant- 
sev [20] reflects established approaches in this field. In the procedure of 
mathematical formulation of the problems, structural and parametric 
identification is emphasized, which implies, respectively, a qualitative 
description of the considered processes by means of differential opera- 
tors and allotting quantitative information to the model. 

Interpretation of physical processes in terms of causality is also 
given. The cause includes boundary and initial conditions with their 
parameters, coefficients of the differential equations and also the do- 
main of the problem. The effect reflects the status of the investigated 
object and represents, mostly, fields of physical quantities of different 
types. 

The restoration of the cause from the information about physical 
fields is considered as an inverse problem. A key consideration is as 
follows (p. 11): "A violation of a natural causal relation that takes 
place in the formulation of the inverse problem can lead to its mathe- 
matical incorrectness, such as, in most cases, instability of the solution. 
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Therefore, inverse problems constitute a typical example of ill-posed 
problems" . 

In connection with the sought function, the following types of in- 
verse problems of the identification of physical processes for partial 
differential equations are singled out: 

1) Retrospective problems: the determination of the prehistory of a 
certain state of the problem. 

2) Boundary problems: the restoration of boundary conditions or 
of the parameters contained therein. 

3) Coefficient problems: the restoration of the coefficients of the 
equations. 

4) Geometrical problems: the determination of geometrical charac- 
teristics of the contour of the domain or of the coordinates of points 
inside. 

A principal difference between inverse problems of identification and 
those of regulation is pointed out, concerning the width of classes of 
possible solutions. Whereas in the former case their increase leads to 
complications in the numerical realization, in the latter case, on the 
contrary, this is a favorable factor. By the way, the algorithmic means 
[20] are almost completely based on the methods of the solution of 
integral equations of the first kind, to which the considered problems 
of heat exchange are reduced. 

In the formulation of inverse problems of mathematical physics, the 
proof of corresponding theorems of existence and uniqueness is of pri- 
mary importance. In this regard, a general approach, outlined schemat- 
ically by A. L. Buchgeim ([22], pp. 133-134) can be mentioned. Thus, 
the following equations are considered: 

Pu^f; Qf^g, (3.15) 

where P is an operator of the direct problem; Q is an "information" 
operator describing the law of the change of the right-hand side; g 
is given, whereas u and / are the sought elements of corresponding 
functional spaces. 

The application of the operator Q to the first equation (3.15) yields 
QPu = g, which is equivalent to 



PQu = [P, Q]u + g, 
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where [P, Q] = PQ — QP is the commutator of the operators P and Q. 

The meaning of the commutation hes in the fact that, as a rule, there 
is no information, except for (3.15), about the function /. Therefore, 
it is easier to study the operator on the solution of the direct problem 
u that satisfies some manifold of boundary conditions. It is important 
that in typical applications the operator Q does not "spoil" the part of 
boundary conditions that reflects the domain of the operator P. As a 
result, one gets a specific factorization of the inverse problem (3.15) as 
a product of two direct problems, induced by the operators P and Q 
under the condition that the commutator is, in a sense, "subordinate" 
to them. 

In the trivial case [P, Q] = 0, the initial problem decomposes into 
two simpler ones: Pv — g; Qu — v. For the description of properties of 
the employed operators, a priori estimates are used. 

Of interest is also a quotation from the introduction to the mono- 
graph by R. Lattes and J.-L. Lions [23]: "In this book, we suggest a 
method of quasiinversion, intended for the numerical solution of some 
classes of ill-posed, according to Hadamard, boundary value problems. 
Practical and theoretical importance of such problems is being more 
and more realized by investigators". And further: "The main idea of 
the method of quasiinversion (universal in numerical analysis!) con- 
sists in an appropriate change of operators entering the problem. This 
change is done by the introduction of additional differential terms that 
are 

i) sufficiently "small" (they can be set equal to zero); 

ii) "degenerate on the boundary" (to prevent, for example, the ap- 
pearance of complicated boundary conditions and of such conditions 
that may contain unknown, sought variables)". 

In particular, the ill-posed problem of thermal conductivity 



dtu - dl 



= 0, 



(3.16) 



ii(0,i)=ii(l,i) = 0; u{x,T)^C{x), 

where C, {x) is an unknown function, is replaced by the following, with 
a small parameter e: 



dtU - dl 



u - edl 



:u = 0; 



(3.17) 
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u — d^u — 0, X — 0; X — 1; u {x,T) — ( {x) . 

The authors point out (p. 36): "In a numerical reahzation, it is 
natural to choose e as the smallest possible one. However, in problems of 
the considered type, one should expect numerical instability for e ^ 0. 
Therefore one can expect at most that for any problem there exists a 
certain optimal value of e equal to eo". The absence of convergence "in 
a usual sense" of the solution of the problem (3.17) to the exact one 
for e — > was pointed out by A. N. Tikhonov and V. Y. Arsenin ([1], 
p. 52). 

3.5 Alternative viewpoints and develop- 
ments 

In Y. I. Liubich's opinion, any more or less general theory of integral 
equations of the first kind is absent, and only in some cases it is possible 
to use special methods. An example is given by known Abel's equation 
([24], p. 83). 

K. I. Babenko's remark ([25], p. 310) is rather typical: "Although 
from the point of view of the loss of information algorithms are not 
estimated, it seems to us that this is an important characteristic and it 
should be taken into account" . In what follows, the lack of optimality of 
the traditional approach to a numerical realization of ill-posed problems 
is concretely demonstrated. 

A profound analysis of methodological aspects of this sphere is given 
by R. P. Fedorenko ([26], sections 40, 41). In particular, he failed to 
establish the value of the regularization parameter a by minimizing the 
functional (3.2), because for small values the sought function began to 
oscillate, whereas with it increase the value of considerably exceeded 
the admissible one. The author arrived at the conclusion that reason lay 
in the inadequacy of the theory [1] to problems of control, characterized 
by discontinuity of solutions. 

By studying the problem (3.16), Fedorenko brought up the following 
consideration: "All the methods of the solution of ill-posed problems 
more or less consist in preventing the appearance in the sought solution 
of higher harmonics with large or even simply finite coefficients. But 
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what is "high frequency"? Beginning with what number n should we 
consider the function sin (mrx) redundant, only spoiling the solution? 
This, of cause, depends on T". It is implied that a hypothetically 
known solution of the corresponding direct problem can be expanded 
into a Fourier series 

oo 

u{x) = u (x, 0) = ^ sin (nTTx) . 

n=l 

It is shown that the use in [23] of the value T = 0.1 and the errors in 
L2 (0, 1) of the satisfaction of the last condition (3.16), with 5 of order 
10~^, imposes the restriction n — 2. In this context, the method of 
P. Lattes and G.-L. Lions came under criticism. These authors, while 
solving the problem (3.17) on a grid with a step of Ax = 0.02, obtained 
an absolutely unacceptable component Uq, namely, 10^ sin (Gvra;). This 
occurred for 6 at the level of 0.05, under the conditions when \( {x)\ < 
1... . 

Note also the remark ([26], p. 360) that, aside the fact of the bound- 
edness of the regularizing operator G (see section 3.1), its norm 
is an exceptionally important characteristic whose value directly influ- 
ences a relation between the accuracy of the given function ( and the 
solution uo — G(.'^ 

Indeed, let us consider Eq. (3.6), written in the canonical form 

1 } 1 
i^ix)^— kix,Oi^iOdC+-fix), xe[a,b]. (3.18) 
a J a 

a 

Let a = 0, 6 = 1, and the kernel k{x,^) be determined by the 
expression (2.8). In this case, for a"^ ^ in-n)'^, its solution is [27]: 

i/j (x) = - f (x) - -Y, — — 7 — ^2 sin (^TTx) , 

1 

Cn = J f (x) sin (mrx) dx. 




^By the way, in most specialized publication this issue is not accentuated. 
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It is not difficult to notice that for small values of a the error in the 
determination of the function / (x) can considerably distort iJj (x) [see 
also a footnote concerning the solution of Eq. (3.18) in section 3.2]. 

In a constructive aspect, Fedorenko recommends to use traditional 
formulations of inverse problems of differential or variational charac- 
ter with an application of additional conditions that rationally restrict 
classes of possible solutions. As the main factor to achieve the de- 
sired efficiency, a comprehensive analysis of qualitative peculiarities of 
solutions to the considered problems, involving elements of numerical 
simulations, is suggested. 

What are the values of the regularization parameter a, typical of 
computational practice? The authors of [28] point out that for problems 
of restoration of time-dependent density of thermal flux on the surface 
from the results of temperature measurements at internal points of the 
samples the corresponding range is rather representative: 10~^ — 10~^. 
The editors of the above-mentioned book have a different point of view: 
" One can give a lot of examples of solutions to inverse problems thermal 
conductivity, when the range of acceptable values of a turn out to be 
rather narrow" (p. 141). 

The main technique of a numerical realization [28] is interpreted 
by the authors as a complement to the method of least squares by 
a procedure that smooths oscillations of the solutions in high order 
approximations. In this regard, they point out a relationship between 
Tikhonov's regularization and algorithms of singular expansions and 
ridge regression (or damping) that are widely used for the suppression 
of the instability of the method of least squares [29] . 

In a number of publications, one can see an orientation towards reg- 
ularization of Eq. (3.1) without the distortion of the original operator 
along the hues of (3.4) or (3.6). Thus, A. P. Petrov [30] suggested a 
formulation of the problem with / (x) G i? (A) by means of the repre- 
sentation / = Alp + UJ, where a) is a random process reflecting errors of 
the data and of the calculations. At the same time, the author failed 
to use his formally achieved correctness to construct an efficient algo- 
rithm of a numerical realization. It seems that the reason lies in the 
insufficiency of the structure of a; from the point of view of adaptive 
compensation of the error of closure of the satisfaction of (3.1). 

A. V. Khovanskii [31] put forward arguments for the regularization 
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of the algorithm of the solution of Eq. (3.1), not the operator A (which 
is the basis of the theory of [1]). The following quotation is of interest: 
"What is more, Tikhonov's regularization contains in an inseparable 
form two completely different notions, accuracy and stability, and there 
is a transformation of one into another. Nevertheless, there exists for a 
long time an idea of the predetermination of the operator [32] , although 
only in the context of conjugate gradients and in a multiplicative form" . 

However, the method of conjugate gradients is, in fact, Fridman's 
iterations of the type (3.14). Note that nonlinearities contained therein 
facilitate the smoothing of a well-known slow-down of the convergence 
of the procedure (3.10) with approaching the solution to Eq. (3.1). 
This effect was demonstrated by A. D. Myshkis [33] with the help of 
the representation of the components of (3.10) by series in terms of the 
eigenfunctions of the kernel k{x,^). This leads to the relations 

Cn+l,m = (1 ~ '^Z '^m) Cn,m + '^fm, TTl — 1,2, . . . , 

where Cn,m and fm are coefficients of the above-mentioned expansion of 
ipn (x) and / (x), respectively. 

When the number of the terms in the representation of the solution 
increases, which seemingly had to improve the accuracy, the coefficient 
of convergence 1 — X/Xm approaches unity and, as a result of the accu- 
mulation of errors, the iterations become "counterproductive". 

Note an effective method of the suppression of instability of the 
procedure of a numerical realization of the Predholm integral equation 
of the second kind 

b 

i;{x) = xJk{x,Oi^{OdC + f{x), xe[a,b], (3.19) 

a 

"positioned on the spectrum", i.e., in the case when A = A„ with 
being a characteristic number, proposed by P. I. Perlin ([34], pp. 105- 
107). 

This problem is ill-posed both with respect to the uniqueness of 

the solution and as a result of the degeneracy of the system of linear 
algebraic equations obtained by discretization. Nevertheless, a pertur- 
bation of the right-hand side of / {x) by a zero (within the limits of the 
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accuracy of calculations) component 



where ip'^ (x) is a normalized eigenfunction of the kernel conjugate to 
k x), allows one to improve radically the situation. 

The essence lies in the fact that, theoretically, the solution to Eq. 
(3.19) is expanded in a power series of A. Provided that computational 
procedures, matching this situation, are identical, one can compensate 
for the errors. 

3.6 A comparison between the main con- 
cepts of A. N. Tikhonov and V. M. 
Fridman 

A. N. Tikhonov's original suggestion (1943) admitting of the consid- 
eration of ill-posed problems by an a priori restriction on the class of 
possible solutions is a kind of refraction of general methodology of in- 
vestigations of the issues of existence and uniqueness into the sphere 
of numerical analysis. Note that A. N. Tikhonov's proof of the well- 
known theorem on the uniqueness of the solution of the inverse problem 
of thermal conductivity in an infinite n-dimensional domain under an 
additional condition of the type \d^u\ < M dates back to 1935. A vivid 
illustration of these considerations is provided by the algorithm of the 
search for a quasisolution (3.7)-(3.9). 

Behind A. N. Tikhonov's method of regularization (1963), there is 
a global idea of a limiting transition to the exact solution with respect 
to a small parameter of the problem, which is unambiguously pointed 
out in ([1], p. 56): "Note that regularizing operators, dependent on a 
parameter, have been employed in mathematics since Newton's times. 
Thus, the classic problem of an approximate calculation of the deriva- 
tive u' (x) by means of approximate (in the metric C) values u (x) can 
be solved with the help of the operator 



b 




a 




u{x + a) — u (a) 



a 
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Then, instead of the exact value of the function u (x) , an approx- 
imate one us{x) = u{x) + Au{x) with |AM(a;)| < 5 is substituted. 
On the basis of these calculations, one makes the statement: "If a = 
S/rj (S), where rj (5) ^ for S ^ 0, then {25/a) = 2r] (S) ^ for 5 ^ 0. 
Thus, for a = ai (S) = S/r] (S), G {us, ai (S)) u' (x)" . 

It should be noted that, using the methodology of a small parameter, 
Tikhonov obtained fundamental results in the field of investigations of 
differential equations with a singular perturbation of the type 

eii^ f {u,v,t); v^g{%v,t), 

where e is a small parameter; f {u,v,t) is a nonlinear function (1948- 
1952)=^. 

The solution of the system of equations does not depend continu- 
ously on the parameter e. Proceeding to the limit e ^ creates a new 
object of investigations with completely different properties. In the first 
place, it imphes the issue of the so-called violation of the stability of 
the root of the equation / {u,v,t) = 0. Nevertheless, Tikhonov man- 
aged to develop a rather constructive theory that served as a basis for 
a number of productive approaches of both fundamental and applied 
character. The importance of Tikhonov's achievements in the sphere 
of system analysis is analyzed in detail by N. N. Moiseev ([36], section 
5). 

However, properties of the integral equation (3.6) for a = also 
change radically. In this regard, generally speaking, a certain analogy 
emerges. One can suggest that Tikhonov undertook an attempt to use 
the techniques of his theory of singular perturbations for the solution 
of ill-posed problems. 

This suggestion is supported by the following quotation from the 
monograph by S. A. Lomov ([37], p. 12): "Now it is becoming clear how 
to isolate in singularly perturbed difi'erential equations small terms that 
can be neglected. It turned out that one needed additional information 
about the solution to do this." 

Note J. Hadamard's remark that an extension of methods of the 
theory of ordinary differential equations to problems of mathematical 
physics should be done with great care ([38], p. 38). At the same 



^See the review by A. B. Vasileva [35]. 
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time, at the turn of the 1950s, the theory of singular perturbations 
became an efficient tool in investigations of complicated problems of 
partial differential equations (publications by M. I. Vishik and L. A. 
Liusternik, O.A. Olejnik, K. O. Fiedrichs, and others). By the way, 
explaining the conceptual basis of their method of quasiinversion, R. 
Lattes and G.-L. Lions ([23], p. 11)^ refer to these authors and A. N. 
Tikhonov. 

Simultaneously, they pointed out that Tikhonov's priority publica- 
tion on the method of regularization [8] (see also [39]) was preceded 
by D. L. Phillips' article [40], whose results with respect to integral 
equations were analogous. In the monograph by F. Natterer [41] this 
regularization figures as Tikhonov-Phillips' method. V. A. Morozov 
estimated the achievements of the latter author in a much more re- 
strained manner ([6], p. 10): "Some recommendations on the use of 
this method are contained in the publications by L. V. Kantorovich 
[42] and D. L. Phillips [40]. There is no theoretical justification of this 
approach in the above-mentioned publications". 

The chronological reference to the most important results in the 
field of the construction of stable algorithms for the solution of integral 
equations of the first kind ([10], p. 234) gives the following informa- 
tion: "1962, Phillips's publication [40], where he suggested a variational 
method of conditional minimization of the functional (with the use of 
restrictions on the smoothness of the solution) and put forward the idea 
... of a choice of the regularization parameter a". 

Turning to V. M. Fridman's achievements, note that it is rather 
difficult to evaluate the premises that form the basis of the iteration 
procedure (3.10). At the first sight, such a computational method has 
a lot of analogs. However, its adequacy, in a sense, to the object of in- 
vestigation, the ill-posed problem of the solution the Fredholm integral 
equation of the first kind, turned out to be rather unexpected. 

Later on, with the aim to improve convergence, Fridman also em- 
ployed the nonlinear algorithm (3.14). In our opinion, different ways 
of the determination of the number of the final iteration and of the in- 
crease of the rate of global convergence (see [14, 19, 20]), despite their 



^Ideological closeness of quasiinversion and Tikhonov regularisation was pointed 
out by M. M. Lavrentiev [23, p. 5]. 
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actuality for practical application, should be interpreted as a technical 
complement to Fridman's methodology. 

Nowadays, the algorithm of conjugated gradients is considered to 
be nearly the most efficient one for the solution of large ill-defined 
sparse systems of linear algebraic equations, obtained by the reduction 
of, apparently, most problems of numerical simulations [32, 43, 44]. 
As is pointed out by J. Ortega [32], this method was proposed by M. 
P. Hestens and E. L. Stiefel (1952). However, for certain reasons, it 
was not employed for a long time. It attracted considerable interest 
at the turn of the 1970s, when one realized the actual sphere of its 
applications, the potential of the above-mentioned predetermination 
and adaptivity with respect to paralleling of computational operations 
in combination with the architecture of modern computers. 

Thus, the priority of the method of conjugated gradients ensured 
its refraction to a class of problems of linear algebra, characterized by 
the instability of the numerical realization, that is, in fact, ill-posed. In 
this regard, we emphasize that Fridman's "methods of the saddle-point 
type" [21] can be interpreted as somewhat simplified representatives of 
the family of the methods of conjugated gradients ([20], section 2.1; 
[43], section 7.1). It seems that V. M. Fridman, who was the first 
to use systematically iterations for the solution of ill-posed problems, 
essentially foresaw the development of computational mathematics that 
followed. 

In hght of the above, the position of M. A. Krasnoselskii and the 
co-authors is worth noting [18]. The role of V. M. Fridman in the devel- 
opment of the iteration procedure (3.11), which is an analog of (3.10), is 
described as follows: "A transition to the equation [ip = [I — i^Ai) ip + z//i] 
was pointed out for some cases by I. P. Natanson [45]. For Fredholm 
integral equations of the second kind, it was already employed by G. 
Wiarda [46]. For integral equations of the first kind, it was, essentially, 
employed in the publication by V. M. Fridman [17]" (p. 73). There 
is no comment on a qualitative difference between the objects of the 
investigation. 

The nontriviality of Fridman's approach is noted in the remark of 
Natanson [45]: "Our method does not apply to the solution of the 
integral equation of the first kind. This could be expected, because the 
use of the method implies complete arbitrariness of the free term of the 
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equation Aip = f, whereas Eq. (3.1) is solvable not for all In 
what follows, the author gives an extended proof of the degeneracy of 
the corresponding discrete problem. 

The gradient algorithm of V. M. Pridman [21] is mentioned by the 
authors of [18] exclusively in the context of the equation Aijj = f, 
where both the operators A and A^^ arc bounded (p. 115). We quote 
the abstract to Fridman's paper [21]: "We present a new proof of the 
convergence of methods of the saddle-point type for a linear operator 
equation. We do not assume, unhke L. V. Kantorovich [47], M. A. 
Krasnoselskii and S. G. Krejn [48], that zero is an isolated point of the 
spectrum of the operator" 

3.7 Ill-defined finite-dimensional problems 
and issues of discretization 

In this subsection, Aip = f denotes a system of linear algebraic equa- 
tions. The conditionality number of the matrix A (see, e.g., [49]) 



where ^0 is a manifold of vectors of the Euclidean space, represents a 
raising coefficient between a relative error of the data and the solution. 
At the same time, cond(A) characterizes the measure of closeness of 
A to a. degenerate matrix, for which the solution of the corresponding 
system of algebraic equations does not exist or is nonunique. 

An algorithm of the solution of a degenerate system of linear alge- 
braic equations, based on the method of least squares, is presented in 
the book by A. N. Malyshev [50]. First, the matrix A is transformed 
to a two-diagonal one by means of a special transformation, and one 
finds its eigenvalues that are subdivided into two groups, ai, (J2,..., Un 
and (Tn+i,---, such that (T„/ ((7„ — (T„+i) is not very large. Then, with 
the help of a rather laborious procedure of the exhaustion of the sec- 
ond group of the eigenvalues, one constructs a matrix An that is stably 
invertible beginning with a certain value n. The accuracy of the thus 



'This is equivalent to the boundedness of the operator A 
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obtained generalized solution ip is determined by the error of closure 
Alp — f / Alp , using heuristic considerations. 

It seems that in a methodological sense this scheme reminds of B. K. 
Ivanov's algorithm [15] that reflects computational relations (3.7)-(3.9). 

L. Hageman and D. Young [43] studied the approach of predetermi- 
nation, employed for the solution of systems of linear algebraic equa- 
tions, close to degenerate ones, to accelerate by the method of conju- 
gated gradients iterations of the type 

1pn+l = Pipn + g, 

where P = I — Q~^A; g = Q~^f. It is assumed that this procedure 
can be symmetrized in the sense that there exists a non-degenerate ma- 
trix W such that the matrix W {I — P) is symmetric and positive 
definite. 

By use of W, the initial problem can be reduced to the solution of 
much better defined systems of algebraic equations Bip = q, where 

B ^W(I - P)W-^; ip^Wjp; q^Wg. 

Formally, a choice of the predeterminer does not pose problems. 
However, in practice, one has to resolve a contradiction between the 
conditions imposed on the matrix W: "closeness" to A~^ to reduce the 
number of iterations; a "rapid" calculation of a product of the type 
Wtp [51]. In the above-mentioned publication, I. E. Kaporin analyzes 
difi'crcnt approaches to the construction of predeterminers for systems 
of linear algebraic equations of a general type. An analogous issue, in 
the interpretation of J. Ortega [32], is oriented mainly towards sparse 
matrices. 

The complexity of problems of linear algebra that arise in the re- 
alization of modern methods of investigations in the field of the me- 
chanics of a continuous medium are characterized as follows [51]: "The 
matrices of corresponding systems are rather large (up to a hundred 
thousand nonzero elements), rather densely filled (up to hundreds or 
even thousands of nonzero elements in each line), have no diagonal pre- 
dominance, are not M-matrices and are rather ill-defined. In general, 
one can expect only symmetry and positive definiteness of the matrix 
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of the system".® 

Note that, for example, in seismic tomography [44], one has to be 
satisfied with a numerical realization of discrete analogs of integral 
equations of the first kind, because their kernels cannot be represented 
analytically and parameters of the considered models are determined 
with the help of natural experiments. 

In light of the above, the considerations of R. W. Hamming ([52], 
p. 360) may seem to be archaic: "A system of linear equations is said 
to be ill-defined, if, roughly speaking, the equations are almost linearly 
dependent. Many efforts were made to investigate the problem of the 
solution of ill-defined systems. However, one may pose the question: Is 
it necessary to solve such systems in practical situations? In what phys- 
ical situation may the solutions prove to be useful, if they depend in 
such a substantial manner on the coefficients of the systems? Usually, 
the following is true: Instead of the solution, one is looking for a system 
of almost linearly independent equations. In light of this information, 
the problem can be better understood and is usually reformulated again 
in a more satisfactory way. It is rather probable that ill-defined sys- 
tems of equations, provided that round-off and measurement errors are 
eliminated, are actually linearly dependent and thus do not refiect the 
physical situation". 

Note that the renowned practitioner adheres to the position of cor- 
rectness according to Hadamard. Let us quote P. S. Outer's preface to 
[52]: "The name of R. W. Hamming, a renowned American scientist, 
former President of the Computer Association, Head of the Mathemat- 
ical Service of Bell Telephone Laboratories, and his works in the field of 
computational mathematics and the theory of information are rather 
well-known and do not need special recommendations. ... The book 
'Numerical Methods for Scientists and Engineers' is without any doubt 
an outstanding phenomenon in mathematical literature" . 

Of special interest is Hamming's opinion about the priority of com- 
putational procedures ([52], p. 90): "It is often believed that the main 
problems of numerical analysis are concentrated on interpolation, but 
this is not the case. They are mostly related to such operations as 



^The non-diagonal elements of an M-matrix are non-positive, and all the ele- 
ments of its inverse are non-negative. 
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integration, differentiation, finding zeros, maximization, etc., in those 
cases when all we have or can compute arc some nodes of functions 
that are usually known not exactly, but approximately, because they 
are spoiled by the round-off error" . 

Thus, the problem should be posed correctly despite an inevitable 
error in the data. It is obvious that such a position witnesses the 
preference of algorithmic efficiency to the quality of initial information. 
Interpolation, mentioned in the above quotation, implies approximate 
representation of the latter for the performance of computer operations 
by means of a finite-dimensional approximation. 

However, in computational mathematics, alternative concepts are 
rather wide-spread, which is reflected in K. I. Babenko's remark [25]: 
"In some spheres of numerical analysis, the theory of approximation 
serves as the foundation for the building of the numerical algorithm" 
(p. 138). "Information, inputted into the algorithm, is characterized, 
in the first place, by its volume... All other characteristics, such as, 
e.g., accuracy, are its derivatives and do not present a true picture of 
the input" (p. 281). 

Here, information is understood in the sense of Kolmogorov's the- 
ory of e-entropy that identifies it with the length of a given table or 
an alphabet, whose words are manipulated by the algorithm. Corre- 
spondingly, the issue of numerical analysis is interpreted in terms of, 
figuratively, the deficiency in the search for necessary words and of the 
deletion of tables in the course of operations. 

Nevertheless, R. W. Hamming's point of view on the relation be- 
tween the method of investigations and the employed information is 
actively developed by a group of specialists with J. Traub and G. 
Wasilkovski at the head. The authors of [53] point out (pp. 9, 6): 
"In this book, we construct a general mathematical theory of optimal 
reduction of uncertainty. We interested in the two main questions: 1) 
Is it possible to reduce uncertainty to a given level? 2) What will it 
cost? The aim of the theory of informational complexity is to provide 
a unified approach to investigations of optimal algorithms and their 
complexity for the problems that involve incomplete, imprecise or paid 
information and to employ the general theory to concrete problems 
from different fields" . 

Here, complexity implies the number of arithmetic operations, the 



3.8. THE CRISIS OF THE TECHNOLOGY OF NUMERICAL SIMULATIONS63 



time of their realization, computer memory resources, etc. By the way, 
the interpretation of the notion of information [52, 53] correlates with 
the expressive statement of R. Bellman and S. Dreyfus ([54], p. 342): 
"Fortunately, in some cases, there is a very simple way to overcome 
this difficulty. Instead of trying to study information as the "smile of 
Cheshire Cat" , we consider the actual physical process, where informa- 
tion is used to work out solutions.^ The value of information can then 
be measured by the efficiency of the solutions. 

Thus, the usefulness of information depends on its application, which 
is the most reasonable concept!" 

It should be noted that the procedure of finite-dimensional approx- 
imation of problems of mathematical physics is, of course, also very 
important, which is accentuated by Babenko. Indeed, the obtained 
discrete model can turn out to be incorrect, and the employed algo- 
rithms of the numerical realization may prove to be divergent even in 
the solution of rather ordinary problems. An example of instability of a 
finite-difference scheme is given by S. K. Godunov and V. S. Ryabenkii 
([55], section 4.9). 

Babenko also emphasized the absence of any general methods of 
the construction of finite-dimensional analogs ([25], p. 622): "...the 
provision of an approximation alone is insufficient"... one has to ensure 
that the discrete problem "retains the type of the original continuous 
problem". In his opinion, to achieve the above goal, "a detailed inves- 
tigation in each concrete case is required, which is the most nontrivial 
part of work" . 

3.8 The crisis of the technology of numer- 
ical simulations 

Of considerable interest is, in fact, a program statement of O. M. Be- 
lotserkovskii and V. V. Stchennikov in the preface to [56] : 

"A rapid development of computers, especially during the last 10-15 
years, with a special acuteness posed the problem of the construction of 

'^The smile of Cheshire Cat, according to L. Carrols "Alice in th Miracle Land", 
existed separately from this cat (editor's note to [54]). 
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a principally new technology of the solution of problems by computers. 
... Historically, the problems of numerical simulations (in this notion, 
we include the actual mathematical simulations related to a numerical 
experiment), being rather advanced already in the "precomputer" pe- 
riod and rapidly developing during the next periods, turned out to be 
the most conservative component of the modern technology of the so- 
lution of problems on the computer. Using, probably, redundant from 
the point of view of a mathematician expressiveness of the description, 
one can characterize the existent situation by two stable tendencies: 

- an increase of the complexity of mathematical models; 

- construction of rather sophisticated mathematical methods. 

Both the tendencies inevitably lead to a technological deadlock, be- 
cause they create comphcations in the solution of the problem of the 
construction of software-hardware means of the support of the whole 
technological chain. ... Without any pretension to profoundness and 
importance of the analogy, we dare say that the present situation in 
numerical simulations is similar to that in mechanics before the ap- 
pearance of main ideas and concepts of quantum mechanics" . 

In the introductory article [56] the same authors emphasize the phe- 
nomenon of the accumulation of the round-off error in the numerical 
realization of algorithms that include up to 10^^ operations and the 
absence of real means to estimate the error of solutions to, in partic- 
ular, evolution problems. In their opinion, "...the following conclusion 
is quite justified: a priori, any evolution problem for large times is nu- 
merically (or computationally) ill-posed in the sense of the absence of 
a practically important solution... 

In the case, when a priori or a posteriori information about the error 
of an approximate solution is absent, it is impossible to claim that the 
solution exists. This conclusion fairly agrees with A. N. Tikhonov's 
theorem that states that the problem with the data on the operator 
and the right-hand side has no solution in the manifold of approximate 
numbers" . 

Belotserkovskii and Stchennikov regard as constructive the idea that 
discrete models of the considered problems should be assembled with 
the aim of increasing the accuracy of information by means of special 
superposition. They also suggest to search for the solution in the class 
of function with a bounded variation, with would endow the difference 
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operator of the problem with smoothing properties. 

As is well-known, N. N. Yanenko paid considerable attention to the 
methodology of mathematical simulations (see [2]). His concept of over- 
coming the above-mentioned crisis is explained by O. M. Belotserkovskii 
([57], p. 106): 

" An investigation of finite-difference schemes, approximating differ- 
ent classes of equations of mathematical physics, led N. N. Yanenko 
to an extension of the notion of the scheme. For the first time, he 
begins to consider the finite-difference scheme as an independent ob- 
ject of the investigation, as a mathematical model, adequate to this or 
that physical model. This fundamental concept is based on profound 
understanding of the fundamentals of differential and integral calculus. 

Indeed, physical and mathematical models, described by differen- 
tial, integral or integrodifferential equations, are obtained from discrete 
models by means of averaging and passing to the limit with respect to 
certain parameters. This is the case, for example, in the model of a con- 
tinuous medium, where for a sufficiently large number of elements in 
the unit volume one comes to the notion of the continuous medium by 
averaging and passing to the limit with respect to the volume. In this 
regard, one can interpret a finite-difference scheme as an independent 
mathematical model with certain properties" . 

Note the fundamental, as it seems, considerations of Yanenko [2]: 
"The objects of modern mathematics, whose theoretical "nucleus" com- 
prises topology, geometry, algebra and functional analysis, are ideal 
logical constructions forming a certain operational system. We will call 
them ideal objects, which underlines, on the one hand, their practi- 
cal inaccessibility and, on the other hand, their excellent operational 
properties that allow one to make operations without loss of informa- 
tion. Ideal objects of mathematics are essentially infinite and require 
an infinite number of operations" (p. 12). 

"The development of the experimental foundation and the tool of 
investigations, the computer, increased interest in such objects as com- 
puter numbers, programs, finite automata. In this regard, the definition 
of mathematics as studies of the infinite, accepted in the 20th century, 
should be replaced by another one, more correctly refiecting its essence, 
i.e., as studies of the relationship between the finite and the infinite" 
(p. 18). 
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Of interest is the following extract from ([58], p. 89): 

"Let us make the following remark about the meaning of mathe- 
matically ill-posed problems. In the old literature [I. G. Petrovskii, 
Lectures on Partial Differential Equations (Fizmatgiz, Moscow, 1961) 
(In Russian)], the above-mentioned lesser value of ill-posed problems 
was even interpreted as their total senselessness. Nowadays it is ac- 
cepted that this is not the case. ... Nevertheless, the fact is, of course, 
that ill-posed problems are substantially sensitive to small errors. A 
misunderstanding of this fact may lead to paradoxes." 

We think that on the basis of the above one can come to a very 
important conclusion: In their construction of the conceptual basis of 
mathematical simulation, the leading specialists were guided by the 
concept of inapplicability of Banach's theorem on the inverse operator. 
Note that N. Dunford and J. Schwartz considered this theorem as one 
of the three principles of linear functional analysis, characterized as 
being rather fruitful ([59], p. 61).^ 

A quotation from K. Maurin's manual ([60], p. 51) reads: "This 
theorem [on the closed graph], in the last years, has gained itself a 
reputation of being the most important theorem of functional analysis, 
if this one is considered from the point of view of applications" . 

An attempt to renew the above-mentioned fundamentals in the con- 
text of the accentuation of peculiarities of computational mathematics 
was made by A. V. Chcchkin [61], who suggested a division of sections of 
mathematics into classical and non-classical ones, respectively: "arith- 
metics, mathematical analysis, algebra, geometry, probability theory, 
etc.; mathematical logic, the theory of information and statistics, the 
theory of fuzzy sets, the theory of algorithms and recursive functions, 
methods of computational mathematics, the theory of finite- difference 
schemes, the theory of cubic formulas, methods of the solution of in- 
correct problems, etc." (p. 8). As a criterion, the authors choose the 
fact of availability of absolutely complete or partial information about 
the considered objects (points, functions etc.). 

Let us quote the abstract of section ([61], p. 78): "We define 
and study a new type of mappings that generalize classical notions. 



^The other two are the principle of Unear boundedness and the Hahn-Banach 
theorem. 
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Classical mappings realize correspondence between the points of a set. 
This implies that the points are known with absolute precision. The 
new mappings, termed ultramappings, realize correspondence between 
pieces of information about points of sets. The main construction of the 
ultramappings, termed ultraoperators, allows one to obtain separate in- 
formation about the image point from separate information about the 
inverse image point. 

Ultracontinuity of ultraoperators is defined, which is a broad gener- 
alization of the notion of the stability of methods. It is found that, for 
an arbitrary base operator, one can construct an ultracontinuous oper- 
ator over it. A class of ultracontinuous operators, termed Tikhonov's 
operators, is singled out. For these operators, the base operators are 
not continuous". Furthermore, "they are related to A. N. Tikhonov's 
ideas and methods of the solution of incorrect mathematical problems" . 

Returning to the question of adequate discretization, we quote the 
abstract of the monograph by A. A. Dezin [62]: "It is devoted to the 
description of the basic structures of multidimensional analysis and 
to the consideration of internally defined discrete problems of analysis 
and mathematical physics. It implies not merely an approximation of a 
given continuous object, but the construction its analog, starting from 
the notion allowing for discrete interpretation". 

Arguments for contradiction to physical sense of differential models 
of certain classes of problems of the mechanics of a continuous medium 
are given by M. A. Zak [63]. In this regard, he developed a general 
approach, wholly based on the concepts of theoretical mechanics with 
a special interpretation of Gauss' principle of least action. 

The position of C. Truesdell is alternative. He thinks that contin- 
uum mechanics of a deformed body "is, in essence, not only subtler, 
more beautiful, majestic than a rather sparse particular case, called 
" analytical mechanics" , but it is much more suitable for the simulation 
of real bodies" ([64], p. 10). 
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Chapter 4 



Comments on the material of 
the previous sections and 
some general considerations 

4.1 The correctness of the formulation of 
problems of mathematical physics 

The conditions of correctness, formulated by J. Hadamard at the turn of 
the 20th century (see [1]) and insistently advocated by him thereafter, 
primarily attract us by their ever-increasing importance for practical 
applications. These conditions deal with the conceptual basis of nu- 
merical simulation of physically meaningful problems, which, in fact, 
is disputed by nobody. At the same time, nowadays, the prevailing 
opinion is that Hadamard's concepts are principally invalid. 

Implied is the basic statement that the properties of existence and 
uniqueness, considered by Hadamard as inherent to mathematical mod- 
els of real processes, lead to the correctness of the formulation of ade- 
quate boundary-value (initial-boundary-value) problems, which implies 
the stability of the employed algorithms of a numerical realization. A 
particular consequence is that the Predholm integral equation of the 
first kind is simply unsuitable for "application" in the problems of 
mathematical simulation. 

A natural course of investigations with the aim to confirm or dis- 
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prove the hypothesis, or, maybe, a prophecy, of Hadamard, seemingly 
had to be conducted from the position of variabihty of formulations of 
the considered problems, which was not the case. The main reason is, 
apparently, a formulation of the belief in a special mission of compu- 
tational means of numerical simulations that lightheartedly neglected 
even one of the main principles of functional analysis, i.e., Banach's 
theorem on inverse operator ([2, 3], section 9, and [4]). 

One can hardly explain the absence in special literature of a consis- 
tently introduced thesis that it is necessary to coordinate constructive 
matching of the formulation of problems of mathematical physics with 
algorithms of their numerical realization. The roots of this situation 
seem to be in systemic character of the giant computer-supply com- 
plex oriented at commercial efficiency at the expense of high costs of 
provided services. 

As a result, the alternative school of A. N. Tikhonov builds up the 
criticism of J. Hadamard ideas according to the following scheme: 

- the solution of the Fredholm integral equation of the first kind 

1 

(Ai;) {x)=Jk {x, (0 d^ = f{x), xe [0, 1] (4.1) 



is, in general, an ill-posed problem (which is undisputable); 

- integral equations of this type are adequate to a variety of real 
phenomena, which is actually supported by a rather transparent in- 
terpretation of corresponding direct problems (calculations of / from 
given k and ip). 

However, what are the grounds for the formulation of the problem, 
inverse to the calculation of / (x), by means of mechanical renaming the 
given and the sought functions in (4.1)? The fact that the procedure of 
the restoration of {x) for given / {x) and k{x,$) is computationally 
incorrect not imply any consequences. 

The reproaches to Hadamard, whose typical elements are repro- 
duced in section 2.1, can be summarized as follows: The great scientist 
slowed down the progress of science by refusing to admit that ill-posed 
problems were adequate to a variety of real processes (see [3, 4, 5]). In- 
deed, the principles formulated by Hadamard do not allow for ill-posed 
problems, but this by no means imply their invalidity. In contrast to 
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Hadamard who put forward convincing arguments in support of his 
concept and, one dares say, rehed on postulates of mathematical reli- 
gion, the "science of ill-posed problems" itself could not provide any 
argument for the very justification of its existence. 

Among supporters of studies of problems of mathematical physics 
exclusively in the correct formulation are: A. Poincare, D. Hilbert, V. 
A. Steklov, I. G. Petrovsky, I. Prigogine [6, 7, 8, 9, 10]. On the other 
hand, the role of the three absolutely independent conditions of the 
correctness (existence, uniqueness and continuous dependence on the 
data of the problem), introduced by R. Courant and D. Hilbert [11], 
can hardly be called positive. 

The potential of the fact that the third condition is a corollary of 
the previous ones could facilitate the activation of research related to 
correct formulation of problems of mathematical physics. When consid- 
ering the Predholm integral equation of the first kind (4.1), one had to 
be more careful with respect to a possibility of performing correspond- 
ing transformations involving f (x) G R{A), as opposed, figuratively, 
to a surrogate of continuous inversion with the use of the rcgularization 
parameter a. 

4.2 A relationship to the theorem on the 
inverse operator 

The above-mentioned fact that the third condition of the correctness 
has the character of a corollary results from Banach's theorem on the 
inverse operator [12] whose optimistic meaning consists in the following: 
If the solution to Eq. (4.1), with D (A) = Bi and R{A) = B2, where 
B2 are Banach spaces, exists and is unique, the inverse operator 
A~^ from B2 into Bi is bounded (see section 2.3). 

Correspondingly, the procedure of evaluation of the function 

00 

"0 (2;) = X] ^nKi^n {x) (4.2) 

n=l 

(this formula follows from the Hilbert-Schmidt theorem [13]), satisfying 
Eq. (4.1) in L2 (0, 1), must be stable with respect to small perturbations 



78CHAPTER 4. COMMENTS ON THE MATERIAL OF THE PREVIOUS SECTIONS 



of k (x, ^) and / (x) under the condition B2 = 1'2- In what follows we 
assume that such a function exists, the kernel k (x, ^) is symmetric and 
closed: we use the notation of sections 2.2 and 2.4. Thus, l'^ is a Hilbert 
space of functions normalized according to (2.7). 

It should be noted that the properties of the Fredholm integral 
equation of the first kind with a symmetric kernel can be easily extended 
to the case when k (x, ^) is an arbitrary function from the space L2 ([13], 
pp. 188-194). 

However, both fulfillment in the course of calculations and a verifi- 
cation of the condition f & I2 are practically infeasible. Therefore, such 
spaces are called "inconvenient" (sec [14, 15]). Hence we are in a prin- 
ciple dilemma as to the choice of the methodology of the investigation: 

- an urge to overcome the difficulties resulting from the use of the 
space I2 related to the boundedness of the operator A'^; 

- the loss of this property in exchange for a possibility of studying 
mathematical models in "convenient" spaces. 

With the beginning of large-scale applications of computational 
methods to mathematical investigations, the second way became dom- 
inant. 

Instructive is the dynamics of the point of view of S. G. Mikhlin, 
reflected in his courses of mathematical physics and the theory of errors 
of 1968, 1977 and 1988 [16, 17, 18]. At the beginning, the author 
considers Eq. (4.1) under the traditional assumption that the operator 
A is completely continuos. In this case, the inverse operator A~^ is 
unbounded. As a result, the problem has no solution in the usual sense 
and one has to turn to the methodology of A. N. Tikhonov. 

Later, Mikhlin drew attention to the fact that if the Fredholm in- 
tegral equation of the first kind (4.1) is interpreted from the point of 
view of a mapping from the space L2 into I2, the operator A is no longer 
completely continuos, the operator A~^ is bounded, and the problem 
of the determination of the function ip (x) becomes well-posed. Simul- 
taneously, the completeness of the conditions of correctness is restored, 
whereas the third condition was initially singled out by the author. 

Thus, the use of the pair of spaces L2 (0, 1) — I2 in a sense transfers 
the canonical ill-posed problem to the mainstream of fundamentals of 
functional analysis. Note the fact that Mikhlin did not devalue the 
importance of his arguments by reasoning in terms of "convenient - 
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inconvenient" or "bad" and "good" spaces. 

Such a position apparently incurred criticism: In his concluding 
monograph, Mikhlin somewhat irritably readdresses actual formulation 
of problems of mathematical physics to specialists in applied sciences, 
including sociologists, who are interested in their solution. Simulta- 
neously, the author has found it reasonable not to consider infinite 
dimensional models with inherent aspects of incorrectness. 

There is a well-known opinion of A. M. Lyapunov that, being posed 
in the framework of initial premises, a problem of mechanics or physics 
should be solved afterwards by means of rigorous methods. Here, im- 
plied is a problem "... that is posed completely definitively from the 
point of view of mathematics" ( [19], p. 26). In other words, this means 
a well-posed problem. 

At the same time, why not consider the procedure of the formula- 
tion of problems of mathematical physics as an additional reserve of in- 
creasing the efficiency of employed techniques of numerical realization? 
Moreover, maybe rigidly predetermined formulations of problems them- 
selves prose artificial comphcations of computational character under 
the conditions when physical considerations admit a small, in a sense, 
variation? In our opinion, the formulation of problems of mathemati- 
cal physics and the algorithm of its numerical realization are essentially 
interrelated categories. 

4.3 The methodology of the solution of ill- 
posed problems 

Ill-posed problems of mathematical physics are deceptively transparent 
from the point of view of the interpretation of considered processes. 
This is stipulated, in reality, by their adequacy to spaces that in the 
computational sense are practically infeasible. If the data of such prob- 
lems are specified in their natural classes of functions, the corresponding 
formulations loose a mathematical sense because of their insolubility. 

In such a nontrivial situation, of crucial importance is, of course, 
a role of general methodological concepts. In other words, one has to 
be guided by a certain system of global principles. From this point of 
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view, if Hadamard's insistence on the correct formulation of problems 
describing physical phenomena [1] still can be interpreted as a kind of 
hypothesis, in fact related Banach's theorem on the inverse operator is a 
universally accepted element of the foundation of modern mathematics 
[20]. 

Nevertheless, there appeared a notion of correctness according to 
Tikhonov that played up a version of a search for the solution of the 
problem (4.1) in a reduced class of functions [14]. Any general recom- 
mendations for finding such a class on the basis of reasonable informa- 
tion were not worked out. 

A shaky conceptual basis led to the failure of the idea of a limiting 
transition with respect to a small parameter in the solution of a family 
of problems that mimicked ill-posed ones (the method of regularization 
[2]). The reason, apparently, lies in the same inadequacy of the use of 
functional spaces. Given that I2 is characterized by an infinite number 
of features that depend on the operator A (a superposition of products 
of squared values that consist of characteristic numbers, integrals over 
free terms and eigenfunctions) , whereas L2 is characterized by only one 
(an integral over the squared function), is it possible, even on a purely 
heuristic basis, to expect to overcome this cardinal disagreement with 
the help of the regularization parameter a? 

The situation in the sphere of activity of numerous followers of 
Tikhonov looks rather deplorable. Actually, the efforts are concen- 
trated on a mathematical object with a small factor a, formed on the 
basis of (4.1): 

1 

ai^{x) + jk{x,i)i^{i)d^^f{x), xe[0,l]. (4.3) 


This is called the Predhohii integral cqiiation of the second kind, with- 
out any mentioning of its insufficiency in this respect. Despite a large 
number of investigations devoted to the determination of the regulariza- 
tion parameter a, any more or less constructive algorithms are absent. 
The main reason seems to be the inconsistency of the idea that implies 
a possibility of efficient matching between the solution and the data of 
ill-posed problems (see, e.g., [2, 21, 22]). 

As a matter of fact, one has to be satisfied only by a comparison 
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of solutions to (4.3) obtained in the range of the decrease of a. One 
can assume that because of great labor input of numerical realization 
for small values of the regularization parameter, a large-scale applica- 
tion of Tikhonov methodology to the practice of scientific investigations 
incTirrcd considerable economic damage. As regards attempts to inves- 
tigate the Fredholm integral equation in functional spaces of its correct 
solvability, they were isolated and were not accompanied by construc- 
tive implementation [23]. 

V. M. Pridman, whose papers [24, 25] are considered in section 3.3, 
approached the solution of (4.1) regardless of its applicability to model- 
ing of concrete processes. From the point of view of our consideration, 
the iterative algorithms of Fridman may be of interest, because they al- 
low one to achieve maximal possible efficiency in the framework of the 
chosen object of investigation, which is indirectly confirmed by their 
simplicity and brevity. In other words, it is hardly possible to obtain 
anything more from the traditional interpretation of Eq. (4.1). Despite 
formally existing convergence, by approaching the solution, the deter- 
mined corrections become small against the background of the values 
of the sought function: 

l/jn+l (x) = I/Jn (x) +X[f (x) - (Al/J^) (x)] . 

In the absence of a timely halt of such a procedure, computational 
"noise" from operations with numbers that differ by order of magni- 
tude can radically distort the solution [5, 15]. It becomes obvious that 
the Fredholm integral equation of the first kind, by virtue of its na- 
ture, contains an inherent defect that principally disagrees with pithy 
formulation of the problem of the determination of the function i/j (x) 
from the kernel and the free term of (4.1). 

In section 3.5, wc have given the argument of K. I. Babenko [26] 
for the necessity to take into account the fact of the loss of information 
when evaluating comparative efficiency of computational algorithms. 
This argument seems to be even more important at the stage of the 
formulation of the problem. Since calculations of / (x) from (4.1) ob- 
jectively delete the information on the function i/) (x), its restoration in 
the framework of the traditional approach quite naturally reduces to 
an ill-posed problem. 
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If we hypothctically assume that for the determination of the func- 
tion ip (x) satisfying (4.1) one can find a different equation that contains 
this function not only under the sign of integration but also in an ex- 
plicit form, all the problems will be removed. Such an appearance of 
ijj (x) can be viewed in the context of modeling of computational errors 
including also the integral component (which yields " zero" in the sum) . 

4.4 Methodological concepts of numerical 
simulations 

The predetermined method of conjugate gradients is considered to be 
one of the most efficient methods for the solution of ill-posed systems 
of linear algebraic equations that appear as a result of discretization of 
different problems of mathematical physics [27]. The predeterminer, a 
non-degenerate matrix, allows one to reduce the procedure of numerical 
realization to a sequence of algebraic problems with desired favorable 
properties. On the other hand, however, the number of necessary it- 
erations and the difficulty intermediate calculations increase (section 
3.7). 

One of the key problems of computational mathematics is the de- 
velopment of the conceptual basis for a relationship between a repre- 
sentation of the data and the efficiency of the employed algorithms. In 
this regard, the ideas of K. I. Babenko [26], completely based on a qual- 
itative interpretation of the notion of information can be estimated as 
rather pessimistic. Indeed, almost all computational operations of this 
guide are accompanied by a "colossal" loss of information, whereas rare 
exceptions correspond only to a special representation of initial tables, 
which, as a riilc, is not realized in practice. 

The position of R. W. Hamming [28], who can be characterized as 
a direct follower of the ideas of J. Hadamard in the field of computa- 
tional mathematics, is alternative. In his opinion, methods of numerical 
realization must be adapted to the available information. As regards 
principal difficulties, such as the incorrectness of the formulation, the 
main attention should be concentrated on a modification of mathe- 
matical models. The arguments of P. Bellman and S. Dreyfus for the 
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expediency of the evaluation of the quahty of information on the basis 
of its efficiency indices [29] are also rather attractive. 

O. M. Belotserkovsky and V. V. Shennikov [30] stated a crisis in 
the sphere of numerical simulations resulting from the complexity of 
both the formulations of practical problems and the techniques of their 
numerical realization (section 3.8). As a reason, they have pointed out 
an inapplicability of methods of "domestic" mathematics to situations, 
when owing to the accumulation of round-off errors actually any al- 
gorithm becomes computationally incorrect. As a matter of fact, the 
authors proposed to develop more intensively approaches in the style of 
Tikhonov, without any mentioning of the alternative way, i.e., match- 
ing the formulations of considered problems with Banach's theorem on 
the inverse operator. 

Note that generations of specialists in different fields of mathemat- 
ical physics were brought up under slogans of the type "all real prob- 
lems of the mechanics of continuum medium are ill-posed" that were 
repeatedly reiterated without any explanations by "greats" at different 
conferences. As a result, we have an implementation at a folklore level 
of the thesis supported only by the practice of scientific research. 

N. N. Yanenko, who, in contrast to some colleagues, was well aware 
of the losses of numerical simulations from the breakup of ties of the 
techniques of numerical realization with the basics of functional anal- 
ysis, can be called a flagship of this ideology. However, he considered 
to be of crucial importance the principal difference between classical 
and computational mathematics consisting in the fact that the former 
dealt with abstract symbols without the loss of information, whereas 
the objects of the latter were numerical arrays whose transformation 
was inevitably accompanied by errors of different kinds (see [3, 31]). 

The arguments of the methodologically oriented works of N. N. Ya- 
nenko allow us to suggest that a certain role in the formation of his 
ideas was played by ambitious motivations of being a co-participant of 
the emergence of "new" mathematics that, while partly employing the 
"old" one, was, in general, substantially superior. A grotesque manifes- 
tation of this position is contained in the materials of the monographs 
[21, 32]. Extracts from these monographs are given in section 3.8. 

It seems that we are facing a distortion of the essence of the prob- 
lem, because Banach's theorem on the inverse operator is an entity of 
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a higher level than numerical operations and, at the same time, is most 
important exactly for them. Indeed, the boundedness of the inverse 
operator yields practically a unique possibility to prevent both inade- 
quate dependence of the solution on the data of the problem and the 
accumulation of computational errors. 

4.5 Ideas of the development of a con- 
structive theory 

Thus, let us suppose that the kernel k{x,^) of the Predholm integral 

equation of the first kind (4.1) is symmetric and closed, and the func- 
tion ip (x) satisfying this equation in L2 (0, 1) exists. Correspondingly, 
/ {x) e I2, i.e., the following condition [13] is fulfilled: 

00 J, 

J2 "n^n <00, f (x) i>n (x) dx, (4.4) 

n=l ^ 

where A„, ipn {x) are the characteristic numbers and the eigenfunctions 
of the kernel k{x,^). Note also that the system of elements ^ipn\ is 
complete in R{A) or in the space ['2 ([33], p. 69). 

In this case, the operator A~^ that maps from the space ['2 into 
L2 (0, 1) is bounded (Banach's theorem). Does it mean that the func- 
tion ijj (x) can be determined from (4.2) without accumulation of errors? 

From this point of view, the Inverse World of Banach is rather 
captivating. However, it does not allow for any differentiation of the 
employed spaces with respect to preference. They are determined by 
the content of the problem, i.e., by the operator A. The dominant 
tendencies in the sphere of computational mathematics are purely al- 
ternative. Therefore, both openly and mainly implicitly, introduced is 
the thesis that Banach's theorem on the inverse operator is useless. 

At the first sight, there is a serious reason for this. Indeed, the 
smallness of the perturbation of the data and of the error admitted in 
computational operations is implied in I2. However, a practical possi- 
bility to satisfy this condition is absent. The space I2 is, in a sense, 
illusive because it deals with an infinite set of features of the data of 
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the problem that, for large values of n in (4.2), in essence, cannot be 
identified. 

One can also not that Eq. (4.1) is, in a sense, nonlinear. Indeed, 
let us represent the function, integrated according to (4.1), in the form 
ip — ipi + ip2- Correspondingly, 

1 



and each of these two equations is solvable in the sense of the fulfillment 
of a condition of the type (4.4). 

However, the function / = /i + /2 can be represented as a sum of 
an infinite number of summands. If we assume that the equation 

1 



where + iIj2 = is solvable for an arbitrary subdivision of / into f[ 
and we arrive at a contradiction. Indeed, the solution of Eq. (4.1) is 
unique, and a condition of the type (4.4) is fulfilled only for fi & R {A). 

Thus, the principle of linear superposition does not apply to the free 
term of Eq. (4.1).^ This situation results from the fact that the range 
of the operator A is not closed, which was mentioned in section 2.3. 

In general, the fact that the function f {x), theoretically, belongs 
to I2, in reality, does not yield anything. However, such a conclusion 
cannot serve as a basis for the neglect of the space ^2 in the consideration 
of the problem (4.1). It seems that constructiveness is possible here 
only in the context of the agreement of, generally speaking, alternative 
aspirations: 

- the function / (x), employed in the calculations, belongs to L2; 

- the operator A maps from L2 into 

The motivation is obvious: to preserve the potential of continuous 
inversion of the operator A for practical realization. At the same time, 
the outlined contradiction is clear, and it cannot be overcome exclu- 
sively in the framework of the Fredholm integral equation of the first 



This point partly overlaps the material of section 7.5. 
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kind (4.1). In this situation, it is quite natural to turn, figuratively 
speaking, to the origin of this equation, that is, to the issues related to 
the formulation of the problem. 

Consider a certain process described by the operator A. The direct 
problem consists in the evaluation of the integral according to (4.1) 
under the substitution of the given function ip (x). This procedure has 
a lot of interpretations and is mathematically correct. 

A key element is the formulation of the inverse problem for the 
same operator A, which is related to the restoration of the function 
i/j (,t) from the realization of the above-mentioned integration, that is, 
/ (x). Correspondingly, implied is the determination of the cause from 
its consequence. Whereas the formulation of the direct problem is trans- 
parent, the status of the inverse problem is diametrically opposed. A 
priority of its solution is the actual algorithmic procedure (on the ba- 
sis of an adequate mathematical model) that is not an analog of the 
process occurring in the regime of real time.^ 

In light of the above, is it possible not to turn to the statement 
of Hadamard that all problems having practical interpretation admit 
a mathematically correct formulation? Prom this point of view, since 
the function (x) entering (4.1) objectively exists, the problem of its 
determination has to be only adequately posed. At the same time, 
Hadamard did mot give corresponding recommendations of practical 
character, and, as already mentioned, his methodology turned out to 
be, in essence, completely rejected. 

Let us try, however, to outline a formulation of the problem, inverse 
of the evaluation of the integral (4.1), that is carried out, in general, 
with a certain error: 

AiP^f + Sf, xe[0,l]. (4.5) 

In the direct formulation, taking into account this error has no 
principal importance. Nevertheless, solutions to the Predholm integral 
equations of the first kind (4.1) and (4.5) can be completely different. 

At this point, it is senseless to pose the question about any quantitative 
interpretation of Sf. One can only assume that the error Sf is small 
compared to the values of the functions ip and /. 



^Indeed, the cause as an outcome of the consequence has no physical sense. 
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By general considerations, the presence of 6f' in (4.5) increases the 
potential of the formulation of the inverse problem, and the question of 
a functional representation of the error arises alongside. In this regard, 
one must take into account that the mechanism of its generation is gov- 
erned by the factor of smoothing of ip {x) by the integration procedure; 
therefore, the structure of 6f' must reflect this situation 

In light of the above, let us use an operator model of the error in 
the form 



where I is the identity operator; S is a certain integral operator; 5 — 
IJ,6'; fi and A are parameters. 

Thus, instead of Eq. (4.1), we propose to consider the following 
problem: 



The aim is to reduce this problem to the solution of the Fredholm 
integral equation of the second kind. The parameter n, like A, in the 
inversion of the operator I — XB serves to prevent this equation from 
positioning itself on the spectrum, which is equivalent to the existence 
and uniqueness of its solution. 

Note that we have just added a function representing "zero" to the 
free term of (4.1). At the same time, the transformation of the ill-posed 
problem (4.1) into the formulation (4.6) creates conditions for a radical 
change of the situation. We can demand, generally speaking, that Sf 
adaptivcly compensate for the errors of numerical operations that take 
/ (x) out of the space /g. As a result, a prospects for a realization of the 
bounded operator emerges. For / + Sff/i e R{A), the negative 
factor of the incorrectness of Eq. (4.1) is fully neutrahzed. 

Let us assume that the operator B in (4.6), for which 5f — in the 
spaces C or L2, can be represented in the form 



(4.6) 



= f^f + Sf; 5f = 0,xe [0, 1] . 



(4.7) 



1 




-1 



under certain conditions on the kernel h{x,^). In this case, the problem 
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(4.7) takes the form 

1 1 
^(x)^X I h{x,i)i;{i)di + ii j k{x,i)i;{i)di-iif{x), (4.8) 

-1 

1 

i;{x)^\jh {x, V' (0 X e [0, 1] . (4.9) 

-1 

Thus, the condition that 5f be equal to zero, which equivalent to Eq. 
(4.9), is supposed to be satisfied with the help of (x) on x e [—1, 0), 
i.e., a new unknown function. 

There exists a well-known opinion that prospects of obtaining new 
substantial results by simple transformation of mathematical relations 
are not great. Indeed, by applying to Eqs. (4.8), (4.9) a subtraction 
operation we again obtain the initial problem which is ill-posed. How- 
ever, first, we are not going to do this, and, second, behind the integral 
equation with the sought function in an explicit form, we intuitively 
feel a constructive potential. 

Prom this point of view, a "refusal" of the well-known example, 
given in a number of references, that illustrates the smooting of infor- 
mation about the function ip by means of integration of (4.1) seems 
to be very significant. Indeed, assuming that the function ip = tp^ {x) 
satisfying the system of equations (4.8), (4.9) is known, we give it a per- 
turbation of the type esin(n7ra;). A substitution into (4.7) shows that 
this perturbation influences the free term / (,t) both via a reduction 
coefficient (smoothing) and without it, at the expense of an integral 
component and of explicit presence of ip {x), respectively. 

What is said does not apply to i/^ix), x e [—1,0). However, the 
determination of this function is beyond the scope of the considered 
problem. We want to emphasize that the latter arguments bear exclu- 
sively heuristic character. 

From the position of a practical realization of the above, an inter- 
relation between the spaces L2, h and I2 seems to be rather significant. 
As is well-known, it is tightest in the pair of the spaces L2 and l2- 
The Riesz-Fischer theorem [34] establishes a one-to-one, continuos and 
linear relationship between functions from L2 and numerical sequences 
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{cji} with a convergent sum of the squares. In other words, there always 
exists a L2-function for which 

oo 
n=l 

is a Fourier series in terms of a system of orthonormal elements {ipn {x)}. 

However, there is also a rather interesting relationship between the 
spaces I2 and Zg, and, correspondingly, L2. Indeed, equation (4.2) repre- 
sents a Fourier series in terms of the orthonormal elements ipn (x) , whose 
convergence condition is given by (4.4). If we assume that A„ = r"**, 
where < r < 1, the space I2 turns into I2 under the condition r ^ 1. 

At the same time, the kernel k{x,^) in (4.1) possesses objectively 
inherent characteristic numbers and, consequently, cannot be used for 
such transformation. However, there appeared the kernel h (x, which 
is independent of the data of the problem: hence a prospect of achieving 
what we set out to do. A considerable part of our consideration below 
will be focused on this issue. 

In conclusion of this section, we want to point out the inconsistency 
of the wide-spread opinion that the formulation of problems of numer- 
ical simulation should be left to specialists in applied sciences, whereas 
pure mathematicians should be concerned exclusively with rigorous an- 
alytical investigations, the development of computational methods and 
participation in their realization. 

It seems that specialists in applied sciences should be concerned with 
the formulation of direct and, generally, well-posed problems. The fac- 
tor of incorrectness is directly related to the procedure of the numerical 
realization. Therefore, the main concern of pure mathematicians should 
be a reduction of formulations of problems describing the considered 
processes and phenomena to the conditions of efficient implementation 
of Banach's theorem on the inverse operator. 
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Chapter 5 



A method of the reduction of 
problems, traditionally 
associated with Fredholm 
integral equations of the first 
kind, to Fredholm integral 
equations of the second kind 

5.1 The formulation of the problem 



In light of the arguments of section 2.4 and 4.5, we proceed with the 
consideration of the Fredholm integral equation of the first kind 



{A^|J) {x) = jk (x, (0 di^fix), xe [0, 1] (5.1) 



under the assumption that its solution exists and is unique, and the 
kernel k {x, ^) and the free term / (x) belong to the space L2. In other 
words, using the terminology of [1], they are L2 - functions: 



1 1 




(x) dx < 00. 
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However, in reality, the determination of the function ip (x) from 
given A and / will carried out not by the use of the solution of the 
Predholm integral equation of the first kind (5.1), but on the basis of 
the following arguments. There is an operator A describing a certain 
phenomenon. This description is expressed in terms of the integration 
of the function (x) G L2 (0, 1) by (5.1). 

The evaluation of / (x) is carried out with an error that we de- 
note as {Sf) {x) / n, where is a constant. In most cases this error, in 
virtue of its smallness, is nonessential or can be reduced to a required 
level. Nevertheless, the computational procedure can be interpreted as 
follows: 

{A^){x)^f{x) + {6f){x)/i,, xe[0,l]. (5.2) 

The situation changes cardinally if, on the contrary, we pose a prob- 
lem of the restoration of the function ip (x) from the information con- 
tained in (5.1), i.e., A and /. Indeed, such a problem is, in general, 
ill-posed, which, in fact, means that Eq. (5.1) is insolvable. 

Prom this point of view, Eq. (5.2) is different because of the presence 
of a potential of the reduction of the problem to a well-posed one. A 
necessary condition of this reduction consists in such a representation 
of the error 5f that, irrespective of the data (5.1) and of the function 

ij{x), 

f{x) + {6f){x)/i2eR{A), (5.3) 

where R (A) is the range of the operator A. In other words, the operator 
S [see (4.6)] must endow the algorithm with adaptive properties. 

Thus, the following problem is posed: Prom given A and /, deter- 
mine constructively the function -0 (x) that, upon substitution in (5.1), 
would satisfy this equation. Here, constructiveness implies a possibil- 
ity to use a stable procedure of the numerical realization as a result of 
the reduction of the problem to the solution of he Predholm integral 
equation of the second kind.^ 

The basis of further transformations will be formed by Eq. (5.2), 
where the central point is the establishment of adequate mutual de- 
pendence of ip and Sf. Equation (5.1) is considered exclusively in the 



^It is supposed that the kernel of this equation does not possess any singularities 
incurred by the method of the realization of corresponding transformations. 
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context of the direct problem of the evaluation of the integral and as a 
source of initial information. 

5.2 The model of the representation of 
the error 

Following the considerations of section 4.5, we present the error of the 
evaluation of / from (5.1) as a difference between the sought function 
■0 and the integral component 

{Sf){x)^i;{x)-X{B^){x), xe[0,l], (5.4) 

where A is a constant; the operator is given by 

1 

B. = J h{x,0.dC] (5.5) 

-1 

■0 (x) = ip{x), X G [—1,0); the kernel h (x,^) will be discussed later. 

However, we intend to construct a stable algorithm of evaluation of 
the function ip (x) satisfying (5.1); hence small variations of the data 
should not substantially influence the solution. In this regard, consider 
a possibility of the fulfillment of the condition 

{Sf){x)^0, xe[0,l], (5.6) 

which means an assumption that the problem posed in section 5.1 can 
be constructively solved (merely) by means of addition to the free term 
of Eq. (5.1) of the "zero" from (5.4) that has the following form:^ 

1 

o = V'W-A / /i (x, V' (0 ^e- 

-1 

^Here, the error 6f or the function dependent on this error are interpreted as a 
component of the free term of the Fredholm integral equation of the second kind, 
employed for the determination of tjj. 
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This equation can be rewritten as 





J hix,C)<f{OdC^g{x), xe[0,l], (5.7) 



where 

g{x) = -^{x)- j h{x,i)^{Odi. (5.8) 



Making the change of variables 

C-27ra;, ^ = 27r(l + e), (5.9) 
we reduce it to the canonical form 

27r 

jh{U)v{e)de^g{C), Ce[0,27r]. (5.10) 



As is obvious, the satisfaction of (5.6) is equivalent to the solvability 
of this equation. Let the kernel h (^, 9) belong to the space L2 and be 
closed. In this case, Eq. (5.10) is a Fredholm integral equation of the 
first kind, whose the solution, if it exists, is unique [1]. By satisfying the 
above conditions, we represent (5.10) in the form of a Poisson integral 
([2], pp. 202-205). Accordingly, the kernel is given by 

h{U)- ^ ri 9 2V 0<r<l; (5.11) 

27r [1 — 2r cos [Q — 6) + 

its characteristic numbers and orthonormal on x e [0, 27r] eigenfunc- 
tions ([1], pp. 187-188) are 

'^O — 1) ^2n-l — ^2n — f ^ ■> U = 1,2, . . . ] 
<?0 (C) = <^2n-l (C) = ^ COS «) , 



(^2n(C) = ^sin(nC), n=l,2,...; (5.12) 
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and, in (5.7), 



^2 



1 — 2r cos [27r [x — ^)\+r^ 
If, in Eq. (5.10), the function 

2 oo 

^ (C) = 7;<^o + J2 COS (^C) + "n sin «) , (5.14) 

^ n=l 

where ckq, and a'^ are the coefficients of its expansion into the Fourier 
series, is absolutely integrable, i.e., 

/ \v(0\dC<oo, 



the function 

-j^ oo 

9 (C) = T^«o + X] ^" K COS «) + a'^ sin «)] (5.15) 

^ n=l 

is the real part of an analytical inside a unity circle function and is 
harmonic ([3], pp. 160-161; [4]): 

where X = r cos {Q, Y = r sin {() arc Cartesian coordinates.^ 

Since the above-mentioned property is independent of a linear change 
of variables, it follows from (5.8) with (5.9) and (5.15) that, under the 
condition (5.6), the function i(j (x) satisfying (5.1) can only be harmonic. 
This means that it belongs to a much narrower class of functions than 
it is supposed in the formulation of the problem in section 5.1. 

Nevertheless, one can conclude that the "zero" error of integration 
by (5.1) of the harmonic function ip (x) is actually representable in the 
form (5.4) with the kernel h{x,$,) from (5.13). This is an important 
point of our consideration. 



^Here, the parameter r is interpreted as a radial coordinate and C is, respectively, 
a polar angle. 
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The components (5.12) satisfy the homogeneous equation 



2n 



ip{C)^xJh{c,e)ip{e)de, c^[oM 



that, by the change of variables 

C = -27rx, 9 = -27r^; C = Sttx, 9 = 27r^ 
is transformed to the following form: 



ip{x)=xj h{x,0^{Od^, a;e[-l,0); 



-1 

1 



^{x) = xJh{x,C)^{C)dC, xe[0,l], (5.16) 



which allows us, taking also account of (5.12) and (5.9), to determine 
the characteristic numbers and the orthonormal on a; £ [—1,0); [0,1] 
eigenfunctions of the kernel (5.13): 

Xo — 1, X2n-l — X2n — f n = 1, 2, . . . ; 

ipo (x) = 1, <^2n-i (x) = \/2 COS {2nnx) , 

'^2n (x) = V2 sin {2nnx) , n=l,2,.... (5-17) 

The solution of the problem (5.1) is unique. Accordingly, by com- 
paring the homogeneous Predholm integral equation of the second kind 
with respect to ip {x) that corresponds to (5.8) (i.e., for g = 0) with 
(5.16), we arrive at the condition 

A^r-", n = 0,l,.... (5.18) 

As the kernel in (5.16) is symmetric, continuos, and all A2„ > 0, by 
Mercer's theorem [1], 

^ ^ <fO jx) (PQ (0 ^ g ^2n-l {X) (P2n-l {Q + <f2n jx) <f2n (Q 
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= 1 + 2 [cos {2n'Kx) cos (2n7rO + sin {2n'Kx) sin {2n7r^)] , (5.19) 

n=l 

where the series can be absolutely and uniformly convergent. 

In what follows, we will need the resolvent of the operator B. Prom 
the bilinear expansion (5.19), by same Mercer's theorem, it follows 
that the characteristic numbers and the orthonormal on x e 1] 
eigenfunctions of its kernel have the form 

Ao = 2' '^2n-l = ^2n ~ 2^ ^ = 1,2, . . . ] 

V'o ^ ^2n-i (x) = COS {27mx) , 

i^2n (x) = sin {27Tnx) , n = 1, 2, . . . ; 
hence a necessity to impose one more condition: 

At^V", n = 0,l,.... (5.20) 

One should take into account that the use of Mercer's theorem is 
different from the former representation of the kernel h {x, ^) be the 
series (5.19). Here, on the contrary, there exists an expansion of the 
kernel h{x,^) into a uniformly convergent bilinear series in terms of an 
orthonormal on — 1 < x < 1 system of elements. Accordingly, these 
elements, under a correction with respect to a normalization factor and 
the value l/2r", are the eigenfunctions and the characteristic numbers 
of the operator B. 

We also note that the functions ■02n-i (x), ■02n (x) are orthogonal 
not only on a; G [—1, 1], but on a; G [—1,0); [0, 1] as well. This point 
will play a rather important role in the context of the simplification of 
the procedure of the numerical realization. 

The resolvent of the kernel (5.5) is represented by the series [1] 

i/(.,e,A) = M!iM) 

AO — A 

_^ g tp2n-l jx) -02n-l {Q + '02n jx) -02n (Q ^ 1 
n=l "^^n ~ ^ 1 — 2A 
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'^^ ~\ — ^r)\ n ^"^^^ (2n7rx) cos (2n7r^) + sin (2n7rx) sin (2n7r^)] (5.21) 

^ 1 X ^ /A / 



n=l 



that, under the condition (5.20), is also absolutely and uniformly con- 
vergent. 

Prom (5.8) and (5.15), taking into account (5.9), we get: 



^ ^ 2 (1 - A) ^ ^ 1 - Ar^ {2n7rx) + a„ sin {2nnx) 



Thus, under the condition (5.6), Eq. (5.1) can be satisfied only in the 
case when 

OqA 



(l-A) 

oo ^^n 



+ ^1 

n=l ^ 



k (x, ^) [an cos (2n7rx) + sin (2n7rx)] d$,. (5.22) 







In what follows, we assume that the function i(j (r, a;) is harmonic 
and the free term of Eq. (5.1) has the form (5.22). As already men- 
tioned, this fact strongly narrows the sphere of practical applications. 
As will be shown below (section 5.6), a solution, obtained for this case, 
by means of the passage to the limit r — > 1 turns into an L2 - function 
■0 (x) that satisfies Eq. (5.1).^ 

5.3 A transformed formulation of the prob- 
lem 

Let us extend Eq. (5.4), under the condition (5.6), in the following 
way: 

1 

^{x) = X J h{x,0^{Odi + >^J h{x,0^lj{Od^ 
-1 

+k{x), a;e[-l,0), (5.23) 



^Simultaneously, Eq. (5.22) takes the form (5.1). 
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where k{x) G Li (—1,0), as a result of (5.14), is a certain undefined 
function. 

We represent the equation that unifies (5.7) and (5.23) in the fol- 
lowing form: 

\fr^ \ ti{x), xe [-1,0), ^ ^ 




I.e., 





b(1']{x) = J h(x,0'p{Od^ + jh{x,0^{Od^- 

^ ^ -1 

Let us introduce also a related and close with respect to its structure 
equation 

^J + l 0, xe[-i,o), ^^-^^^ 

where ip' {x) and % {x) are two more undefined functions (like -i/^, they 
are harmonic). The expediency of this step will be clear from what 

follows. 

It is not difficult to represent the procedure of the construction 
of Eqs. (5.24) and (5.25) from the practical point of view. There is 
a harmonic function ip [x) that is integrated according to (5.1). As is 
shown above, there exist the kernel h{x,^) and an absolutely integrable 
function (p{x) for which Eq. (5.24) is satisfied on x G [0,1]. One can 
assume that the function (/? (x) is specified in a certain way. Now both 
Ip (x) and (x) are given functions. Equation (5.24) is satisfied by 
means of the function k,{x) on a; G [—1, 0) and on the whole. 

The function ip (x) is again given. The function (p' (x) is determined 
from Eq. (5.25) on x G [-1,0): 



^' {x)^\j h{x,i)^'{Odi + g'{x), 



1 

g'{x)^xj h{x,i)i^{i)di. 





(5.26) 
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This is a Fredholm integral equation of the second kind with respect 
to ip' {x). According to the foundations of the general theory [1], under 
the condition (5.18), the solution of (5.26) exists and is unique. The 
functions ip (x) and (/?' (x) are given, and Eq. (5.25) is satisfied by means 
X {x) on a; G [0,1] and on the whole. 

In terms of the notation 



^' (x) 



^ (x) , X e [0, 1] ; 
ip{x), a;e[-l,0). 



Eqs. (5.24), (5.25) are Predholm integral equations of the second kind 
with respect to ^ and with the free terms 



P{x)^\ 
P' (x) = 



0, xe [0, 1] ; 
k{x) , X e [—1, 0) , 



X (x) , X e [0, 1] ; 

0, xe[-i,0), 

respectively. 

Under the condition (5.20), the solutions of these equations are given 

by 



^l;{x)^xl H{x,^,X)K{Od^, xe[0,l]; (5.28) 
-1 



^{x) = K{x) + xjH{x,^,X)K{Odi: a;e[-l,0) (5.29) 



-1 

and 



J. 

,l;{x)^x{x) + xjH{x,C,X)x{OdC, ^e[0,l]; (5.30) 





1 



ip'ix)^xlH{x,C,X)x{OdC, xe[-l,0), (5.31) 
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where H {x, ^, A) is the resolvent of the operator B that has the form 
(5.21). 

By subtracting (5.25) from Eq. (5.24), we get 



x(x)^xJh(x,0[^(0-'P'md^, ^e[0,l]; (5.32) 





ip{x)-ip'{x) = X J h{x,0[^{0-v^'iO]d^ + ^^ix) ^e[-i,o). 

-1 

(5.33) 

Prom these relations, it follows that the function x can be construc- 
tively expressed via k, i.e., by means of the solution of the Fredholm 
integral equation of the second kind. Indeed, under the condition (5.18), 
ip — ip' is determined via the resolvent of the kernel h{x,^) in (5.33). 
However, the inverse procedure, i.e., a representation of the function k 
via X, would be related to the solution of the Fredholm integral equation 
of the first kind. 

Let us add to Eqs. (5.24), (5.25) the "zero" from (5.1), i.e., fiAip — 
nf with the free term of the form (5.22). As a result, we obtain, 
respectively. 



^i^) \_^J.(i^\(,^,^^^iA^|;){x)-|lf{x), x G [0, 1] ; 

(5.34) 

^'{x) j - [ ^' ) ^''> + \ 0, xe[-l,0). 

(5.35) 

Thus, instead of the ill-poscd problem (5.1), in what follows wc will 
consider the two systems of integral equations (5.24), (5.34) and (5.25), 
(5.35).^ 



^Note that (5.34), (5.35) do not constitute Fredholm integral equations of the 
second kind with respect to the functions (5.27). 
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5.4 A constructive algorithm of practical 
realization 

A further orientation of transformations is, in a sense, opposed to the 
previous one. Indeed, above, in a fact, we have done our best [beginning 
with the model of the error (5.4)] to ensure that the sought function 
ip (x), as well as (/? (x) and (/?' (x), appear in specially constructed equa- 
tions not only under the sign of integration but also in an explicit form. 
As a consequence, we have obtained (5.30), a representation of the so- 
lution ijj {x) with the function x i^) also in an explicit form. 

It would be highly desirable to derive a different representation of 
i/j (,t) that would apparently contain the data of the problem (5.1) and 
where the function x (x) would appear only under sign of integration. 
Upon elimination of the function tjj (x) both from this representation 
and from (5.30), we could obtain a Predholm integral equation of the 
second kind with respect to x {^)- 

Another way of achieving the same goal consists in the determina- 
tion of the integrand (5.32) via x{^)- Since the function ^p' (x) is, in 
this sense, known [see (5.31)], it is necessary to establish a relationship 
between (p, x and the data of the problem. 

The realization of each of the two outlined versions can be repre- 
sented in the context of the reduction of (5.35) to the form (5.34). The 
grounds for this reduction lie in the fact that the function ip (x) enters 
both the equations and that their structure is analogous. These are 
heuristic arguments. 

In order to eliminate the function x (x) from (5.35), we use the 
equation 

fi (Ai^o) (x) + X (x) , xe [0, 1] ; 
0, xe[-l,0). 

(5.36) 




By subtracting this equation, we get 



(x) - ipfi (x) \ _ .p, ^ - ^0 \ /X 

V.' (x) - (x) j """"[v'- V'o ) ^ ^ 
+ 1 0, a;e[-l,0), ^^-^^^ 
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or 



tpi {x) 
ip[ (x) 



XB 



^1 



(x) + 



fi{A,j;,){x)-txf{x), xe[0,i]; 
0, xe[-l,0), 

(5.38) 



where 

ipi (x) ^ijj{x) - ipo (x) ; (p[ (x) = (x) - (p'o (x) . (5.39) 
If introduce the notation 



equation (5.38) takes the form 
^i{x) = X{B^i){x) + \ 



(x), xe [0, 1] ; 

ip[{x), xe[-i,o), 



II (AVi) (x) - l^f (x) , X e [0, 1] ; 

0, xe[-i,o). 



This is a Predholm integral equation of the second kind with respect 
to The inversion of the operator I — XB under the condition (5.20), 
taking into account (5.1), yields: 

1 

i^,{x)^xj K{x,^)Mi)di + h{x), xe[0,l]; (5.40) 



where 



^\ {x) ^f,xj ^1 (0 / H {x, C, A) k (C, dC 

L 

-H{x,tX)f{0]d^, xe[-l,0), 
1 

K {x, = k{x,0 + xjH {x, C, A) k (C, dC; 



J. 

f{x) + xj H{x,i,x)f{i)di 



(5.41) 
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Thus, the function ipi (x) is determined from the Fredholm integral 
equation of the second kind (5.40) and depends only on the data (5.1) 
and on the chosen kernel h{x,^).^ Here, we assume that 

li^l^n, n=l,2, 

where /i„ are the characteristic numbers of the homogeneous equation 
obtained from (5.40) in the case /i = 0. The values of as well as the 
solution of (5.40), should be found by means of approximate methods 
[5]. After that, the function Lp[ (x) is evaluated from the formula (5.41). 

However, Eq. (5.37) can be regarded as Eq. (5.34). Indeed, the 
elimination of x i^) from Eq. (5.35) is, figuratively, equivalent to a flow 
of this function to k,{x) with the appearance of Eq. (5.34). Conse- 
quently, what is needed is an identification of the functions (p (x) and 
K (x) on the basis of (5.34) in the structure of Eq. (5.37). 

To this end, we use Eq. (5.37) on x e [— 1, 0), 



(x) -^',{x)^\jh (x, y (0 - v'o m 
-1 

1 

+A j /i (x, [V' (0 - V'o (01 (5.42) 



paying attention to the method of its derivation. It consists in the 
elimination from Eq. (5.35) of the part of the solution that depends 
on the component of the free term x {x)- However, in this procedure 
the functions satisfying this equation both on x e [— 1,0) and on x e 
[0, 1] have changed. In other words, both the functions </?' and have 
undergone change. 

At the same time, the structure of Eqs. (5.34), (5.35) implies a 
transformation of one of these equations into the other by means of a 
change of the contained functions only on x e [— 1, 0), that is, of </?' and 
</7.^ Therefore, we will correct </?' {x) in Eq. (5.42) in order to eliminate 

^As a matter of fact, Vi represents the part of the solution (5.35) on a; G [0, 1] 
that is stipulated by the component of the free term — /x/. 

^What was in position of V' in (5.35) must remain unchanged. 
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the term with the function ipQ (x) . Accordingly, we must include in 
K {x) the terms of Eq. (5.42) that contain the function ipQ (x). 
As a result, there appear the relations 



K{x) = ip',{x)-Xj h{x,0^o{Od^, a;e[-l,0); (5.43) 
-1 

(p{x) ^ (f' {x) - ipo{x) , xe[-l,0). (5.44) 

Here, (/9o (x) is the solution of the Predholm integral equation of the 
second kind 



{x)^\jh (x, i) ipo iO + fo{x), xe [-1, 0) , (5.45) 

where 



-1 



/o (x) --^Jh (x, V'o (0 d^, 



under the condition (5.18). 

Subtracting (5.45) from Eq. (5.42), we get 



^' (x) - (x) -^o{x)^\jh (x, y (0 - ¥^0 (0 - (01 di 

-1 



+\j h{x,Oi^{i)di, xe[-l,0). (5.46) 
-1 

Equation (5.35) on x e [—1, 0) has the form 

1 

^'{x)^\jh (x, {i)di + \jh (x, V' (0 di. (5.47) 



Its comparison with (5.46) yields: 

(^0 {x) = -^'q {x) , (5.48) 
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that is, we have, in fact, returned from (5.42) to Eq. (5.35) on x G 
[—1,0) in such a way that allows us to establish this relation. 

It should be noted that relations (5.43), (5.44) transform (5.46) into 
(5.23). Now we will show that relations (5.43), (5.44) and (5.48) indeed 
reduce Eq. (5.35) to the form (5.34). To this end, we turn to Eq. (5.35) 
on x e [0, 1]: 

1 
-1 

+fi{Aij;){x)-fif{x)+x{x). (5.49) 
Using (5.44) and (5.32), we get: 



A I h{x,0 v'iOd^^xj h(x,0v(0d^ + X I h(x,0^o(Od^, 
-1 -1 -1 

where 



A I h{x,0'Po(Od^--X I h(x,0 [^(^)-^'(^)]d^^-x(x), 
-1 -1 

which, by means of substitution of the above expressions (5.49), is 
transformed into Eq. (5.34) on a; G [0, 1]. 

The substitution of the function (/?' (x) from (5.44) into (5.47), with 
the use of (5.48) and (5.43), leads to Eq. (5.34) on a; G [-1,0).^ Thus, 
by means of the established relations, Eq. (5.35), both on x e [—1,0) 
and on a; G [0, 1], is transformed into Eq. (5.34). 

By (5.44), (5.48) and (5.39), 

If (x) - If' (x) = -ipo (x) = ip'o (x) = If' (x) - ip[ {x) , (5.50) 
and, as a result, expression (5.32) takes the form 



X{x) = \jh (x, i) y (e) - (^i (0] di. (5.51) 



'What was in position of V' in (5.23) must remain unchanged. 
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The derivation of relations (5.43), (5.44) as well as (5.48) and, fi- 
nally, (5.50) is the main link in the construction of the algorithm. 

The substitution of expression (5.51) into (5.31) leads to a Predholm 
integral equation of the second kind: 

1 

x{x)^Allix,OxiOdC + Qix), ^e[0,l], (5.52) 



where A = A^; 



l(x,0 = J h{x,C)H{C,^,X)dC; (5.53) 
-1 



q(x)^-X I h(x,0v[(0d^- (5-54) 

-1 

Expression (5.53), after the substitution of (5.19) and (5.21), takes 
the form 

Z (x, — r + 2 ; — [cos (2n7rx) cos (2n7rO 

V >sy 1 _ 2A ^1 1 - 2Ar" ^ ^ J y 

+ sm {2mrx) sm (2n7r^)J = — — — 

1 — 2A 

°° /l - 2Ar"\~^ 

+ I ) ['P2n-l (x) (p2n-l (0 + V2n (x) (p2n (01 > xe[0,l], 



n=l 



(5.55) 

where <^2n-i (x) and (p2n (x) are the eigenfunctions of (5.17), orthonor- 
mal on X e [0, 1]. This fact allows us to determine the resolvent of the 
kernel l{x,^). Indeed, its characteristic numbers are 

1 — 2Ar" 

Ao = l-2A; A2n-1 = ^2n = ' n = 1, 2, . . . , 

and, because of the property < r < 1 for a bounded A, which is 
assumed, only a limited number of these can take on negative values. 
By Mercer's theorem [1], expression (5.55) is a bilinear expansion of the 
symmetric continuos kernel I (x,^), < x,^ < 1. Under the condition 

At^Ao; A7^A2„, n=l,2,.... 
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which is equivalent to (5.20), its resolvent takes the form 



j.2n 



1 - 2A - A2 

[cos {2m:x) cos (2n7r^) + sin {2m:x) sin (2n7r^)] . 



^1 1 - 2Ar" - 

(5.56) 

As a result, the solution of (5.52) can be represented as follows: 

1 

X (x) = g (x) + A I L {x, e. A) q (0 d^. (5.57) 



Obviously, for the convergence of the scries (5.56), in addition to 
(5.18) and (5.20), it is necessary that the following condition be fulfilled: 

A (-1 ± V2) r-", n = 0, 1, . . . . (5.58) 

The substitution of expression (5.57) into (5.30), by use of (5.21), 
allows us to evaluate the function ip {x) , which is the solution of the 
considered problem. 

The procedure of the numerical realization includes the following 
stages: 

- concretization of the parameter < r < 1; 

- determination of the parameter A from the conditions (5.18), (5.20) 
and (5.58), taking also account of (5.4), that is, 

A^O, A^r-", A^^r-", A^(-l±v^)r-^ n = l,2,...; 

(5.59) 

- determination of the parameter /i in (5.40), so that the equation 

1 

V'(x)=/i|x(x,OV'(0^e, xe[0,l] (5.60) 


possess only the trivial solution; 

- determination of the function ipi from Eq. (5.40); 

- evaluation of the function ip[ by formula (5.41); 
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- evaluation of the function q by formula (5.54); 

- evaluation of the function x by formula (5.57); 

- evaluation of the sought function ip by formula (5.30). 

Note that the reahzation of the algorithm is related with the use 
of quadrature and cubature formulas on a two-dimensional domain [6] . 
Simultaneously, one can apply the technique of the improvement of 
convergence of trigonometric series ([7], pp. 187-193) and the methods 
of integration of oscillating functions ([8], pp. 112-115). 

5.5 The reliability of the obtained results 

Thus, the function ip (x), i.e., the solution of the problem (5.1) in its 
restricted formulation (see section 5.2), is determined by formula (5.30). 
At the same time, expressions (5.28)-(5.31) that represent the solution 
of Eqs. (5.24), (5.25) satisfy these equations identically, irrespective of 
the form of k (x) and x (^)- 

Therefore, one cannot argue on the basis of simple subtraction of 
Eq. (5.25) from (5.35) that the solution of the former equation also 
satisfies Eq. (5.1). In the general case, solutions of these equations can 
be completely different. 

Accordingly, one has to show that the function ip (x), determined by 
expression (5.30), that together with ip' (x) satisfies Eq. (5.25) is also 
the solution of Eq. (5.35) on a; e [0, 1]. To this end, it is reasonable 
to introduce new notation for the functions ip (x), ip (x), ip' (x) entering 
Eqs. (5.24), (5.25) and (5.34), (5.35), namely, ip{x), (fi{x), (f'{x) and 
Ip (x), (f (x), (p' {x\ respectively. 

By use of the above-mentioned pairs of equations, respectively, the 
following relations have been obtained in section 5.4: 



X{x)^\\h (x, i) {(p (0 - (p' di (5.61) 
-1 

and 
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[see (5.32) and (5.39), (5.44), (5.48)] or 

(f (x) - (p' {x) = -(^0 {x) = ^0 (3^) = - 
[see (5.50)]. Finally, in a short form, 

(p {x) - (p' {x) = (p' {x) - (^'1 {x) . (5.62) 
However, only the relation 

^ {x) - {x) = {x) - 'P'l {x) , (5-63) 

where 

1 

^'{x) = xJh {x, e. A) X (0 d^, 


has been used [see (5.31)]. Upon substitution into (5.61), this leads to 
the relation 



X{x) = xj h{x,0 [^' (0 - ^'1 (0] 
-1 

[see (5.51)], which has, as a result, Eq. (5.52). 

In other words, we have substituted the function (f' (x) into the 
right-hand side of (5.62), in place of (p' (x). In general, the derivation 
of (5.63) has been as follows: 

(p (x) - (f' (x) = ip{x)- (f' (x) = (f' (x) - (p\ {x) = (p' {x) - (p\ {x) , 

i.e., the two premises 

^ {x) - ^' {x) = (p{x)- (p' {x) , <p)' {x) = (p' {x) 

and relation (5.62) have been used. 

Indeed, when these identities are satisfied, relation (5.63) turns into 
(5.62). Thus, one can conclude that the above-mentioned premises, i.e., 
the identities 

(p{x)=(p [x] ; (p' {x) = (p' [x) , (5.64) 

constitute sufficient conditions for the reduction of the problem (5.24), 
(5.34) and (5.25), (5.35) to the solution of Eq. (5.52). 



5.5. THE RELIABILITY OF THE OBTAINED RESULTS 



115 



At the same time, they are also necessary. Indeed, by (5.62), (5.63), 

(fi{x)-(fi{x)^2[(fi'{x)-(p'{x)]. (5.65) 

Analogously, that is, by subtraction of Eqs. (5.34), (5.35) from (5.24), 
(5.25), respectively, we get 



^ix)-^{x)^X I h(x,0 [<^(0-<^(0]rfe 



+A / h {x, ^ (0 - ^ (0 dC, 



^' {x) -ip'{x) = \ h (x, (0 - ^' m 



-1 



+XJ h{x,o [i'iO-i^iOl dC- 



As a result of the subtraction, with the use of (5.65), there arise 
homogeneous equations: 



ip{x)-ip{x) = X h {x, [ifi (0 - ^ (0] dC; 



^' (x) -^'{x) = \jh (x, W (0 - ^' m di, X e [-1, 0) , 

-1 

whose solution under the condition (5.18) is trivial. Thus, relations 
(5.62), (5.63) automatically result in the identities (5.64). In other 
words, the existence of the above-mentioned relations imply that the 
functions </?(x), (x) in Eqs. (5.24), (5.25) and (5.34), (5.35), respec- 
tively, are the same. 

The procedure of subtraction in each of the pairs of the equations 
yields: 



h{x,i) i^(0-^(0 di^Q, [-1,0), 
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This is a homogeneous Fredholm integral equation of the first kind 
with the kernel (5.11). As it is closed, we can conclude that 



Consequently, in order that the functions ip (x), if' (x), determined 
by formulas (5.30), (5.31), satisfy both Eq. (5.25) and Eq. (5.35) as 
well as their difference, Eq. (5.1), the function x (x) must represent the 
solution of the Fredholm integral equation of the second kind (5.52). 
This is a very important point of the whole consideration. 

Note that, instead of (5.64), one could employ a single identity 



However, in this Fredholm integral equation of the second kind, 

obtained by the substitution of expression (5.63) into (5.61), would be 
more cumbersome. 

5.6 An arbitrary function from L2 as the 
solution 

Beginning from section 5.2 and up to the present point, we have as- 
sumed that the function %b (x) satisfying Eq. (5.1) can be only har- 
monic. Accordingly, its free term / (x) is determined by expression 
(5.22). Here, we present a generalization of the algorithm of section 
5.4. To this end, we will employ an approach which is analogous to 
Abel-Poisson's method of the summation of Fourier series [2, 3]: 

- execution of the transformation in an analytical form with a har- 
monic function ip (r, x) that is represented by a well convergent series 
for < r < 1 in (5.19); 



C = 27r(l + x); 9^27T^, 



we get 




ijj{x) = ip (x) = -00 {x) + i>\ {x) = -00 (x) + i>\ {x) =ip{x) . 



(f (x) - (f' (x) =^{x)~ ^' {x) . 
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- a passage to the limit r ^ 1 in the expression for ijj (r, x) via 
the data of the problem that is represented by a series whose terms 
expUcitly depend on the parameter r. 

In this way, we will obtained the solution of the problem posed in 
section 5.1: namely, the restoration of the L2 - function (x) from the 
results of integration according to formula (5.1) or from a given related 
expression for / [x) . 

In Eq. (5.40), we use the following representations: 

-j^ 00 

■01 (x) — -So + X/ '^^ {2mrx) + sin {2mrx) , (5.66) 

2 n=l 

where Sq, s„ and s'^ are undefined coefficients; 

k {x, e) = -ko (e) + E (0 cos {2n7ix) + k', (0 sin {2n7ix) , (5.67) 

^ n=l 

00 

/ (^) — 0*^0 + X] {2n'Kx) + sin {2m:x) , (5.68) 

2 n=l 

where the Fourier coefficients are given by 

1 1 
/co (0 = 2 j k {x, ^) dx; /c„ (^) = 2 J k {x, ^) cos {2n7rx) dx; 


1 

k'^{$) ^2 j k{x,$)sn\{2m:x)dx, n = l,2,... (5.69) 



(note that explicit evaluation of these functions is unnecessary); 

1 1 
Co = 2 y / (0 dC; c„ = 2 y / (0 cos {2m:i) d^; 


1 

d, = 2j f{i)siTi{2mii)di, n=l,2,.... (5.70) 


Accordingly, 



nSCHAPTER 5. A METHOD OF THE REDUCTION OF PROBLEMS, TRADITION. 



°^ 1 — Ar" 

+ E 1 _ 2Xrn ''''^ (2'''''') + ^" "^"^ (2''''^)] ' (^-^^^ 



n=l 



- E I _ 2Xrn (SnTTx) + < sin (2n7ra;)] . (5.72) 

On substitution of expressions (5.66)-(5.68) into Eq. (5.40) and 
reduction of the factors multiplying cos {2mTx), sin {2n7cx), the evalua- 
tion of the coefficients So, s„, s'^ reduce to the solution of the following 
linear algebraic equations: 

oo 

[2 (1 - 2A) - /i (1 - A) Poo] So = 2n{l- X)J2 POmSm+P'om^'m-'^l^ (1 " ^) ^0] 

m=l 

oo 

2 [1 - 2Ar" - /X (1 - Ar") p J s„ = (1 - Ar") p„oSo+2/x (1 - Ar") ^ p^mSm 

m=l, rriT^n 

oo 

+2/. (1 - Ar") ^ - 2// (1 - Ar") c„; 

m=l 

oo 

2 [1 - 2Ar" - // (1 - Ar")pl] s', = (1 - Ar'^)p;o4+2/^ (1 - Xr^ ^ p:;^^^ 



m=l 



+2/.(l-Ar") ^ p:;'„4-2/.(l-Ar")<; n = 1, 2, . . . , (5.73) 

m=l, m^n 



where, by (5.69), 



1 1 

Poo = 2 y y /c (x, ^) dxd^; 



1 1 

Pom = 2 y j k{x,^) COS {2m7rx) dxd^; 



1 1 

Pom — 2 y y A;(x,Osin (2m7rx) dxdC; 







5. 6. AN ARBITRARY FUNCTION FROM L2 AS THE SOL UTIONl 19 



1 1 

Pon = 2 / k (x,^) COS (2n7rx) dxd^; 





1 1 



2 J J k{x,0 sin(2n7rx) dxdC; 


1 1 

Pnm = 2 J j k (x,^) COS {2n7rx) cos (2m7r^) dxd^; 

1 1 

— J j k{x,^) cos (2n7rx) sin (2m7r^) dxd^; 


1 1 

Pnm^'^j j k sin (2n7rx) cos (2m7rO o?a;c?^; 


1 1 

= 2 y J k{x,^) sin (2n7ra;) sin (2m7r^) dxd^, n, m = 1, 2, 



(5.74) 

Obviously, to ensure the solvability of the system of equations (5.73), 
the parameter must be such that, as in section 5.4, Eq. (5.60) would 
admit only of the trivial solution. Note that for A = r~", n = 0, 1, . . ., 
that is, in the case when the condition (5.18) is not fulfilled, the ele- 
ments of the column of the free terms (5.73) tend to zero. 

In expression (5.41), 

J H {x, C, A) k (C, dC = 2^i_2X) ^' 



00 j.n 



+ E I _ 2Xrn + "^"^ ' 

and, accordingly, 

-j^ 00 ^n 

V'l (^) = 2 "° ^ E I _ 2Xrn (2n7rx) + sin (2n7rx)] , 
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where 



"0 ^ 77i ^ ^oO"^o + 2 Pom-Sm + Pompom - 2co ; 



4(1- 2A) 



m=l 



(5.75) 

The substitution of this function into (5.54) and subsequent substi- 
tution of q (x) into (5.57) lead to the expression 

(1-2A)A 
^^"^ = - 2(l-2A-P) "° 



n=l ^ 



(1 -2Ar") A 



[a„ (r) cos (2n7ra;) + a'^ (r) sin (2n7ra;)] 



2Ar" - AV2" 
As a result, by formula (5.30), we obtain 

2 oo 

■0 (x) = -to + '^tn COS (2n7rx) + sin (2n7rx) , (5.76) 

2 n=l 

where 

. _ (l-A)A (1 - Ar") A 

i_2A-A2''°' i-2Ar--A2r2-''"^'"^' 

= (1 - Ar") A n = 1 2 

l_2Ar"-AV2-'^"^ ^' i>A---- 

The passage to the limit r ^ 1 yields the following coefficients of 
the series (5.76): 

to = crbo; tn = crbn, t'„ = ab'„, n = 1, 2, . . . , (5.77) 

where, by (5.75), 
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^ m=l 
2 oo 

^ m=l 



(1 - A) 



(1 - 2A) (1 - 2A - A2) 



is a constant. Moreover, 



1-2A 1-2A i-:2A 



jjiX 



jjiX 



As an example clarifying the mechanism of the performed transfor- 
mations, let us consider the determination of the functions (5.54) and 
(5.30): 



q{x) = -\ j h{x,0^'i{Od^ 
-1 

= -A J 1 1 + 2 2 r" [cos {2mTx) cos {2mT^) + sin {2mTx) sin {2mT^)] I cp[ (0 



n=l 



[ cp[ (0 + 2 ^ cos (2n7rx) / ^[ (0 cos (2n7rO 
-1 "=i L -1 



+ sin {2mTx) j (p[ {^) sin (2n7r{) (i{ 
-1 

1 

V'(a;)=x(^) + A / i/(a;,e,A)x(Oo?e 





"=i L 

1 

+ sin {2mTx) j x (0 sin {2mT^) 



cos {2mTx) / X (^) cos (2n7r^) li^ 
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Here, 



-2A J if[ iO d^, -2Ar" J ip[ (0 cos {2n7rO d^, 
-1 -1 



-2Ar"y' <^;({)sin(2n7rOdC 



and 



2(1 - A) }^ ,^^^^_ 2(1-Ar") ^ 
1-2A 



■ / X (0 dC; / X (0 cos (2n7rO d^; 



1 



yx(e)sin (2n7rOc^e 



2 (1 - Ar") 
1 - 2Ar'* 



are the Fourier coefficients of the functions q {x) and ip (x), respectively. 

In other words, in the limit r ^ 1, in the first case there occurs 
a redefinition of the Fourier coefficients by the factor —A {(p'l and q 
are determined on a; G [—1,0) and on x G [0, 1], respectively), whereas 
in the second case the function ip (x) is expressed via x (^) by simple 
multiplication by the factor (1 — A) / (1 — 2A). 

The system of the algebraic equations (5.73) in the hmit r ^ 1 
takes the form 

CXD 

[2 (1 - 2A) - (1 - A) Poo] So = 2/x (1 - A) ^ pomSm+p'omC-'^l^ (1 - ^) Co; 

m=l 

oo 

2 [1 - 2A - (1 - A) Pnn] S„ = // (1 - A) PnoSo+2ll (1 - A) PnmSm 

m=l, m^n 



+2// (1 - A) £ -2/^(1- A) c„; 

m=l 

oo 

2 [1 _ 2A - (1 - A)p:J 4 = (1 - A)p;oSo + 2/. (1 - A) ^ p'^^s. 



m=l 



+2/i(l-A) ^ - 2/i (1 - A) 4 n = l,2,..., (5.78) 

m=l, mT^^n 
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The elements of its matrix predominate on the diagonal owing to 
the component 1 — 2A that does not depend on n and m. As a result, 
for the determination of the coefficients Sq, s'^, contained in (5.66), 
various methods prove to be efficient [9]. 

If k {x, and / (x) are - functions, the corresponding Fourier 
series (5.67), (5.68) converge in the mean. From (5.71) and (5.72), for 
r = 1, we get 

K{x,0 = ^^k{x,0, h{x) = -^^^^f{x), (5.79) 

and the factors contained herein are bounded. Hence, the series ob- 
tained by the substitution of expressions (5.67), (5.68) are analogously 
convergent. 

The solvabihty condition for the system of equations (5.78) is equiv- 
alent to the absence of nontrivial solutions to Eq. (5.60) with the kernel 

and the free term (5.79). At the same time, the class of functions that 
may contain the solution of the Fredholm integral equation of the sec- 
ond kind (5.40) is extended to the space L2. 

Accordingly, the series (5.66) approximating the function ipi (x) con- 
verges in the mean, and by Parseval's relation 

00 

n=l 

From (5.77), it follows that the series (5.76) is analogously convergent: 

00 

n=l 

By the Riesz-Fischer theorem [1] and on the basis of the previous 

consideration, we can conclude that it represents an expansion of the 
L2 - function ip (x) satisfying Eq. (5.1) into a Fourier series in terms of 
the elements {cos {2mix) , sin (2n7r,x)}. 

It should be noted that the parameter < r < 1 plays here an ex- 
clusively important role, because in its absence it would be impossible: 

- to construct the algorithm that lead to Eqs. (5.40) and (5.52); 

- to perform transformations of integrals whose kernels have the 
form of the series (5.19), (5.21), (5.55) and (5.56) that diverge for r — > 1. 
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Thus, the values of to- and in (5.76) are determined by the 
Fourier coefficients of the data of the problem by means of a stable 
procedure of the numerical realization that include the following stages: 

- determination of the parameter A from the condition (5.59) with 
r = 1, i.e., 

A^O, A 7^1, A 7^ 1/2, \y^-l±V2; 

- determination of the parameter fi from the condition that Eq. 
(5.60) with the data (5.79) admit only of the trivial solution; 

- evaluation of the coefficients cq, c^, and poo, Pnm, ■ ■ ■ , Pnm using, 
respectively, formulas (5.70) and (5.74); 

- determination of the coefficients sq, Sn and s[j from the system of 
linear algebraic equations (5.78); 

- evaluation of the coefficients to, tn and t'^ using formulas (5.77). 
The fulfillment of the condition (5.6), after substitution into (5.8) 

of expressions (5.76) with the coefficients (5.77) and 

{x) — -aQ -\- an cos {2n'Kx) + a'^ sin {2n'Kx) , 

2 n=l 

reduces to redefinition of the Fourier coefficients: 

«o ^ ^ ^o; Qin = ^ tn] ^ t^, n= 1,2, .... (5.80) 

Accordingly, a limit procedure with respect to r transforms also (f (x) 
into a L2-function. The condition (5.6) is now understood in the sense 
that 

ll<^/lk(o,i) = llV'-ASV'|L,(o,i)=0. (5.81) 

Thus, for t/j (x) e L2(0, 1), one can find the Fourier coefficients 
of the function if (x) that allow for the fulfillment of the condition 
(5.81). However, this discretization done at the very beginning, i.e., 
without the transformation with the parameter < r < 1, as already 
mentioned, would completely exclude any possibility of the construction 
of the algorithm permitting the determination of the function ijj (x) 
satisfying (5.1). 

In the limit r ^ 1, expression (5.22) turns into it by (5.80) and 
(5.76). Accordingly, the restriction on the form of the free term / (x), 
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imposed by the "harmonic" case of the solution of the problem, is no 
longer in force. 

It is important to note that the conclusion of section 5.5 that the 
function t/j (x) actually satisfies Eq. (5.1) still holds for r — > 1. Relation 
(5.61) is fulfilled in this case by analogy with (5.81), that is, in virtue 
of mutual dependence between the Fourier coefficients of the functions 
X{x) and if' (x), (p[{x). 

In section 6.3, we present a method of the solution of the problem 
(5.1) without proceeding to the limit with respect to the parameter r. 
This is achieved at the expense of satisfaction of the condition (5.6) in 
the sense of generalized functions. The general orientation of transfor- 
mations remain unchanged and the results of section 5.4 will be used 
to a full extent. 



126CHAPTER 5. A METHOD OF THE REDUCTION OF PROBLEMS, TRADITION. 



Bibliography 



[1] F. G. Tricomi, Integral Equations (Dover, New York, 1957) (the 
Russian edition: IIL, Moscow, 1960). 

[2] G. V. Tolstov, Fourier Series (Nauka, Moscow, 1980) (in Russian). 

[3] N. K. Bari, Trigonometric Series (Fizmatgiz, Moscow, 1961) (in 
Russian) . 

[4] A. F. Timan and V. N. Trofimov, An Introduction to the Theory of 
Harmonic Functions (Nauka, Moscow, 1968) (in Russian). 

[5] The approximate Solution of Operator Equations: M. A. Krasnosel- 
skii, G. M. Vinikko, P. P. Zabrciko et ai, Restoration of Dependencies 
from Empirical Data (Nauka, Moscow, 1969) (in Russian). 

[6] V. I. Krylov and L. T. Shulgina, A Handbook on Numerical Inte- 
gration (Nauka, Moscow, 1966) (in Russian). 

[7] Mathematical Analysis (Functions, Limits, Series, Continued Frac- 
tions). Edited by L. A.Liusternik and A. R. Yanopolskii (Fizmatgiz, 
Moscow, 1961) (in Russian). 

[8] N. S. Bakhvalov, N. R Zhidkov, and G. M. Kobelkov, Numerical 
Methods (Nauka, Moscow, 1966) (in Russian). 

[9] D. K. Faddeev and V. N. Faddeeva, Computational Methods of Lin- 
ear Algebra (Fizmatgiz, Moscow, 1963) (in Russian). 



127 



BIBLIOGRAPHY 



Chapter 6 

An analysis of the material of 
the previous section and 
some additions 



6.1 Comments on the material of the sec- 
tions 

In section 5.1, we have developed the previous arguments that the 
restoration of the function (x) from the results of integration of / (x) 
cannot be considered in terms of the solution of the Predholm integral 
equation of the first kind (5.1). Wc have formulated the problem of 
the determination of from the data A and / taking account of an in- 
evitable error of the calculations. To this end, we have proposed to use 
a functional relationship between the error of integration, (5/) (x), and 
■0 (x) in order to compensate adaptively for a small mismatch between 
R{A) and {Aip) {x) that are actually known [sec (5. 2), (5. 3)]. 

Further, it is shown in section 5.2 that a functional model of the 
error of evaluation of the integral (5.1), sec section 4.5, can indeed be 
represented by expression (5.4). The latter is a difference between the 
sought function ip {x) and an integral over this function as well as one 
more unknown function ip (x), with the kernel h{x,^) that has the form 
(5.13). In this case, the fulfillment of (5.6), the condition that reflects 
the smallness of {Sf) {x), requires that the function ip {x) be harmonic. 
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Such an assumption is, apparently, applicable to the type of prob- 
lems that are concerned with the determination of the heat transfer 
(described by the Laplace equation) from the result / (r, x) of its ac- 
tion on a system characterized by k{x,^). At the same time, it is 
desirable that the function ip (x) satisfying Eq. (5.1) be more or less 
arbitrary and, ideally, belong to the space L2. 

Obviously, the above-mentioned harmonicity is stipulated by the 
presence in the expression h {x,^) of the parameter < r < 1. More- 
over, the use, instead of (5.13), of a different, also bounded, kernel, in 
practice, does not yield anything new, because the range of a completely 
continuous operator is not closed.^ 

A rather important point is the extension of (5.4) to a; e [—1,0) 
under the condition (5.6), carried out in section 5.3, which led to Eq. 

(5.24) . In contrast to this equation, Eq. (5.25) is more abstractly 
related to the problem (5.1). This equation arises as a result of the 
suggestion that the efficiency of the transformations will be facilitated 
by the use, together with (5.24), of an analogous equation that is distin- 
guished by its free term going to zero on the other part of the interval of 
definition, x e [—1,0). By means of simple transformations, it proved 
to be possible to obtain the key, in this case, relations, i.e., (5.32), 
further, (5.43), (5.44), (5.48), and, finally, (5.50). 

Equations (5.24) and (5.25) are rather specific. Obviously, on the 
subtraction of 



ip'^{x) = xj h{x,0^o{Od^ + i^(^), ^e[-l,0) (6.1) 
-1 

from Eq. (5.23), on this part of the interval of definition appears Eq. 

(5.25) , and, accordingly, taking into account also (5.44) and (5.48), 

(f' {x) = (p{x) - {x) = if' (x) - (^0 {x) - v^o {x) 

= V?' (x) + (j-) - v?o {^) ■ 

Hence, ipQ = ipQ, i.e., Eqs. (6.1) and (5.36) are identical on a; G [—1,0). 

Simultaneously, the free term of Eq. (5.24), k,{x), "flows" (it is 
difficult to characterize this procedure otherwise) to the free term of 



^Any closed subspace of R{A) is finite-dimensional ([1], p. 96). 



6.1. COMMENTS ON THE MATERIAL OF THE SECTIONS 131 



Eq. (5.25), x{^')- Indeed, by (5.32), equation (5.24) on x G [0,1] 
undergoes the transformation 



-1 -1 



^Xj h{x,0^'{Od^ + x{x). 
-1 

However, whatever one might say about the premises of the con- 
struction of (5.24), (5.25), these equations are, both formally and actu- 
ally, Predholm integral equations of the second kind, whose the solution 
has the form (5.28)-(5.31). On one part of the interval of definition their 
free terms are contained in an explicit form, whereas on the other part 
they enter under the sign of integration. This issue, being absolutely 
nonessential from the point of view of both general theory of this type 
of equations as well as methods of their numerical realization, is a very 
important factor of the realization of further transformations. 

In section 5.3, we have presented a scheme of the construction of 
Eqs. (5.24) and (5.25) starting from a hypothetically given function 
ip (x). In other words, the structure of these equations does not contain 
contradictions. 

A trivial, at the first sight, addition of (5.1) to (5.24) and (5.25), 
which led to Eqs. (5.34), (5.35), has rather substantial meaning of 
embedding the model of the error in the procedure of the determination 
of the function ip (x) . 

Turning to section 5.4, we note that, with the help of (5.45) and 
(5.48), relations (5.43), (5.44) reduce to the following: 

1 

K{x) = xJh{x,Oi^o{Od^, a;e[-l,0); (6.2) 



ip{x) = ip'{x)+ip'o{x), a;e[-l,0). (6.3) 

(It seems that, irrespective of the above reduction, this result is by no 
means obvious.) 
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Let us demonstrate the reduction of (5.34) to Eq. (5.35).^ The 
substitution of ip (x) from (6.3) into (5.34) leads to the equations 

1 

^{x) = xjh {x, [ip' iO + ^0 m d^ + xjh{x,0i^ iO 

-1 

+l^iA'ijj){x)-fifix), xe[0,l]; (6.4) 



^' (x) +ip',{x)^xjh {x, (0 + ^0 m 

-1 

1 

+A y /i {x, ^ iO d^ + K{x), xe [-1, 0) . (6.5) 



Prom (6.3), (5.32) and (5.43), it follows that in (6.4) and (6.5) we 
have, respectively, 



A j h{x,i)v'^{i)di^\j h(x,o b(0-¥''(0]rfe 
-1 -1 

= x{x): xe[0,l]; (6.6) 



p,'{x) + K{x)^\j h{x,Ov'{i)di 
-1 

1 

+\jh (x, V' (0 d^ + i^ix), xe [-1, 0) . (6.7) 



The fact that (6.4), (6.5) are identical to Eq. (5.35) is obvious [the 
function k, is eliminated from (6.7)]. Analogously, vice versa, equation 
(5.35), by the use of relations (6.2), (6.3) and (6.6), is reduced to Eq. 
(5.34). 

Prom (5.39), it follows that 

i/j {x) = iJo {x) + ^1 {x) ; (6.8) 



^This procedure is inverse to that of section 5.4. 
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that is, the function satisfying (5.1) is a sum of the solutions of the 
Fredholm integral equations of the second kind (5.36) and (5.38) that 
are stipulated by the components of the free term of (5.35), i.e., x ^-iid 
—/if, respectively. 

The function ijji (x) depends on the data of the problem and, as such, 
represents the solution of the modified Eq. (5.1), artificially "shifted" 
into the plane of the stability of the procedures of numerical realization. 

This is the solution of a problem that is completely different from 
the considered one, and, quite naturally, the function ipi (x) does not 
satisfy Eq. (5.1). 

In its turn, the function i/jq (x) depends on ipi (x) , which follows 
from Eqs. (5.36), (5.52) and expressions (5.54), (5.41). The addition of 
i/jq and ipi in (6.8) compensates adaptively for the effect of the above- 
mentioned "shift", which makes the function ip {x) satisfy Eq. (5.1). 

Here, it should be emphasized that, at every stage of the solution, 
the transformations, i.e., the "shift" and "compensation for the shift" 
are carried out in association with a well-posed problem. The procedure 
(6.8) can be interpreted as discarding a part of the function ipi (x) that 
prevents satisfaction of Eq. (5.1). 

Let us employ the relation cp — 2(pQ + (p[ that follows from (6.3) and 

(5.39) . Accordingly, 9? — 9?q = <^q -|- (f[ and, in virtue of + ^p'l = (/?', 
we get (p = (f' + ipQ, that is, we return to (6.3).^ This situation is 
completely in line with the logic of the "flow" of the functions k and 
X from one to another. Indeed, by "giving away" ip^, the function ip 
turns into (p', and, instead of k, there appears x- Equation (5.34) takes 
the form (5.35). The inverse procedure, i.e., a transformation of (5.35) 
into (5.34), is, naturally, related to the "return" of ipQ. 

Thus, under the assumption that the function ip (x) is harmonic, 
the problem has been reduced to the solution of Eq. (5.52). Its free 
term depends on the function ip[ (x) that, in turn, is also determined 
by the solution of the Fredholm integral equation of the second kind 

(5.40) and by expression (5.41). 

The above-mentioned results of the transformations (they can be 
characterized as equivalent) should be interpreted in the following way. 

^In other words, the solution of Eq. (5.34) on a; e [—1,0) is the sum of the 
solutions of Eqs. (5.35) and (5.36) on this interval. 
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There is a harmonic fTinction ip {x). After integration according to (5.1), 
it is determined by expression (5.30). The latter, in virtue of < r < 1, 
is a Predholm integral equation of the second kind with respect to the 
function x(x). Under the condition (5.18), its solution, x* i^), can be 
determined in a certain way. From this point of view, the substitution 
of X = X* (^) iiito Eq. (5.52), irrespective of the form of the kernel 
I (x,^), allows us to evaluate the free term q = q* {x). Hence, equation 
(5.52) has every right to exist. 

In other words, for any given kernel /(x,^) that, specifically, has 
the form (5.55) and for a corresponding value of the parameter A, there 
exists a free term q {x) such that the solution of Fredholm integral equa- 
tion of the second kind (5.52), x* (^)) after substitution into expression 
(5.30), allows us to determine the function ip = {x) that satisfies Eq. 
(5.1). 

The above transformations consisted, in essence, both in the de- 
termination of Eq. (5.52) itself and in effective determination its free 
term q. Here, the kernel / (x, ^) docs not depend on the data of the 
problem and is stipulated exclusively by the interests of a constructive 
side of the transformations. Implied is a possibility to make use of the 
techniques of the theory of Predholm integral equations of the second 
kind with symmetric kernels resulting from the model of the error (5.4), 
condition (5.6), the kernel (5.13) and the way of further extension of 
the problem to a; G [—1, 0). 

Carrying out the transformations in an analytical form, including 
finding the resolvent (5.56), was substantially facilitated by the proper- 
ties of the kernel h (x, ^).^ At the same time, for this purpose, instead 
of (5.19), we could use other convergent series in terms of the elements 
<f2n-i (x), <f2n (x) from (5.17). 

However, the kernel (5.13) has an inherent special property that 
consists in the fact that, for r — > 1, the integral 

} I CO 

/ h {x, ^) -0 (0 — T^^o + X] ^" (2n7rx) -|- 1'^ sin {2n7rx) , x e [0, 1] 



list of these properties is given in the next section. 
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[see (5. 19), (5. 76)], where 

1 1 
to = 2 y (e) c^C; tn = 2j\l: (e) cos (2n7re) d^; 



1 

4 = 2y^(Osin(2n7r0^ie, n=l,2,..., 



is a Fourier series of the function ip (x) in terms of the elements (5.17). 
As is known (see, e.g., [2], pp. 110-116), using such a series, one can 
approach in the mean an arbitrary function from the space ^2.^ 

Here, a one-to-one, continuous and linear correspondence between 
the spaces I2 and L2, resulting from the Riesz-Fischer theorem ( [3], pp. 
116-119), manifests itself to a full extent. At the same time, a passage 
to the Hmit r ^ 1 can be regarded as a reahzation of the objective to 
transform I2 into the space I2 (see section 4.5). 

It should be noted that one can draw a conclusion about the stabil- 
ity of the computational procedure of section 5.6 using the passage to 
the limit r — > 1 from the linear dependence of the Fourier coefficients 
to, tn, t'^', sq, Sn, s'„ and Co, c„, of the sought function ijj{x), the 
function ipi (x) satisfying Eq. (5.40) and of the free term / (x) from 
(5.1), respectively [see (5.76), (5.77) and (5.78)]. 

The following point seems to be characteristic. Upon the substitu- 
tion of expressions (5.71) and (5.72) with r — 1, that is, (5.79), equation 
(5.40) docs not change its status as a Fredholm integral equation of the 
second kind. In this regard, it should be noted that the expansion of 
ipi (x) into the series (5.66) is merely one of possible ways of its so- 
lution. If one carries out the same substitution into (5.41), evaluates 
numerically ipi {x) from Eq. (5.40) and, after that, the function (f[ (x), 
the function (x) is determined by means of multiplication by the co- 
efficient 

A(l-A) 
1 - 2A - A2 

[see (5.77)]. At the same time, this fact became clear only as a result 
of the transformations with the parameter r and letting it go to 1. 

^In this sense, an alternative is given by the kernel (5.19) for r = 1, which is the 
series h{x,^) = 1 + 2 X^^i cos [2n7r {x — ^)] whose sum is not bounded. 
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The proof that ijj (x) satisfies (5.f) (see section 5.5) is a very impor- 
tant point whose meaning lies in the following. In the derivation of Eq. 
(5.52), the condition concerning the identity of the solutions of Eqs. 
(5.25) and (5.35), although in an implicit form, has been employed. An 
analysis of the actual transformations has allowed us to draw a conclu- 
sion that this condition is indeed fulfilled and that the function ip (x), 
determined by means of the sohition of (5.52), satisfies Eq. (5.1). 

In this way, we have essentially confirmed a possibility to realize in 
(5.52) the free term q (x) that is adequate to the substitution for (x) 
of a function whose integration by (5.1) yields, as a result, / (x). 



There exist a number of works concerned with the issue of the per- 
turbation of linear operators ([4], [5] section 7, and others). Therein, 
mostly completely continuous perturbations as well as perturbations of 
the spectrum arc studied. The zero error (5.4) is an incompletely con- 
tinuos perturbation. As shown in section 5.4, such a perturbation (in 
contrast to a completely continuous one) can qualitatively change the 
formulation of the problem and introduce principally new possibilities 
of its numerical realization. 

In this regard, condition (5.6) that subsequently terns into (5.81) is 
necessary. Indeed, there arises (5.8), a Fredholm integral equation of 
the second kind with respect to the sought function ip {x), that creates 
the premises of far-reaching transformations. Taken together, equations 
(5.4) and (5.6) can be characterized as the main factor of the construc- 
tion of a stable algorithm of numerical realization of the problem (5.1). 

Nonetheless, the above does not suffice to carry out the transforma- 
tions of Chapter 5. Let us set in (5.4) A = 1 and, instead of (5.5), let 
the operator be 



6.2 Additional arguments 



X 




(6.9) 



-1 



For unique solvability of Eq. (5.7), it is necessary here to have a 
kernel h {x, ^) that possesses the property of being closed. Therefore, 
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it can be taken in the form (5.13). Instead of (5.8), we now have 

X 

g{x) = i;{x)- J h{x,0i^{0dC- 
-1 

Taking into account this point, by extending Eq. (5.7) to x G 
[—1,0), analogously to section 5.3, we obtain 

X 

iP{x) = Jh {x, ^ (0 d^ + Jh {x, (0 d^, X e [0, 1] ; (6.10) 

-1 

X 

V{x)^ jh {x, V (0 di + K{x), xe [-1, 0) , (6.11) 

-1 

where k (x) is an undefined function. 

The solution of Eq. (6.11) is expressed via the resolvent of the kernel 
h{x,^). Its substitution into (6.10) leads to an equation of the form 

X 

iP{x) = Jh{x,0i^{0d^ + x{x), xe[0,l], 



where the function x depends on k. 

However, this equation cannot be related to Eq. (6.11), that is, the 
procedure of extension to x e [—1, 0) does not yield anything in reality. 
The reason lies in the absence of the function ip in Eq. (6.11). If the 
extension of (6.10) to a; G [—1,0) is done with the use of a definite 
integral over ijj, we get the algorithm of section 5.4 in a complicated 
form. 

At the same time, the actual reason for the invalidity of the opera- 
tor (6.9) for application in (5.5) is rooted deeper. The essence lies in a 
qualitative mismatch between the ranges of the Fredholm and Volterra 
integral operators of the first kind. Whereas in the first case the so- 
lution of the corresponding equation exists only under the conditions 
of Picard's theorem, in the second case, it is sufficient for its definition 
that the kernel and the free term be continuous.^ 

^Implied is a reduction to the Volterra integral equation of the second kind by 
differentiation. 



138CHAPTER 6. AN ANALYSIS OF THE MATERIAL OF THE PREVIOUS SECTIOl 



In light of the above, the second factor of the achieved efficiency 
should be noted. It is related essentially with the extension of Eq. 
(5.4), where the operator B has the form (5.5), under the condition 
(5.6), to X e [-1, 0) by (5.23). In this case, the solution of Eq. (5.24) on 
X G [— 1, 0); a; e [0, 1] contains the function k (x) only in an explicit form 
and under the sign of integration, respectively. This point constitutes 
an important prerequisite of obtaining a Fredholm integral equation of 
the second kind for the function n {x) . 

The third factor consists in the use of Eqs. (5.25) and (5.35) along 
with (5.24), (5.34). With the help of these equations, the construction 
of the algorithm moves into the plane of practical realization. In the 
process of the reduction of (5.35) to Eq. (5.34) that has the same form 
of the solution on x e [0, 1], we have obtained the basic computational 
relations. 

And, finally, the fourth factor is related, in fact, to the choice of the 
kernel h {x, ^) that allowed us to do the following: 

- carry out the transformations in an analytical form up to their 
final stage; 

- determine the function t/j (x) for the data of (5.1) from the space 
L2 by means of a passage to the limit in the solution obtained for the 
case when < r < 1. 

In addition, the kernel (5.13) has a whole spectrum of positive prop- 
erties: namely, it is closed, symmetric and positive definite; it depends 
on the difference of the arguments, and the eigenfunctions of the opera- 
tor B are orthogonal both on the interval x e [—1, 1] and on x e [—1, 0); 



Let us turn to the question that is related to Eq. (5.40). For r = 1 
in (5.71), (5.72), we get (5.79). Accordingly, equation (5.40) takes the 
form 



[0,1]. 



1 - A 
1-2A 



1 



1-2A 



/ 



k {x, ^1 (0 



fix), xe[0,l], 







(6.12) 



or 



1-2A 



1 



/.(I -A) 



V'i(x) + |fc(x,OV'i(0^e = /(^), ^e[0,l], (6.13) 
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which makes it rather interesting. As a matter of fact, instead of an ill- 
posed problem, as a basic object of investigation, there actually arises a 
Predholm integral equation of the second kind obtained just by adding 
to (5.1) the sought function with a coefficient whose set of admissible 
values is practically unlimited. 

Indeed, for A ^: 0,1,1/2,-1 ± ^2, it is not difficult to choose 
the parameter /x in such a way that the solution of the homogeneous 
equation (7.12), i.e.. 



where 



be trivial. 

Note that in the process of the evaluation of the function (p\ (x), 
information about the data of the problem contained in ipi (x) under- 
goes substantial changes that involve the kernel and the free term of Eq. 
(5.1). Simultaneously, the next stage of the calculations concerned with 
the determination of q (x) is transferred from x e [— 1, 0) to a; G [0, 1]. 

After that, i.e., in the process of the evaluation of the Fourier co- 
efficients of the functions q{x), x (x) and ip{x), no new information 
about the data of the problem is introduced. As a matter of fact, the 
Fourier coefficients of the function ip[ {x) are triply multiplied by the 
corresponding constants. At the same time, by turning to the system of 
equations (5.78), we can notice that a relationship between the Fourier 
coefficients of the functions ip (x) and ipi (x), i.e., to, tn, t'^ and Sq; •Sn; 
s^, respectively, has, by (5.70) and (5.74), rather substantial meaning. 

Upon the substitution of the function ipi [x) from (5.39) into (6.12), 
taking account of (5.1) and (6.14), we get 



1 




X e [0, 1] , 







1 

V'o {x) ^l^xjk (x, V'o (0 + fo (^) , (6-15) 
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where 

1 

{x) =ij{x)-iixjk {x, V iO + /^a/ (x), xe [0, 1] , 



which leads to rather interesting, as it seems, conclusions: 

1) There exists a Frcdholm integral equation of the second kind with 
the kernel k {x, ^) from (5.1) and the parameter (6.14) whose free term, 
on the subtraction of fi\f {x) , is the same as in a completely identical 
equation for the sought function ip (x). Moreover, 

ipo (x) =ip{x)- ipi (x) , 

where ipi (x) is the solution of (6.12) that is also a Fredholm integral 
equation of the second kind. 

2) And, vice versa, the function ip (x) satisfying (5.1) is expressed 
from (6.15) via the sohition of the Fredholm integral equation of the 
second kind (5.36) and the data of the problem. Note that, as a result 
of the subtraction of Eqs. (5.36) and (6.15), the function (x) can also 
be represented in terms of integral dependence on ipo (x), <^q (x). 

3) For A = 1/2, equation (6.13) turns into (5.1). At the same time, 
for a different value of the parameter A, the solution of this equation is 
a well-posed problem, and, as shown above, it serves for the determina- 
tion of ip (x) by means of a stable procedure of numerical realization. In 
other words, equation (5.1) corresponds to a set of well-posed problems 
for the critical values of the parameters contained therein. 

Thus, the functions ip (x), ipo (x) and ■01 (x) in (5.39) are mutually 
related by means of the Fredholm integral operator (5.5). It seems that 
we have outlined an important point that deserves further interpreta- 
tion. 

The next issue is related to a possibility of the reahzation of other 
methods of the determination of the function x {x) or k. {x) that allow 
us to find the solution of the problem from (5.28), (5.30). We outline 
schematically one of these methods: substitution into (6.2) from (5.30) 
and (5.40) of 

ifjo (x) =iIj{x) - ipi {x) 
1 

^X{x) + \j H{x, e, A) X (0 di - V'l {x) ; (6.16) 
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substitution of k{x) into (5.28); elimination of tp (x) from expressions 
(5.28) and (5.30). 

However, in this situation, the kernel of the integral equation that 
serves for the determination of x (^) depends analytically on the pa- 
rameter A. As a result, there appear unnecessary complications. As 
a matter of fact, generally speaking, such equations may prove to be 
insolvable irrespective of the value of the parameter A ([6], pp. 130-132; 
[7]). In general, an approach to the solution of the problem based on 
the use of (6.16) seems to be less efficient 

An analogue of Eq. (5.52) can be constructed also for the function 
K,{x). To this end, we substitute ipo{x) from (5.39) into (6.2), where 
the function ip (x) has the form (5.28), i.e., 



V^o {x) = xl H{x,^,X)K{Od^-^i (x). 
-1 

As a result, 



K(x)^Ajl'(x,Oi^(Od^ + (l'(x), xe[-l,0), (6.17) 
-1 

where A = A^; 

1 

I' {x,0= J h{x,OH{c,tx)dC. 



A comparison with (5.53), by (5.19), (5.21), shows that I' (x,^) = 

1 

q'{x) = -Xjh{x,Oi^i{OdC- (6.18) 



Compared to the algorithm of section 5.4, in this case, the necessity 
of intermediate determination of the function (pg (x) drops out, which, 
by the way, does not simplify substantially the transformations. 

Note that if the free term / (x) has discontinuities or other sin- 
gularities stipulated by the kernel k{x,^), its explicit presence in the 
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solution may prove to be desirable. To this end, the function ip (x) 
should be expressed via k{x) with the help of Eq. (5.34). As regards 
Eq. (5.24), it is satisfied as a result of the use of (5.31) in the process 
of the construction of Eq. (5.52).^ 

Wc want to conclude this section by returning to the condition (5.3). 
Let us draw attention to an adaptive connection of / (x) with the space 
I2 and the logic of its practical realization. Preliminarily, we note main 
facts concerning the mapping (5.1) from one of the spaces L2, I2 into 
another, irrespective of supposed use of 5f. 

Thus, the range R{A) — is not closed. As a result, the operator 
A is incompletely continuos, and its inverse A^^ acting from V2 into 
L2 (0, 1) is bounded. However, in the framework of the traditional 
object of investigation, (5.1), there is no possibility to use somehow 
this fact for the construction of the operator A'^. 

At the same time, the following sequence of arguments arises: 

- objectively, a bounded operator A"^ from ^2 does exist, that is, 
the solution of Eq. (5.1) in the pair of spaces (1/2,^2) is a well-posed 
problem; 

- the property of limited inversion is directly associated with the 
Predholm integral operator of the second kind, which in our case is 
I-\B] 

- insertion in the scheme of transformation of the identity operator 
/ organically combines with the modelling of the error of integration 
induced [together with the prime cause, i.e., non-closed character of 
it! (^4)] by the incorrectness of the problem (5.1); 

- the adaptation of / {x) to the space ^2 ^-nd a functional represen- 
tation of the error are, thus, closely related; 

- the fact that R {A) is non-closed does not prevent the determi- 
nation of the function ipi [x) from Eq. (6.12). As shown above, this 
function serves for finding the solution of the problem, the function 
il}{x)] 

- from this point of view, equation (5.7) that follows from (5.4), 
(5.6) realizes a connection of the Predholm integral equations of the 
first and the second kind via their common range, which creates a 
serious prerequisite of the solution of the problem of the determination 



''As a matter of fact, this issue is discussed in section 5.5. 
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of the function ip (x) satisfying (5.1) in the correct formulation. 
The outlined orientation will be developed in the next section. 



6.3 The second version of the solution of 
the problem 

Here, the transformations of sections 5.2-5.4 are extended to the case 
when the data of Eq. (5.1), i.e., k {x, ^) and / {x) as well as the func- 
tion ip (x) satisfying this equation, from the very beginning belong to 
the space L2. The kernel h{x,$,) has the previous form (5.13); the 
parameter < r < 1 is fixed. ^ 
Let us return to Eq. (5.24): 

-XB(^']ix) + l 0' (6.19) 

\ (p p ' [ ii{x), xe [-1,0), 




or 



U 1 

i;{x)^\jh{x,i)^{i)di + \jh{x,i)i:{i)di, xe[0,l]; (6.20) 



^{x)^xjh (x, V (0 d^+^ j h {x, V' (0 d^+K (x) , xe [-1, 0) , 
-1 

(6.21) 

which follows from the representation of the error (5.4), under the con- 
dition (5.6) and the extension of (6.20) to x & [—1, 0). 

A possibility to satisfy this equation on the interval x e [0, 1] by 
(p{x) was considered in section 5.2 [the function ip {x) was given]. In 
order to satisfy the condition (5.6) that implied fitting in the space 
C [0, 1], we had to restrict the class of admissible functions for ip [x) by 
harmonic functions. 

In this case, the issue essentially reduced to an investigation into 
the solvability of the Fredholm integral equation of the first kind (5.10) 



^Below, we point out that, in the framework of the present version of the solution, 
expression (5.13) is, generally speaking, has alternatives. 
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stipulated by the conditions of Picard's theorem [2]: 

oo 27r 

Y: alXl < oo, [g (C) <fin (C) dC (6.22) 

n=l ^ 

where A„ and {() are, respectively, the characteristic numbers and 
the eigenfunctions of the kernel (5.11) numbered in order of their se- 
quence order [see (5.12)]. Moreover, the system of the eigenfunctions 
(fn (C) must be complete on the interval ( G [0, 27r], and the kernel must 
be real and symmetric, which, in this case, is certainly fulfilled. 

And, nevertheless, as shown by E. Goursat ([8], pp. 141-143), even if 
the condition (6.22) is not fulfilled, one can always find such a function 
ip {Q that the difference between the integral 

J h{c,e)ip{e)de 



and the function g {() from (5.10) is arbitrarily small. It should be 
noted that h{(,9), in this case, is assumed to be a much more general 
kernel than (5.11). 

A proof is based on the fact that 

9n iO = Jh (C, 0) (6) de, 



where 

n 

^(«)(^) = ^a,\^,(^), (6.23) 

i=l 

coincides with the sum of the first terms of the Fourier series of the 
functions g {() in terms of the elements ipn{d)- Therefore, one can 
establish the number n for which the integral 

27r 


is smaller than a certain £ > 0. 
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However, the series (6.23) that satisfies in this way Eq. (5.10) di- 
verges in the space L2 with the growth of n. Accordingly, it is not 
possible to regard (5.19) as a Predholm integral equation of the second 
kind for the function 

and to invert the operator / — \B. In terms of the terminology that 
became predominant later, 93 (x) is interpreted as a generalized function 
(distribution) ([9], section 2.1.5; [10], section 12). 

Thus, there exists a function (p{x) satisfying Eq. (6.20), or (5.7), in 
the sense that 

jdxi^jhix, io - A (B^) m c^ej - 0, 

where the operator B has the form (5.5); ijj (x) = (p{x), x G [—1,0) 
and the change of variables (5.9) is made. 

At the same time, the kernel (5.13) that is actually used is infinitely 
differentiable and depends periodically on ^, which allows us to inter- 
pret (f (x) as a generalized function in a less restrictive sense of the 
convergence of the series (6.23) ([11], pp. 17-18): 


Jim J h {x, (^("^ (0 d^^Jh {x, ^ (0 d^, X e [0, 1] (6.25) 
-1 



(the variable x plays the role of a parameter). 

Indeed, the substitution of expressions (6.23) and (5.19) into (6.25), 
with the use of (5.17) and under a redefinition of the coefficients, yields 
the following: 

lim J 1 1 + 2 ^ r* [cos {2mx) cos i^m^) + sin {2m x) sin (2i7r^)] | 
+ E COS (2i7rO + a% sin (2i7rO] i d^ 
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1 " 
= -Oq + lim ^ (i2i-i cos (2i7rx) + a'g'j sin (2i7rx) 



n— »oo • 

i=l 



— -Oo + ^ a2n-i COS (2n7ra;) + sin (2n7ra;) , a; e [0, 1] . 

2 n=l 

Further, the function k (x) will be determined from an equation 
constructed on the basis of relation (6.2).^ By (5.39), it takes the form 

1 

K{x)^xJh{x,ObP{0-i^iiO]dC, ^e[-l,0), (6.26) 



where ipi (x) is the solution of the Predholm integral equation of the 
second kind (5.40). 

In this case, it is by no means necessary to substitute the function 
^ (x) from (5.28) into (6.26).^° It is sufficient to express the integral 

1 







via the function k (x), which is equivalent to interpreting also ip (x) as 
a generalized function. 

Realizing objectively in this sense ip {x), we integrate with the kernel 
h{x,^) of Eqs. (6.20), (6.21) in the limits 0, 1 and —1,0, respectively. 
We get: 



1 



i^{x)^\Jh{x,0^{OdC + ^Jh{x,Oi^{OdC, ^e[0,l]; (6.27) 









0{x)^\jh (x, (0 di+X I h (x, V' (0 d^+k (x), xe [-1, 0) , 

-1 -1 



(6.28) 



^This is in contrast to section 5.4, where for an analogous purpose relation (5.32) 
was used, which led to an equation for the function x (a;)- 

^°In this way, we would arrive at Eq. (6.17) obtained under the assumption that 
the function ip {x) is harmonic. 
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1 

-1 


k{x) = J h{x,0f^{0dt (6.29) 
-1 

Given that 

1 

jh{x,i)^{i)di^ j h{x,i)^{i)d^- 

-1 
1 

j h {x, ^{Od^^Jh {x, ^ (0 d^, X e [-1, 0) ; [0, 1] , 

-1 

by 

1 

J h(x,0h(C,0dC^ J h(x,Oh(C,OdC^h(x,0 

-1 

[see (5.19)], equations (6.27) and (6.28) are equivalent to the following 



ones: 



u ± 

^{x)=xjh {x, O0iOd^ + xJh {x, ^ (0 d^, X G [0, 1] ; 



-1 



1 

0{x) = \Jh {x, iO d^ + \ J h {x, ^ (^) d^+k (x) , a; e [-1, 0) , 

-1 

or 

* (X) = A {B^) . A / Mx, * (0 d« + { , _ ; _ 

(6.30) 



-1 



where 



U8CHAPTER 6. AN ANALYSIS OF THE MATERIAL OF THE PREVIOUS SECTIOl 



[Note that the choice of the kernel h {x. ^) again played a positive role.] 
Thus we have obtained an exact analogue of Eq. (6.19), where 
the generahzed function (6.31) substitutes for the function (6.24). The 
inversion of the operator I — XB in (6.30) yields 







[an analogue of (5.28)]. As a result, expression (6.26), by (6.29) and 
(6.18), takes the form 



k{x)^X^ j H {x, e. A) k {i) di + (x) , a; e [-1, 0) . 
-1 

The substitution of k {x) from (6.29), in virtue of 



j H{x,C,X)h{C,OdC = J h{x,OH{C,^,X)dC 
-1 -1 

[see (5.53), (5.19) and (5.21)], leads to the same equation (6.17). 

Thus, the use of generalized functions, in this case, has ensured the 
legitimacy of the transformations and has not affected the final result. 

Further, let us turn to Eqs. (5.34) and (5.35): 

ip{x) /"''''U i^{x), a;e[-l,0), 

(6.32) 

V'(a^) 1_xr/^ '^\(^^Al^{M){x)-^lf{x)+x{x), x G [0, 1] ; 
^'{x) / - 1^ J ^^^+1 0, xe[-l,0). 

(6.33) 

where B is the operator (5.5). 

As shown above, the transformation of (6.33) into Eq. (6.32), or vice 
versa, with the help of relations (5.32), (6.2), (6.3) and others allows 
us to obtain the Fredholm integral equations of the second kind (5.52), 
(6.17) whose solutions, for corresponding values of the parameters yU, A 
exist and are unique. For a specific choice of the kernel h{x,^), either 
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in the form (5.13) or in any other one, they depend exclusively on the 
data of the problem. 

In this sense, the solutions of Eqs. (6.32), (6.33) exist and are unique 
under the assumption that their free terms including the functions 
X (x) arc given. Indeed, the inversion of the operator I — \B in (6.32) 
leads to a Fredholm integral equation of the second kind:^^ 



^l^{x)^l,l K{x,O^I^{OdC + F{x), [0,1]. 



(6.34) 



Here, 



1 

K {x, = k{x,0 + xjH {x, C, A) k (C, dC- 



(6.35) 



F{x)=X J H{x,^,X)K{Od^- 
-1 



-fx 



± 

f{x) + J H{x,^,X)f{Od^ 



(6.36) 



where H{x,^,X) is the resolvent (5.21). This equation differs from 
(5.40) only by a component of the free term, and its solution can be 
found analogously. If ip (x) is known, we can obtain 



^(x)^k(x)+X I H(x,^,X)K(Od^ 



/xA / H {x, e. A) [{Aij) (0 - / m dC 



[the same can be easily done with the use of Eq. (6.33)]. 



11 



They are mentioned in the Introduction. 
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The function k, (x) entering expression (6.36) can be found from Eq. 
(6.17), whose solution is 



K{x)^q'{x) + Aj L{x,^,A)q'{Od^, x e [-1,0) , (6.37) 
-1 

where A = A^; L{x,^,A) is the resolvent (5.56); the function q' (x) is 
determined by expression (6.18). 

The latter depends on the function ipi (x) that, in its turn, repre- 
sents the solution of the Predholm integral equation of the second kind 
(5.40). In this case, it is assumed that the conditions on the parameter 
A, (5.59), are fulfilled. The value of is chosen from the condition that 
Eq. (5.60) possess no nontrivial solutions. 

Provided the function ipi (x) is found, a general sequence of compu- 
tational procedures consists in the determination of the following: 

- the function q' by (6.18); 

- the function k by (6.37);^^ 

- the kernel and the free term of Eq. (6.34) by (6.35) and (6.36) ; 

- the sought function ip by (6.34). 

A proof that the so obtained solution satisfies (5.1) is analogous to 
that of section 5.5. 

Let us touch on numerical realization of Eqs. (5.40) and (6.34). As 
is known, there are a number of stable algorithms for the solution of 
the Predholm integral equation of the second kind, including both the 
evaluation of spectral characteristics and of quadratures. Along with 
the handbook referred to in Chapter 5 as [5], they are presented in [12], 
[13] and in a number of other books. Generally speaking, analogous 
approaches are discussed. 

We can point out S. G. Mikhlin's algorithm of numerical realization 
of the resolvent ([14], section 12) based on a finite-element approxima- 
tion of the kernel of the integral equation. Of special interest is the 
method of G. N. Polozhij [15] that transforms the Predholm integral 
equation of the second kind in such a way that its solution, irrespective 
of the value of the parameter, is achieved by means of simple iterations. 

^^Note that, irrespective of the data of (5.1), the function k is infinitely 
differentiable. 
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In this regard, the initial approximation is established and the second 
iterated kernel is used. 

For convergence in the mean of simple iterations to the solution of 
Eq. (6.34), it is necessary that 

Ci < 1, 

where 

1 1 

ci^l I \K(x,Ofdxd^. 



If 

1 

I \Kix,0fdC<C2, 



where C2 is a constant and condition (6.38) is fulfilled, corresponding 
Neumann's series is absolutely and uniformly convergent [16]. 

Note that a possibility of varying the parameter < r < 1 creates 
an additional reserve of increasing the efficiency of the procedure of 
numerical realization. 

Instead of (5.13), a different kernel h{x,$^) could be used in the 
transformations of the present section. At the same time, the above- 
mentioned numerous advantages of Poisson's kernel make such a change 
absolutely unnecessary. Indeed, if h{x,^) satisfies the conditions of 
Mercer's theorem, which would be unreasonable to give up, then ([2], 
p. 166) 

oo 
n=l 

and, as characteristic numbers of an alternative kernel, one could ac- 
cept, for example, the inverse of the terms of an arithmetic (or geomet- 
ric) progression. Such a preference is rather problematic, whereas some 
of the available advantages could be lost. 

Let us now discuss the issue of the adaptation of Eq. (5.1) to the 
space I2 [defined by the condition (6.22)], which was more than once 
mentioned above. Thus, R{A), the range of the operator A, is not 
closed, which is the origin of the difficulties of the determination of 
the function ip {x) satisfying Eq. (5.1) in terms of the solution of the 
Fredholm integral equation of the first kind. 



(6.38) 
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As regards (6.32), this issue is unimportant, because, by inverting 
the operator / — \B, we get a Fredholm integral equation of the second 
kind, namely, (6.34). At the same time, equation (6.32) is constructed 
by a very interesting subject, namely, (5.7). This is a Fredholm integral 
equation of both the first and the second kind for the functions (p (x) 
and ip (,t) respectively. 

However, as regards the first of the above-mentioned properties, this 
equation is principally different from (5.1) in the following respect: the 
free term g (x), determined by expression (5.8), depends on the sought 
function i/j (x) while not taking on any concrete values. Therefore, one 
can assume that 

g (x) = (x) -jh {x, (0 d^eR {B') , 



where the operator is given by 



B: = j h{x,^),d^, xe[0,l]; 
-1 

and it seems that an element of the above-mentioned adaptivity is ob- 
vious. 

After that, as a result of the extension of (5.7) to a; G [—1, 0), there 
appeared the Fredholm integral equation of the second kind (5.24) for 
^ (x) [see (6.24)] with the function k (x) representing its free term. 
Indeed, for a given k, (x), by inverting the operator / — XB, one can de- 
termine the function ip (x) as well as (/? (x). Note, however: the problem 
of finding the function ip (x) , contained under the sigh of integration 
in (5.1), has transformed into the problem of the determination of the 
function k (x) that enters explicitly! 

Is it possible not to return to the arguments of J. Hadamard for the 
existence of correct formulation of physically substantial problems as 
well as to our repeated suggestion to consider (5.1) as a rule for carrying 
out the integration of the function -0 {x)l As a matter of fact, being 
substituted into a certain Fredholm integral equation of the second 
kind, this function generates a corresponding free term. Consequently, 
the problem reduces to a procedure of its determination, which has 



6.3. THE SECOND VERSION OF THE SOLUTION OF THE PROBLEM153 

a powerful resource, namely, a possibility of an arbitrary choice of the 
kernel of the integral equation (in reality, constructed in a certain way). 

We may assume that here we realize only one of a number of existing 
approaches of the outlined orientation. Given such a favorable object 
as the Predholm integral equation of the second kind and the fact that 
a free term that allows the function ip (x) to satisfy this equation exists 
objectively, there is no alternative to correct formulation of problems 
of mathematical physics! 

Returning to the algorithm, let us follow the way this formulation is 
actually realized. Thus, we have found the function ipi (x), part of the 
solution (6.34), stipulated by —/if {x). From a computational point of 
view, the properties of the employed Eq. (5.40) are excellent. 

Elimination from (6.33) of the component of the solution depending 
on X (x) has led to relation (6.2). With the help of this relation, by 
means of the resolvent, the function k (x) has been determined. This 
function is exactly the free term to whose determination the problem is 
reduced. In this regard, also taking account of (6.37), let us again draw 
attention to (5.28), (5.30), a representation of ip{x), which is, respec- 
tively, "integral" and with the function x (x) entering explicitly. Their 
variation in the process of transformations ensuring the realization of 
a stable procedure of the evaluation of the function k (x) has been one 
of the decisive factors. 

Thus, the problem (5.1) has turned into the one of the solution 
of the Predholm integral equation of the second kind (6.34). By the 
subtraction of (5.40) from this equation, we get 



1 




X e [0, 1] , 



(6.39) 







where 








-1 



the function k{x) is determined by expressions (6.37), (6.18). Indeed, 
ipo {x) is the same function that enters (5.36). This follows from Eqs. 
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(5.34), (5.35) and the representations of the solution (5.28), (5.31), i.e., 

1 
-1 

Thus the function ip (x) can be determined by adding the solutions 
of (5.40) and (6.39): 

ip (x) = -00 (x) +ipi{x), xe [0, 1] . 

After the determination of the function ipi (x) , the problem can be 
solved also by using the substitution of expression (6.26) into (5.28). 
In this case, the function ip (x) satisfies the Predholm integral equation 
of the second kind 

1 

i;{x) = X'Jl {x, (0 + /' {x) , (6.40) 



where 

1 

f{x)^-xll{x,Oi^{OdC 



[see also (5.55)]. 

Accordingly, the solution of (6.40) is expressed via the resolvent 
(5.56): 

1 

ij{x)=f{x)+AjL{x,tA)f{OdC, A^A^; 



and, as can be noticed, in this case, only one Predholm integral equation 
of the second kind is solved numerically, namely, (5.40). After that, only 
procedures of integration are carried out. 

Equations (5.40) and (6.39), or (5.40) with the use of (6.40) whose 
data are stipulated above, embody the problem (5.1) in its correct 
formulation! 

By comparing the algorithms of sections 5.4 and 5.6, we find that 
the latter one is, obviously, more formalized, which may prove to be, 
in a sense, more advantageous. 
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6.4 A summary of computational relations 
(to section 6.3) 

As the necessary formulas are scattered throughout the text, we find it 
reasonable to present them in a consecrative and the most convenient 
for computation form. 

Under the assumption that / (x) G R (A) and the kernel k {x, ^) is 
closed, the function satisfying the equation 

1 

{Aij){x)^ fk{x,^,X)ij{Od^ + f{x), xe[0,l] 



is defined as 

(x) = (x) + ^1 (x) , (6.41) 
where ipo (x) and ipi (x) are the solutions of the Predholm integral equa- 



tion of the second kind 

1 



i^{x)^yij K{x,0ip{0d^ + F{x), x e [0, 1] , 



(6.42) 



with the free term, respectively 



u 

F{x) = Fo{x) = X J H{x,^,X)K{Od^; 
-1 



F(x) = Fi(x) = -// 



u 

f{x)+X I H{x,^,X)f{Od^ 



[in relation to Eq. (5.40), Fi = /i]. 
The kernel of Eq. (6.42) is given by 

1 

K {x, 0=k {x, i) + xj H (x, C, A) k (C, dC 



Here and above, 

1 



H{x,^,X) 



oo 



1-2A 



n=l ^ 



2Xr' 



(6.43) 



(6.44) 



(6.45) 



cos [2n7r (x - 0] , (6.46) 
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with the parameter < r < 1; 



k{x)^p{x)+aJ L {x, e, A) p (0 (6.47) 



-1 



[in relation to Eq. (6.18), p = q']- 
In this expression, A = A^; 

1 

p{x)^-Xjh (x, V'l (0 di, (6.48) 



where 

1 

h{x,i)^ 



1 — 2r cos [2mT (x — ^)] + 

oo 

= 1 + 2 ^ r" COS [2n7r {x - 0] ; (6.49) 



n=l 

h 2 V 

1 - 2A - A ^1 1 - 2Ar" - Ar^" 

(6.50) 

Parameter 



L A) = , , + 2 g ^ COS [2n7r (x - 0] ■ 



A ^: 0; r-"; ^r""; (-1 ± x/2) r"", n = 1, 2, . . . ; 

parameter 

n=l,2,..., 

where Pn ^-re characteristic numbers of the homogeneous equation 

1 

ij{x)=pj K{x,0i^{0d^, xe[0,l]. 



When the values of p, A and r arc chosen (they can be corrected 
later, depending on various arguments), the sequence of computational 
procedures is as follows: 

- determination of the kernels K of Eq. (6.42) from (6.45), using 
expression (6.46); 
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- determination of the free term Fi from (6.44); 

- determination of the function tjji from (6.42), with F = Fi] 

- determination of the function p from (6.48), using expression 

(6.49) ; 

- determination of the function k from (6.47), using expression 

(6.50) ; 

- determination of the free term Fq from (6.43); 

- determination of the function from (6.42), with F = Fq; 

- determination of the sought function ip from (6.41). 
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Chapter 7 



A reduction of linear 
boundary- value and 
initial-boundary- value 
problems to Predholm 
integral equations of the first 
kind 

7.1 Ordinary differential equations 



Consider, for example. 



u" — a{x)u = f (x) , a; e [0, 1] ; 



(7.1) 



u' (0) = (1) = 0, 

where a (x) and / (x) are given L2-functions. 
Prom the notation 



(7.2) 



u" {x) = ip {x) , 



(7.3) 



it follows: 




(7.4) 
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u 



X 

{x) = y (x - V' (0 + cix + Co, (7.5) 



where cq, ci are the constants of integration. 

The substitution of expressions (7.3) and (7.5) into (7.1) leads to a 
Volterra integral equations of the second kind: 

X 

iP{x)^a (x) Jix-C)ip (0 dC + {cix + Co) a{x) + f{x), xe [0, 1] , 



(7.6) 

whose the solution is 



^ (x) = {c^x + Co) a{x) + f{x) + J Q {x, ^ [(ciC + co) a (C) + / (C)] d^, 



(7.7) 

where Q {x, C,) is the resolvent of the kernel a (x) {x — C,)- 

Taking into account (7.4), (7.5) and (7.7), we find from the bound- 
ary conditions (7.2): ci = 0; 



Co = 



/(l-o 




f{0 + IQ{^X)f{C)dC 




di 


i + /(i-0 




a(0+IQ(^,0a(0dC 




di 



(7.8) 



One can act in a different way: Namely, upon the substitution of 
expressions (7.4), (7.5) into (7.2), we get ci = 0; 



Co 



± 



and, as a result, 

u {x) 



X -L 

/(^-0-/(i-0 



^ (0 di. 



(7.9) 
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In contrast to (7.6), the problem reduces to the Fredholm integral 
equation of the second kind 



ip {x) = a {x) 



1 



/ (^-0-/(1-0 



LO 



(7.10) 



whose the solution is 



1 

^{x)^fix)+lQ{x,OfiOdC, (7.11) 



where Q (x,^) is the resolvent of the kernel^ 

The substitution of (7.7) into (7.5), taking into account (7.8), or the 
substitution of (7.11) into (7.9), allows us to find the solution of the 
problem (7.1), (7.2). Note that the outlined approach is substantially 
indifferent to the order of the differential equations, the form of initial 
or boundary conditions and to the data of the problem. 

Analogous transformations are traditionally discussed in courses of 
the theory of integral equations (see, e.g., [1, 2]). At the same time, as 
far as the solution of applied problems is concerned, the construction of 
integral equations of the second kind did not gain sufficient popularity, 
which can be characterized as a kind of a paradox. It is rather surpris- 
ing in light of rather active attempts of its popularization: see, e.g., 
publications of S. E. Mikeladze, I. A. Birger, and A. N. Golubentsev 
[3, 4, 5]. 

It seems that the reasons for this situation are, on the one hand, 
inefficiency of technical means of numerical realization of integral equa- 
tions before wide-spread implementation of computers, and, on the 
other hand, insufficient popularity of the techniques of the theory of 
integral equations among specialists in applied science. 



^It is assumed that the homogeneous equation (7.10) has only a trivial solution. 
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Nonetheless, here is an opinion of G. Wiarda ([6], p. 5): "... an inte- 
gral equation substitutes for a corresponding differential equation with 
its boundary conditions that, as far as a concrete physical phenomenon 
is concerned, necessarily arise with any differential equation. An inte- 
gral equation already contains all the elements specifying the physical 
problem. One more advantage of integral equations lie in the fact that, 
in most cases, we arrive at equations of the same type..., whereas the 
types of differential equations, even in closely related problems, often 
turn out to be rather different." 



7.2 An illustration of the procedure of re- 
duction 

Let us turn to the problem of bending of a membrane stretched along 
a contour by a uniform load: 

+ (7.12) 

«(0,t/) = u(l,y) = 0, (7.13) 
M (a;, 0) = M (x, 1) = 0. (7.14) 



From the notation 



dlu{x,y)^'il;{x,y), (7.15) 



it follows: 



u 



(x, y) = / (x - (e, y) + x^ii {y) + gi2 (y) , (7.16) 



where gij (y) are functions of integration. 

In view of (7.15), equation (7.12) takes the form 

d^u {x,y) = -l-'i/j{x,y), 

and, respectively. 



1 ^ 

u {x, y) = --y^ - j (y-v)'^ {x, v) drj + yg2i (x) + 922 (x) , (7.17) 
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where g2j (y) are also functions of integration. 

The substitution of expressions (7.16), (7.17) into the boundary 
conditions (7.13) and (7.14), respectively, allows us to determine gi2 — 
922 = 0; 

yii(y) = -/ {l-0^{ty)d^; 92i{x) = 2 +/ {^-v)'^l^{x:V)dv- 



As a result, expressions (7.16) and (7.17), respectively, take the form 



u {x, y) 



X i 

.0 



uix,y) = -y{l-y) 



^P{U)di; (7.18) 
ip{x,r])dr]. (7.19) 



J {y-v)-y j (1 - ?7) 

.0 

Eliminating u from these expressions, we get a Predholm integral 
equation of the first kind: 

1 



J{x-0-xJ{l-0 



^{^,y)d^ 



+ 



J {y-v)-y j {'^--v) 



^ix,v)dv ^^yi'^-y)- 



(7.20) 



Thus, a principal difference from the one-dimensional case consists 
in the reduction of the problem (7. 12)- (7. 14) to an ill-posed one. How- 
ever, here we will be interested not in the determination of the function 
■0 satisfying Eq. (7.20) (just note that the algorithms of sections 5.4, 
5.6 and 6.3 apply to it as well) but in the universality of the procedure 
of transformation. 

Indeed, let the domain of the problem be different from the canonical 
one, and let, for example the second condition (7.13) have the form 
(T) y) — 0) where x — ^ (y) a certain single- valued function. Instead 
of (7.18), we have 



u{x,y) = 



y{y) 



j{x-i)-x I [7(2/) -e] 



^{^,y)dC, 
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and, from a computational point of view, any differences are absent. For 
the transition to an ordinary procedure of the evaluation of the integral 
on a rectangular domain, it is sufficient to employ a non-orthogonal 
mapping of the type x — ^x, y — y. 

It is not difficult to notice that each of expressions (7.18) and (7.19) 
satisfy identically the pair of boundary conditions (7.13) and (7.14), 
respectively. The rest of the conditions are fulfilled approximately, 
depending on the accuracy of the determination of ■0. At the same time, 
the solution can be represented in the form that satisfies identically 
both the conditions (7.13) and (7.14): 

Ui (x, y) = ui (x, y)- {1- y) ui (x, 0) - yui (x, 1) ; 

U2 {x, y) = U2 (x, y)-{l- x) U2 (0, y) - xu2 (1, y) . 

Here, the functions Ui, U2 are determined by (7.18) and (7.19), 
respectively. 

The norm of the error of closure of the values of Ui {x, y) or Ui (x, y) 
allows us to estimate the error of the approximate solution: 




d^u (0,1/) = d^u (1,1/) = 



would be imposed, they could not be satisfied by the expression for the 
derivative 



X 








that follows from (7.15). 

Nevertheless, this complication can be easily overcome by the use, 
in particular, of the relation 



diu + (3u = ip^ 



where /? is a constant that allows us to retain both the functions of 
integration g^j {y). 
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Let us turn to an equivalent formulation of the problem (7.12)- 
(7.14): 

dlui + dyU2 = -1; ui{x,y) = U2{x,y) , (7.21) 

m (0, y) = u, (1, y) = U2 {x, 0) = {x, 1) = 0, (7.22) 
using a representation of the solution of the type 

ui {x, y)= ki {x, y, ^1 (C, y)d^ + Yl /^li (^) 9ij iy) ; 
J=i 

y 2 
U2 {x, y)^ k2 {x, y, r]) V'2 {x, ri)dri + J2 f^2j (y) g2j (x) . 
i=i 
We assume that the kernels are given and satisfy the conditions 

ki {x,y,x) = kx {x,y,y) = 0; 

d^ki{x,y,x) ^0; dyk2 ix,y,y) ^ 0, x,ye[0,l]; (7.23) 

fiij (x), fi2j (y) ai'e also given; gij (y), g2j (x) are to be determined from 
the boundary conditions as discussed above. 

Let us set /lu — x, ^21 — y, A*i2 = A*22 = 1- In this case, under the 
conditions (7.22), we get: 



ui {x,y) 



X 1 

j h {x,y,^)-x J ki{l,y,^) 
.0 

■/M..,)-./M.M) 



i^UC,y)dC; (7.24) 



iP2{x,v)dv, (7.25) 



U2 {x, y) 
and, respectively, 

X 

dxh {x, y, x) ipi (x, y) + j dlki {x, y, ipi (C, y) = <9^Mi {x, y) ; 

° (7.26) 



(7.27) 
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Let, in addition to the conditions (7.23), d^ki y, ^) and dyk2 {x^ y, r}) 
be L2-kernels. Here, (7.26), (7.27) are Volterra integral equations of the 
second kind with respect to the functions -01, ^/'2, whose solutions, by 
general theory, exist and are unique. Therefore, the representations 
(7.24) and (7.25) correspond to the physical content of the problem 
(7.21), (7.22). 

In (7.23), we can set 

ki {x, y, ^)^{x- K {^^ ; h {x, y, r))^{y- rj) k'^ {x, y, rj) , 
where 

K {x, y, x) ^ 0; k'^ {x, y, y) 7^ 0, x,y e [0, 1] , 

using these expressions to refract the a priori information about the 
solution in order to smooth the sought functions ipi and, in general, 
to simplify the procedure of calculations. It is clear that this point 
is important for more complicated problems with different kinds of 
singularities of the behavior of the solutions, and we just outline it 
here. 

The substitution of d^ui and S^wa from (7.26), (7.27) into (7.21) 
leads to a system of integral equations 

y 



{x, y) 



1 



dyk2 {x, y, y) 



dyk2 {x, y, rj) -02 {x, r))dr) + F {x, y, , 







(7.28) 



where 



F{x,y,ipi) 



dyk2 {x,y,y) 

X 

+ j dlki (a;,!/, 0^1 {^:y)d^ 


X 1 

j h {x,y,^)-x j ki{l,y,^) 



[1 + d^ki (x, y, x) i/ji (x, y) 



M...,)-./M..M) 



ip2 {x, 77) d?7 = 0, x,y e [0, 1] . 

(7.29) 
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Prom Eq. (7.28), we find 

y 

^2 {x, y)^F (x, y, ^i) + j Q {x, y, rj) F {x, r], ^i) dr], (7.30) 



where Q {x, y, rj) is tiie resolvent of tlie kernel —dyk2 {x, y, rj) /dyk2 {x, y, y). 

The substitution of expression (7.30) into (7.29) allows us to obtain 
a Fredholm integral equation of the first kind with respect to the func- 
tion -01 (x, y). Clearly, the above reduction scheme is more cumbersome 
compared to that based on the formulation of the problem in the stan- 
dard interpretation (7.12)-(7.14). At the same time, one may discern 
in it some iteration elements that result from the fact that (7.28) is a 
Volterra integral equation of the second kind with respect to both tjji 
and ■02- 

The reduction procedure apphes also to differential equations of 
other types. As an illustration we consider the simplest problem of 
thermal conductivity: 

dtu - dlu2 = 0, (7.31) 
u{x,0) =Uo{x); u{0,t) =u{l,t) = 0. (7.32) 
Prom ip — d^u, equation (7.31) and conditions (7.32), we get 



u {x, t) 



X 1 

J (x-0-xJ{l-0 



(e, t) dc, 



u{x,t) = j ijj {x, rj) dfj + Uq (x) 



Accordingly, 



X 1 

/ {x-0-xJ{l-0 

.0 



'tp (C, t) d^- / 'fp [x, Tj) dr] = Uq {x) ; x,y e [0, 1] . 



In order to make an analogous reduction of the problem of bending 
of a rectangular plate of variable stiffness D, fixed along a contour [7], 



DAAu + 2d^Dd^Au + 2dyDdyAu + ADAu 
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- (1 - I/) {dlDdlu - 2d,yDd,yU + dlDdlu) = g, (7.33) 

a> (0, y) = 9> (a, y) = 9> (x, 0) = 9> (x, 6) = 0, n = 0, 1, 

(7.34) 

where A = 9^ + 9,^; z/ is the Poisson coefficient; q (x, y) is the intensity 
of the transverse load, we can set 

u{x,y)= / k {x, y, ^ ({, y)d^ + Y^ x^'^gij (y) . (7.35) 
J=i 

Here, 

9^A;i(x,y,x) =0, n = 0,l,2; 9^A;2 (x, y, y) ^ 0, x e [0, a] ; y e [0, 6] , 

and the functions gij (y) are intended to satisfy the conditions (7.34) 
for X = 0, X = a. The second representation of the solution via ip {x, y) 
is determined by means of the substitution of (7.35) into Eq. (7.33) 
and four-fold integration over the variable y. The appearing functions 
g2j {x) allow us to satisfy the conditions (7.34) for y = 0, y = b. After 
that, u {x,y) is eliminated from the representation of the solution.. 

Note that with the help of k {x, y, rj) one can easily satisfy conditions 
at isolated points inside the considered domain, e.g., u {xi, yi) — 0. The 
procedure of the reduction also applies to mixed boundary conditions 
(a change of the type along a side) and to the case of a connection 
of plates. Analogously, three-dimensional problems of mathematical 
physics can also be reduced to Predholm integral equations of the first 
kind. 

7.3 Universality and analogous approaches 

Thus, a comparatively elementary method of the reduction of linear 
boundary-value and initial-boundary-value problems to Fredholm in- 
tegral equations of the first kind is rather universal from the point of 
view of its realizations as far as the following aspects are concerned: 

- the order and structure of differential equations; 

- the form of boundary conditions; 

- the availability of variable coefficients; 
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- the form of the domain; 

- the dimensionahty of the problem. 

In this sitnation. all the information about a concrete problem is 
transferred into a functional equation, whose solution does not require 
any conditions on the contour of the domain, which poses a substantial 
advantage. Thus, its solution can be sought in the form of a series 
of a system of coordinate elements intended exclusively to ensure the 
efficiency of the procedure of the numerical realization. 

However, the problem obtained as a result of transformations is 
incorrect, hence its numerical realization requires adequate methods. 
At the same time, in applications, the solution of such a problem can 
be acceptably approximated by a series with the number of terms that 
does not affect the stability of the numerical algorithms. Therefore, 
one can hardly explain the absence of interest to a systematic use of 
this procedure, especially in the period before the general orientation 
at the discretization of problems of mathematical modeling. 

One may state that special literature did not point out the existence 
of a formalized method of the reduction of practically arbitrary initial- 
boundary- value problems to Fredholm integral equations of the first 
kind. At the same time, there a number of examples of applications 
of analogous transformations in rather particular situations. As a rule, 
they were a given physical interpretation that considerably disguised 
the generality of this approach. 

Thus, Yu. V. Repman used as a function closely related to bound- 
ary forces of a plate of a canonical configuration that allowed one to 
satisfy conditions on an internal contour of complex configuration [8]. 
L. A. Rozin has developed a method of separation that admits a re- 
duction of the problems of calculations of membranes to systems of 
Fredholm integral equations of the first kind for the forces of interac- 
tion of isolated bars ([9], section 9). Some publications point out the 
advantages of the approximation of higher-order derivative of differen- 
tial equations with respect to one of the variables that, compared to 
numerical differentiation, are much more accurate. However, as a rule, 
no comments were made on an actual transition to ill-posed problems 
(see, e.g., [10]). 

Some problems for differential equations, and, in particular, the 
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following one: 

d'xyU = adxU + bdyU + cu + f, 

where a, b, c and / are given functions of the variables x and y, can be 
reduced directly to Volterra and Fredholm integral equations of the sec- 
ond kind with respect to the higher-order derivative {ip — dxyu) . These 
issues are studied in detail by G. Miintz [11]. Of considerable interest 
is the fact, established by this author, that analogous transformations 
cannot be extended to the case of the simplest equation of the elliptic 
type. 



7.4 A connection to the algorithm of sec- 
tion 6.4 

The Fredholm integral equation of the first kind that arises as a re- 
sult of the reduction of two-dimensional boundary- value (initial- value) 
problems, can be represented in the form 

1 1 

/ n {x, y, V' it y) j T2 {x, y, 77) V' {x, rj)dri^ f (x, y) , x,y E [0, 1] , 


(7.36) 

where Ti{x,y,^), T2{x,y,r)) and f{x,y) are given functions; ip{x,y) 
has to be determined. 

Under the assunption that the function satisfying (7.36) exists and 
unique, it is represented in the form 



ip (x, y) = -00 {x, y) + -01 {x, y) , 



(7.37) 



where ipQ (x, y) and tpi (x, y) are solutions of the Fredholm integral equa- 
tion og the second kind 



V' {x, y) = /i 



1 

N {x, y, V' (C, y)d^ + j M {x, y, 77) V {x, rj) dr] 



1 1 



+ 



jd^j T{x,y,^,v,)i^{^,v,)d'n 







+ F(x,y), X, ye [0,1], (7.38) 
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with the free term, respectively, 







F {x, y) = F^{x,y)^xj H (x, A) k (C, y) d^; (7.39) 

-1 



F {x, y) = Fi {x, y) = 



1 



. (7.40) 



f{x,y) + jH{x,^,X)f{^,y)d^ 



In Eq. (7.38), the kernels are given by 

1 

N {x, y, = Ti {x, y,i) + xj H (x, C, A) n (C, 2/, dC,] 



M {x, y, T)) = T2 {x, y,r]); T {x, y, v) = A// {x, A) T2 (C, y, v) ■ 

(7.41) 

Here and above. 



^ 00 n 

H {x, e, A) = ^— ^ + 2 X: 3^:^ cos [2n7r (x - 0] , (7.42) 



where the parameter is < r < 1; 



K {x, y)= p{x,y)+A j L (x, A) p (C, v) d^. (7.43) 

-1 

In this expansion, A = A^; 

1 

p (x, y)^-xjh (x, V'l (e, y) ^^e, (7.44) 


where 

1 -r2 

/i(x,0 



1 — 2r cos [2n7r (a; — ^)] + 



= 1 + 2 ^ r" cos [2n7r {x - 0] ; (7.45) 

n=l 
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j,2n 



L (x, ^, A) = — + 2 y — — - cos [2nn (x - 01 ■ 

^ ^ 1-2A-A 1 - 2Ar" - Ar2" ^ ^ 

(7.46) 

The parameter 

A ^: 0; r"; ^r""; (-1 ± v^) r'^, n = 1, 2, . . . ; 
the parameter 

H ^ Hn, n= 1,2,..., 
where are the characteristic numbers of the homogenous equation 

1 1 



LO 
1 1 



J N(x,y,^)^lj{^,y)d^ + I M(x,y,ri)^lj(x,ri)dri 



+ 





x,ye [0, 1] 



After the choice of the values of fi, A and r has been made (note that 
they can be corrected afterwards, depending on different situayions), 
the sequence of the computational procedures is as follows: 

- determination of the kernels of Eq. (7.38), A'^, M and T, from 
(7.41), with the use of expression (7.42); 

- determination of the free term Fi from (7.40); 

- determination of the function ipi from (7.38), with F = Fi, 

- determination of the function p from (7.44), with the use of ex- 
pression (7.45); 

- determination of the function k from (7.43), with the use of ex- 
pression (7.46); 

- determination of the free term Fq from (7.39); 

- determination of the function ipQ from (7.38), with F = Fq; 

- determination of the sought function ip from (7.37). 

As can be seen, the algorithm of section 6.4 applies to the solution 
of the two-dimensional Fredholm integral equation of the first kind 
without any substantial changes. In this case, the variable y plays the 
role of a parameter. 
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7.5 A verification of the solvability of boundary- 
value problems 

In the above consideration, we have assumed that the function ip {x, y) 
satisfying the Fredholm integral equation of the firs kind (7.36) in the 
space L2 exists and is unique. Nonetheless, by formal use of the com- 
putational relations of section 7.4, one can "find" ip {x, y) also in those 
cases when Eq. (7.36) has no solution at all or has a variety of so- 
lutions. In the first case, the function ip {x, y) being substituted into 
(7.36) cannot satisfy this equation. 

Indeed, the function thus obtained is senseless because the construc- 
tion of the algorithm (see section 6.3) was based on the assumption 
that the function satisfying Eq. (5.1) existed, and what is more, the 
free term / {x) was interpreted as a result of previously performed inte- 
gration. However, on the other hand, if the function {x, y) found by 
means of the algorithm of section 7.4 does not satisfy Eq. (7.36) upon 
substitution, it implies that this equation is insolvable. 

So, what have we got in the end? The unpleasant properties of 
the Fredholm integral equations of the first kind, expounded on above, 
can be rather efficiently employed to verify the solvability of boundary- 
value (initiall-value) problems. Indeed, they are easily reduced to two- 
dimensional (or of higher dimension) Fredholm integral equations of 
the first kind, which was discussed in section 7.2. Consequently, after 
the realization of the algorithm of section 7.4, we are left only with a 
necessity to verify whether the obtained solution satisfies an equation 
of the type (7.36). 

For comparatively simple problems of the previous sections, such a 
verification is not very important; however, a lot of investigations are 
concerned with the adequacy of the problem (7.33), (7.34) with regard 
to the description of the bending u {x, y) at the corners of a rectangular 
plate. Of interest is another issue: one of the most important prob- 
lems of numerical simulations is, as a matter of fact, a formulation of 
problems that implies construction of differential or integro-differential 
equations. In this regard, the Fredholm integral equation of the first 
kind (after a reduction to it of a certain posed problem) may serve as 
a filter discarding invalid versions! 
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This short subsection seems to be important. Its brevity results 
from the fact that it is based on the material given above. 
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Chapter 8 



Other classes of problems 



8.1 The initial-boundary- value problem for 
the Korteweg-de Vries equation 

Let us assume that the problem 

dtu - 6ud^u + dlu = 0, (8.1) 

-u (x, 0) = Mo (x) ; u{Q,t) = ui{t)\ dxu{0,t) = U2{t)] ti (1, t) = M3 (t) , 

(8.2) 

has a unique solution in the space L2 for given functions Uq (x); Ui (t), 
i = 1,2,3. 

There is no general theory that would allow us to make a priori 
judgements about the solvability of the problems of this type. Results 
of numerical simulations as well as solutions of specially simplified equa- 
tions near the boundary (see [1], section 10) may prove to be the main 
tool of refinement on physical models.^ 

Using the procedure of the previous section, we can reduce the prob- 
lem (8.1), (8.2) to an integral equation of the first kind with respect 
to 

il) (x, t) = dlu {x, t) , 
^In this regard, the arguments of section 7.5 may prove to be useful. 
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which yields 

o,t) = ^J{x- 0' ^ (e, t) + ^x^gs (t) + xg2 (t) + gi (t) , (8.3) 



U [X, 



with the functions determined from the boundary conditions: 

gi (t) = ui (t) ; g2 (t) = U2 (t) ; 

1 

(73 (t) = 2 [us {t) - U2 {t) - ui {t)] - j {l-if^ (e, t) di. 



Substitution into (8.3) leads to the expression 

X 1 



U (x, t) 



^ (e, t) di 



.0 
+x'^U3 (t) +X{1-X)U2 (t) + (l - x"^) Ui (t) . 

Now we rewrite Eq. (8.1) in the form 

dtu — 6udxU — d^u. 



(8.4) 



(8.5) 



The substitution of (8.4) into the right-hand side of (8.5) and in- 
tegration from to i under the initial condition (8.2) allows us to 
determine 



u {x,t) — 6 / 6u (x, Tj) dxU (x, rf) — d^u (x, rj) drj + Uq. 



(8.6) 



The elimination of u {x,t) from (8.4), (8.6) leads to an equation of 
the form 



{A^l)){x,t) = f {x,t) , x,ten:0<x,t<l, 



where A is a nonlinear integral operator, and the function / depends 
on the data of the problem. 

In order to determine the function tjj [x, t), we can employ the algo- 
rithm of section 7.4 (the variable y is replaced by t). This function will 
satisfy a nonlinear integral equation of the second kind. By the con- 
traction mapping theorem, for small absolute values of the parameter 
H, its solution can be found by means of simple iterations [2]. 



8.2. A BOUNDARY-VALUE PROBLEM FOR A SUBSTANTIALLY NONLINEAR DIFFERENTIAL 



8.2 A boundary- value problem for a sub- 
stantially nonlinear differential equa- 
tion 

Here, we discuss nonlinearity related to higher-order derivatives. As an 
example, consider Monge-Ampere's equation : 

dludyU — {dxyuf = sxdlu + S2dyU + s^d^yU + q, (8.8) 

where Sj, i = 1, 2, 3 and g, in general, depend on the variables x, the 
sought function u ( x, y) and its first derivatives d^u, dyU [3] . 
Let us assume that Sj — Si{x,y), q — q {x, y) and 



u (0, y) = u (1, y) = u (a;, 0) = -u (x, 1) = 0. 



L9) 



We also assume that the solution of the problem in L2 exists and is 
uniqe. Using the notation 

dlu (x, y) = -01 {x, y) ; (x, y) = -02 (a^, y) ; 

d^yu{x,y) = ipix^y) , x,y e Q : < x,y < 1, 
taking into account (8.9), we get 



u {x, y) 



X i 

J (x-O-xjil-O 



i^i (C, y) d^; 



u{x,y) 



y 



1 



{y -v) - y J C^-v) 


X y 



i)2 (x, 7]) di]; 



u{x,y) ^ J J ^ (C, rj) dr]. 



Upon the substitution of these expressions into (8.8) and the elimi- 
nation of the function u, we reduce the problem to the following system 
of equations: 



-01 {x, y) -02 {x, y) - V'^ (x, y) = si {x, y) ipi {x, y) 
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+S2 (x, y) -02 {x, y) + S3 {.X, y) -0 (x, y) + q (x, y) ; (8.10) 

' X 1 "I X y 

/(x-0-^/(l-0 ^i{U)di- j dij ^{i,'q)drj^Q- (8.11) 







X y 







(8.12) 

Two-fold differentiation of Eqs. (8.11), (8.12) with respect to x and 
y yields, respectively. 



tpi {x, y) = j 9x'4' {x, v) dm tp2 {x, y) = j ^yi) (C, y) d^. 



Equation (8.10) takes the form 

' y 1 r ^ 1 ^ 

j d^ip {x, Tj) drj j dy^ (^, y) d^ --0^ {x, y) = Si (x, y) j d^jp {x, rj) dr] 
.0 J Lo J 

X 

+S2 {x, y) J dyi/j (^, y) d^ + sg {x, y) i/j (x, y) + q (x, y) , x,y efl, 



and after integration in the limits 0, x and 0, y reduces to the following: 



{Alp) {x, y) = f {x, y), x,y efl, 
where is a corresponding nonlinear operator; 

X y 

f {.x,y) = j di j q (C, Tj) drj. 



(8.13) 







The above implies the boundedness of the derivatives dxSi^ dyS2- A 
possible way of the solution of this equation is discussed in section 8.1. 



8.3. NONLINEARITY OF THE BOUNDARY CONDITION 
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8.3 Nonlinearity of the boundary condi- 
tion 

Consider a typical problem of the irradiation of an infinite plate with a 
thermally insulated surface into a medium whose absolute temperature 
is equal to zero [4] : 

dtu - adlu = 0, (8.14) 
u{x,0) =uo{x); pd^u{0,t)+u"'{0,t) = 0; d^u{l,t) = 0. (8.15) 

Here, u{x,t) is the temperature gradient; Uo{x) is a given func- 
tion; a is the temperature conductivity; p = A/a, with A, a being the 
thermal-conductivity and the heat-transfer coefficients, respectively; m 
is a parameter. 

Introduce the notation 



9lu {x, y) = ip {x, t) , 



^.16) 



which leads to 



u 



X 

= / {x-^)ij{^,t)d^ + xgi {t)+g2{t), 



where gi (t) are functions of integration. 
The boundary conditions (8.15) yield 

1 

g,{t)^-Ji; (e, t) dC; pgi it) + gT (t) = 0, 



and, accordingly, 

X 11 r 1 



u (x, t) 



X — ^) — X 





p 



Using (8.14), (8.16) and taking into account the initial condition 
^.15), we get 



u{x,t) — a J ip {x, rj) dr] + uq (x) , 
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and the problem is reduced to the solution of the nonlinear integral 
equation of the first kind (8.7), where 



A. 



11 



j{x-^)-xj d^+ipj.dn -aj ,dr); 



f {x, t) = uo (x) . 



8.4 A small parameter by the highest-order 
derivative of the differential equation 
of the problem 

As an illustration of general considerations, we consider the problem 
of heat transport induced by the processes of thermal conduction and 
convection (the first and the second terms of the equation, respectively) 
[5]: 



dfU = ed^u + (3dxU. 
Here, /? > is a constant; e is a small parameter, 

u{x,0)=Q; u{l,t) = 0; u {0,t) = Ui (t) , 



1.17) 



with Ui (t) being a given L2-function. 

The notation (8.16) under the boundary conditions (8.18) leads to 



u {x, t) 



X 1 

J {x-0-xJ{l-0 

.0 



(8.19) 



The integration of (8.17) in the limits 0, t with the use of (8.19) and 
of the initial condition (8.18) yields 



u {x, t) 



X 

# (x,77) + /3 j i){i,r})di 



drj + ui (t) 



8.4. A SMALL PARAMETER BY THE HIGHEST-ORDER DERIVATIVE OF THE DIFFERENTIAl 



The elimination of the function u from these relations leads to Eq. 
^.7), where 

A = A' + A", 



Z XIX 1 

A: = eJ.dv; A: = f3jd^J.dr]- j {x - i) - x j {I - i) 



LO 

/(x,t) = -«i {t) 



The algorithm of section 7.4 allows us to reduce the problem to a 
Predholm integral equation of the second kind of the form 

s {x, t)^n [{eRi + R2) s] (x, + / (x, t) , (8.20) 

where Ri and R2 are corresponding integral operators; s {x, t) = ip {x, t). 
As a result of the expansion [6] 

00 

m=0 

we get a sequence of recursion relations 

So {x, t)= II {R2S0) {x, t)+ f {x, t) ; 
si {x, t)^ II (-R2S1) {x, t)+n {RiSq) {x, t) ; 

Sm+l {X, t)^ H {R2Sm+l) {x, t) + li (RiS^) {x, t) 

that are canonical Predholm integral equations of the second kind. 

It follows from the above that the proposed approach is rather effi- 
cient in the problems of mathematical physics with a singular pertur- 
bation, whose numerical realization, as a rule, meets with considerable 
difficulties (see, in particular, [7]). Indeed, we managed to transform 
the singular perturbation (8.17) into the regular one (8.20), which fa- 
cilitated a radical simplification of the problem.^ 



^Singular and regular perturbations affect, respectively, main and dependent 
terms of the operators. 
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8.5 Equations of a mixed type 



Boundary- value problems for equations of this type are characterized by 
complexity of the investigation into the issues of existence and unique- 
ness (see [8]). As a consequence, one has to consider such equations 
on rather special domains, which restricts the field of practical appli- 
cations. 

Leaving this issue be, only for the sake of an illustration of the 
procedure of reduction, we turn to well-known Tricomi's equation 



ydlu + dyU = 



(8.21) 



that belongs both to the hyperbolic and elliptical types for y < and 
y > 0, respectively. As an for example, we employ the following bound- 
ary conditions: 



u (0, y)=u (1, y) = u {x, -1) = 0; u {x, 1) = u {x) , 

where the function u (x) is such that (0) —1/(1) — 0. 
Prom the notation 



by (8.22), it follows: 



u{x,y) 



dlu (x, y) = -0 {x, y) , 



X 1 



(8.22) 



(8.23) 



Two-fold integration of Eq. (8.21) in the limits —1, y under the 
conditions (8.23) and (8.22) yields the expression 



u{x,y) = 



y ^ _^ 1 



L-1 



■qil) {x,r])dr]+-{l + y)u{x) . 



The problem reducea to a Predholm integral equation of the second 
kind (8.13) on the domain Q:0<a;<l,— l<y<l with the operator 



A. 



jix-O-xjii-O l{y-v)-^la-v) 



r],dr] 



8. 6. THE INVERSE PROBLEM OF THE RESTORATION OF THE COEFFICIENT OF THE DIFFI 



and the free term ^ 

fix,y) = 2 (l + ^)^(^)■ 
Notethat the so-called condition of "matching" on the line of parabolic 
degeneracy y = 0, imposed on the solution of Eq. (8.21) ([8], p. 27), is 
fulfilled in a natural way: 

y}^o ^ ^ y}^o ^ ^) ' ^ ^ ^] ' 

lim^ dyU {x, y) = hm^ dyU (,x, y) , x E [0, 1] . 

As in the previous subsection, this situation results from the fact 
that the singularity of the problem is transferred from the main term 
of the relevant operator to the dependent one. 



8.6 The inverse problem of the restora- 
tion of the coefficient of the differen- 
tial equation 

Small oscillations in the transverse direction of a stretched string of 
variable density are described by the equation 

d'^u = a (x) dlu. (8.24) 

Here, x, t are dimensionless coordinates; 

a(x) ^ NT^/p{x) 

with N being the tension, p (x) the density of the material, 21 the length 
of the string, T the time interval. 

We assume that the ends of the string are fixed, whereas its density 
and the oscillations arc symmetric with respect to the coordinate x — 0. 
The corresponding boundary conditions have the form 

d^u{0,t)^u{l,t)^0. (8.25) 

We also employ the following initial conditions: 

u {x, 0) = -uo {x) ; dtu {x, 0) = 0. (8.26) 
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The coefficient a (x) is to be determined from (8.24)-(8.26) for given 
Uo{x), N, I, T and additionai information on tlie oscillations of the 
middle cross-section of the string: 

u(0,t)^u{t). (8.27) 

We assume that the conditions ensuring the existence and uniqueness 
of the solution of the considered problem ([9], section 4) arc fulfilled. 

By analogy with what was done many times before, using the no- 
tation (8.16) and (8.24)-(8.26), we find 



u {x, y) 



X 1 

/(^-O- / (1-0 

.0 



^P{U)di; (8.28) 



z 

u {x, y) = a (x) j {t — r])ijj {x, rj) dr] + uq (x) . 



By eliminating u {x, t), we obtain an equation of the type (8.7). The 
substitution of (8.28) into (8.27) leads to the integral equation 

(^V) it) ^f{t), te [0, 1] , 

where 

1 



The procedure of the so posed system of equations can be viewed in 
the context of supplementing the algorithm of section 7.4 by iterations 
with the function a. 



8.7 The problem of the Stefan type 

Consider the classical model [10]: 

dtu = dlu, < X < 7 (t) ; < i < 1, (8.29) 



u{x,Q) =Uq{x)] u{0,t) = u{j{t),t)=0, 'Uo(0) = 0. 



(8.30) 



8. 7. THE PROBLEM OF THE STEFAN TYPE 
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On the moving boundary that separates the phases an additional 
condition is imposed: 



ad,u(^(t),t)^^'(t), 7(0) =70, 



(8.31) 



where 70 > 0; the constant a can be both positive and negative; 7' (t) = 
^7 (t) /dt. 

Thus, the data of the problem are Uq {x), a and 70; the functions 
u {x, t) and 7 (t) are to be determined. 

In Eqs. (8.29)- (8.31), we make a non-orthogonal mapping 



X — x/^ (t) , t — t 

on a canonical domain fl : < x,i < 1. We get: 

dtu - [x-f' (t) /7 (t)] d^u = d^u, 

u{x,0) = 0; u{0,t) =u{l,i) = 0] 
ad^u{l,t)^i{t}, 7(0) = 7o- 
By analogy with the above, the notation 

d^u {x, t)=i) (x, i) , 

conditions (8.34) and equation (8.33) lead to 



(8.32) 

(8.33) 

(8.34) 
(8.35) 



u 



(x,t} 



X J. 

j{x-i)-xj{l-i) 



^lj{^,t}d^; (8.36) 



X 1 



U [X 



,t)^jU{x,r)) + [xj'{t)/j{t)] J-J{l-0 



LO 



The substitution of (8.36) into (8.35) yields 

1 

7' {t}^al (e, I) di, 



V' (C, "n) di i d-q 
i.37) 
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from which we get 



1 t 



(8.38) 





Accordingly, in the expression (8.37), we have 



xV (t) 



1 It 
«^ / ^V' it t)dC/ a j idi j ^ (^, 77) dr] + 70 



The ehmination of u{x,t) from (8.36), (8.37) leads to the integral 
equation of the first kind (8.7). The function t/S determined from this 
equation should be approximated by an analytical dependence on x 
in order to make an inverse change of variables. The sought separa- 
tion boundary 7 it) is determined from the nonlinear integral equation 
(8.38). Then, by (8.36), using (8.32), we can calculate the function 
u (x, t). 
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Chapter 9 
Conclusions 



Let us summarize the main points of the above consideration. Thus, 
the solution of the Fredholm integral equation of the first kind 

1 

(Aij) {x) = Jk {x, iO di = f{x), xe [0, 1] (9.1) 



in the " convenient" for the numerical realization space L2 is an ill-posed 
problem. In the case of the space Z2 that is adequate to the range of 
the operator A, the situation is different: the data of Eq. (9.1) may, 
theoretically, satisfy the conditions of its correctness, but, nevertheless, 
the solution will constitute a series that diverges as a result of the 
accumulation of errors of calculations. 

It should be noted that even a verification of whether the data of 
(9.1) belongs to the space I2 is, in general, infeasible. At the same 
time, as an objective factor of incorrectness, there appear the error of 
experimental determination of / (x) and, sometimes, inaccurate infor- 
mation about the function k {x, ^) that characterizes the system under 
consideration. 

The basis of our work is formed by the suggestion to connect adap- 
tively Eq. (9.1) to the space ^2 by means of a modelling of the er- 
ror (Sf) {x) that arises owing to the smoothing of information by the 
procedure of integration. We assume that the function satisfying this 
equation exists, is unique, and the condition 

{Sf){x)=Q, xe[0,l], (9.2) 
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reflecting objective smallness of tfiis error compared to xjj {x) and the 
data of tlie problem, is employed. 

Starting by heuristic considerations that were later supported by 
more firm arguments, we demonstrated the expediency of the represen- 
tation of the error as a difference between the sought function in an 
explicit form and of the integral component: 

{df){x)=ij{x)-\{Bij){x), xe[0,l], (9.3) 



where 



i){x), xe [0, 1] ; 
ip{x), a;e[-l,0). 



Of exclusive importance for the whole complex of the transforma- 
tions, especially in the first version of their realization, was the repre- 
sentation 

oo 

/i (x, ^) = 1 + 2 ^ r" cos [2n7r (a; - 0] , < r < 1, (9.4) 

n=l 

i.e., in the form of Poisson's kernel, that allowed us to satisfy the con- 
dition (9.2) for the case when the function ip (x) was harmonic. 

By the use of (9.2), (9.3), the formulation of the problem (9.1) was 
transformed: 



where 



(x) 



i^ix), xe [0, 1] ; 
xe[-i,o), 



*W = A(M)M + { [0,1]; 
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(9.8) 

Prom the point of view of constructiveness of what followed on the 
basis of (9.2)-(9.4), a key role was played by the following factors: 

- an extension of Eq. (9.3) under the condition (9.2) to a; G [—1,0). 
As a result, there arose the Fredholm integral equation of the second 
kind (9.5) with an undefined function k,{x), a component of its free 
term; 

- the use of the equation of analogous structure, Eq. (9.6), whose 
the free term goes to zero on the second part of the interval of definition; 

- an incompletely continuous perturbation of the operator A by 
consecutive addition of (9.5) and (9.6) to (9.1). As a result, Eqs. (9.7) 
and (9.8) arose. 

The elimination of the function x i^) from (9.8) using the equation 
*; (.) = A (B%) (.) + {" W \^^fj^ ^ '1 • (9.9) 
where 

in' _j i^o{x), xe [0, 1] ; 
^"^^ - \ (^'o (x) , xe[-l,0), 

allowed us to obtain 

*;(x) = A(B*^W + { (9.10) 
Here, 

yT.'(^\-j i^i (x) = ij {x) - {x) , a;e[0,l]; , . 

^'^''^-\^[{x) = ip'{x)-ipUx), xe[-i,0). ^^-''^ 

The inversion of the operator I ~\B in (9.10) leads to the Fredholm 
integral equation of the second kind 

1 

^|J^{x)^^lJK{x,0^liOdC + fl{x), xe[0,l], (9.12) 
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where the kernel and the free term are determined by the data of (9.1); 
the parameter /i , as in analogous cases A, must satisfy a solvability 
condition. After the determination of ipi {x) , the function ip[ {x) is 
given by quadratures. 

Prom (9.5) and (9.6), it follows: 



A further orientation of the transformations was concentrated on the 
determination of the difference in the right-hand side of this equation 
in order to turn it into a Predholm integral equation of the second kind 
with respect to x{^)- The fact that such a result could be achieved 
was by no means obvious, because, although the function (/?' (x) was 
expressed in a simple way via x {^) [from Eq. (9.6)], such a possibility 
was absent for the function (x). 

In the process of attaining the set objective, a stress was put on 
obtaining a relation between the sohitions of (9.7) and (9.8) that would 
allow one of this equation to turn into another. The form of the free 
terms of Eqs. (9.7) and (9.8) implies a possibility of such transfor- 
mations, that is, a possibility of a "flow" of their nonzero components 
from one part of the interval of definition to the other, which opens 
up a prospect of the representation of ip {x) in two ways, i.e., with and 
without the function x (^) in s-rL explicit form. 

The actual realization of the above arguments showed that the so- 
lutions of Eqs. (9.5)-(9.8) on x e [—1,0) are mutually related via the 
functions entering (9.9): 








(9.13) 



-1 



ip (x) = ip' (x) + ipQ (x) ; 



(9.14) 



1 




(9.15) 







Prom (9.11), (9.14), 



if (x) - if' (x) = (p'q (x) = if' (x) - Lp\ {x) , 



(9.16) 



197 



and, after substitution into (9.13), the problem reduced to the solution 
of the Fredholm integral equation of the second kind 

1 

x{x)=Ajl{x,OxiOdi + q{x), xe[0,l]. (9.17) 



Here, A = A^; 



l(x,0^ I h(x,OH(C,^,X)dC, 
-1 

with H (x, ^, A) being the resolvent of the kernel h (x, ^) on x E [—1, 1]; 



q{x) = -X J h{x,0^'i{Odt 
-1 

The solution to Eq. (9.17) has the form 

1 

X{x) = q{x)+A J L(x,e,A)g(0 d^, 



where 

L (x, A) = ^ + 2 y ^ — - cos [2n7r (x - C)] 

^ ' ^' ^ 1 - 2A - A ^1 1 - 2Ar" - Ar^" ^ ^ ^'^ 

is the resolvent of the kernel / {x, ^). 

The inversion of the operator I — XB in Eq. (9.5) allowed us to 
represent the function as a Fourier series in terms of the elements 
{cos {2mTx) , sin {2mTx)} whose coefficients were expressed via the data 
of (9.1) and depended on the parameter r. Note that on the previous 
stage of the calculations by (9.12) the function ipi (x) was determined 
in an analogous form. 

The solution so obtained was restricted only by the case when the 
function ip (x) satisfying Eq. (9.1) was harmonic. However, 
to the limit r — > 1 easily removes the problems by transferring the free 
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term of Eq. (9.1), as well as the function ip {x) satisfying this equation, 
into the space L2. 

This is the main point in achieving the final objective of the trans- 
formations. Accordingly, condition (9.2) takes the form 

II^/IIl.(o,i) = 0. 

In general, the transformations seem to be rather transparent. Thus, 
the determination of the function tp (x) is transformed into the prob- 
lem (9.5)-(9.8). The obtained function ^' (x) is a part of the solution of 
Eq. (9.8) that depends on -/if (x). From (9.14) and (9.11) the function 
(P'q (x) is determined in two ways, which is reflected by relation (9.16). 
Hence (9.13) turns into Eq. (9.17). Finally, in the obtained solution a 
passage to the limit with respect to r was made. 

The following interpretation of the algorithm of the reduction is pos- 
sible. First, the transformed formulation of the problem is "deformed" 
by eliminating the function x i^) from Eq. (9.8). Then this "deforma- 
tion" is adaptivcly smoothed out by means, which is very important, 
of the solution of the Fredholm integral equations of the second kind 
(9.12), (9.17) and (9.6). 

We have discussed the first version of the solution of the problem. 
The second version of its solution is also based on the relations given 
above. By the use of (9.15), an analogue of Eq. (9.17) for the function 
K {x) was obtained. Its solution has the form 



K{x)=q'{x)+AjL{x,^,A)q'{Od^, ^e[-l,0), (9.18) 
-1 

where 

1 

q'{x)^-X J h{x,0i^i{0d^- (9-19) 



The problem (9.1) is reduced to a numarical realization of, consecu- 
tively, two Fredholm integral equations of the second kind, namely, Eq. 
(9.12) and 

1 

V'o (^) = / {x, ^0 (0 + Fo{x), xe [0, 1] , (9.20) 
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where 

1 



As a result, its solution is sought in the form 

ip (x) = ipo (x) + ipi (x) , 

see also (9.9). 

Note that the kernels of these equations are the same. It is also 
shown that, by the use of (9.5), (9.11), (9.18) and (9.19), the problem is 
reduced to a numerical realization of a single Frcdhom integral equation 
of the second kind with respect tot he function ip (x). 

In contrast to the previous version of the solution, there is no need 
here to evaluate the Fourier coefficients of the functions k {x, {) and 
/ (x) and perfor the summation of infinite series, which may be regarded 
as an advantage. At the same time, universal algorithms are available 
for the solution of Eqs. (9.12) and (9.20). In general, the second version 
of the solution of the problem is more formalized. As an advantage of 
the first version, one should point out a possibility of obtaining the 
function ip (x) in a convenient, as a rule, form of a Fourier series. 

The principal difference between the two versions lies in the way of 
satisfying (9.2), or the equation 

1 

\jh{x,i)^{i)di^i^{x)-xjh{x,i)i^{i)di, xe[0,l], (9.21) 

-1 

where, for ip {x) G L2 (0, 1), the function ^p{x) can only be a general- 
ized function. In contrast to the first version, where, for this reason, 
the transformations were performed with the function ip (x) that was 
assumed harmonic up to the final stages, in the second version, it was 
imphed that Eq. (9.21) was satisfied in the sense of generahzed func- 
tions. 

Specifically, we employed an equation obtained by applying to (9.21) 
the operator 

1 

jh{x,i),di 
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with respect to to the generahzed functions^ 

1 
-1 

At the same time, exactly condition (9.2) appears to be absolutely 
necessary for the realization of both the first and the second versions 
of the solution of the problem. Indeed, equation (9.21) that may be 
called "free-lance" cardinally changes the problem (9.1) with regard to 
the solvability of the Predholm integral equation of the first kind. With 
the help of (9.2), one essetially removes an inherently insurmountable 
problem of an objective mismatch of f {x) and R{A), which is the 
reason for the problem (9.1) being ill-poscd. 

The liberation of / (,x) from formal association with R [A) by means 
of (9.2) and (9.7), (9.8) simultaneously results in the fact that the free 
term of Eq. (9.21), when considered as 



where 



{B'^){x)=f{x), a:e[0,l]. 



1 

: = lh{x,0.d^, f{x)^J^^;{x)- J h (x,^) ^^; (0 



-1 

becomes functional. 

As a consequence, the condition f {x) e R {A) (that is actually 
infeasible) is replaced by the following: 

~f{x)eR{B'), (9.22) 

which, in fact, is equivalent to 

/ (x) + (5/) /^^eR {A) . 



^Note the following characteristic feature: the final result of the transformations 
appears to be the same as if, without these transformations, one postulated the 
applicability of the theory of Fredholm integral equations of the second kind to the 
case when the function ^ (x) from Eq. (9.5) is a generalized function. 
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Thus, it proves to be possible to go over from a numerical compari- 
son between / (x) and R (A) just to the question of the existence of 
the function function ip (x) allowing for the fulfillment of the condition 
(9.22). 

Moreover, given that (9-21) is a Fredholm integral equation of the 
second kind with respect to i/j {x), in the course of subsequent transfor- 
mations, there occurs, in a sense, a readdressing of the status between 
/ (x) and R {B'). Namely, the range of the operator B' manifests itself 
as the free term, and the problem essentially reduces to finding the 
function ^0 [x) from it. Here, the fact that R {A) is not closed does not 
play any role. 

Note the following: as a result of (9.21), the determination of the 
function %l)[x) satisfying Eq. (9.1) was carried out, figuratively, by 
" materialized pressing" with regard to the validity of seemingly abstract 
Banach's theorem on the inverse operator. Specifically, this is done by 
the identity operator from I — XB by ensuring the entering of the above- 
mentioned function in an explicit form. 

It is shown that wide classes of problems of numerical simulation are 
easily reduced to Fredholm integral equations of the first kind. After 
that, the procedure of correct formulation and of constructive realiza- 
tion, discussed for a one-dimensional case, is directly extended to them. 
Therefore, a differentiation between direct and inverse formulations of 
problems of mathematical physics to a certain extent loses significance. 
We have also proposed a method of verification of the solvability of 
problems formulated in terms of partial differential equations. 

In light of the above, we can draw a conclusion that, if the phe- 
nomenon (process) admits an adequate description by methods of nu- 
merical simulations, the restoration of its underlying cause or of dif- 
ferent parameters from an objectively sufficient volume of additional 
information does not pose principal difficulties, because the correspond- 
ing problems can be well-posed. From this point of view, an analysis 
of actually observed events, including multi-factor social-economic and 
ecological processes, can be done with much larger efficiency. 

Maybe, it would be reasonable to suggest that, in general, the pro- 
cess of the understanding of the World is much simpler than a wide 
audience usually supposes it to be under the influence of the sphere of 
applied science that, at present, armed with means of electronic pro- 
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cessing of information, constitutes, in fact, a natural monopoly with an 
almost dominant role of commercial component and, correspondingly, 
a systematic drive for investment? 

Thus, colossal means are invested in problems of the restoration of 
dependencies from empirical data and, in particular, in remote probing 
of the surface of the Earth by spacecraft. What is actually realized is a 
search for minimally and maximally acceptable values of the parameter 
a in the integral equation of the type 

1 

aij{x) + J k{x,0i^{0d^ = fix), xe[0,l]. 



The essence lies in the necessity to establish a balance between 
computational and, respectively, financial abilities of the solution of 
an almost degenerate algebraic problem and an approximation to the 
"exact" formulation that is associated with the factor of incorrectness 
for a = 0. 

In this regard, we note that, of course, it would be incorrect to 
suppose that problems in science are altogether absent or that one can 
develop, irrespective of the circumstances, efficient means to overcome 
these problems. However, in our opinion, complications of principal 
character are inherent, in the first place, to direct formulations of some 
problems, that is, to the construction of mathematical models of insuf- 
ficiently studied processes and phenomena. 

It is clear that the solution of some classes of inverse problems of 
numerical simulation may also pose substantial difficulties, but, never- 
theless, the wide-spread dogma that the procedure of the restoration 
of the cause from the consequence is ill-posed, in general, seems to be 
manifestly erroneous. 

J. Hadamard's statement that the problems that adequately de- 
scribe real processes are well-poscd is an ingenious idea, whose con- 
structive development allows one to attain a qualitatively higher level 
of the potential of methods of numerical simulations. 
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